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Abstract 

This  report  offers  details  regarding  the  theory  and  mathematical 
formulations  implemented  in  the  newly  developed  general  constituent 
simulation  module  (GC),  contaminant  simulation  module  (CSM),  and 
mercury  simulation  module  (HgSM).  The  GC  models  simple  kinetics  for 
multiple  size  classes  of  solids  and  user-defined  constituents  in  the  water 
column.  The  CSM  models  the  kinetics  of  user-defined  contaminants  in 
aquatic  systems.  The  kinetic  processes  modeled  in  CSM  include:  ionization, 
multi-phase  partitioning,  degradation,  photolysis,  hydrolysis,  volatilization, 
generalized  second-order  reaction,  and  transformations  where  one  chemical 
species  undergoes  a  reaction  and  is  transformed  to  a  daughter  product.  The 
HgSM  models  mercury  species  (elemental  mercury,  inorganic  mercury,  and 
methylmercury)  and  their  cycling  in  aquatic  systems.  The  three  water 
quality  modules  only  compute  internal  sources  and  sinks  of  each  state 
variable  for  both  water  column  and  an  underlying  sediment  layer.  The  GC, 
CSM  and  HgSM  modules  are  written  as  “plug  in”  dynamic  link  libraries  and 
compiled  as  GC.dll,  CSM. dll  and  HgSM.dll  respectively.  These  modules 
have  been  integrated  into  the  HEC-RAS  model  and  will  also  be  integrated 
into  other  hydrologic  and  hydraulic  models  (e.g.  AdH  and  SRH-2D). 
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1  Introduction 

1.1  Background 

Contaminants  are  introduced  into  aquatic  systems  as  a  result  of  the 
weathering  of  soils  and  rocks,  and  from  a  variety  of  human  activities 
involving  the  mining,  processing,  or  use  of  substances  that  contain  toxics. 
There  are  different  types  of  sources  of  contaminants.  The  first  type  are  point 
sources,  where  contaminants  come  from  a  single  identifiable  source.  The 
second  type  of  contaminant  sources  are  non-point  sources,  where 
contaminants  come  from  distributed  sources.  Once  dissolved  in  water, 
contaminant  concentrations  may  exceed  water  quality  standards  and  impair 
the  designated  water  uses.  Exposure  of  aquatic  organisms  to  contaminants 
can  cause  toxicity,  or  result  in  bioconcentration  through  the  food  chain. 
Human  health  may  also  be  affected  by  ingesting  contaminated  water  or  fish. 
As  a  result,  the  criteria  for  protecting  human  health  and  indigenous  aquatic 
communities  have  been  promulgated  by  the  U.S.  Environmental  Protection 
Agency  (EPA)  for  toxic  chemicals  (Ambrose  et  al.  1993).  Of  these  toxics, 
some  organic  chemicals  and  metals  are  of  particular  concern  within  the 
aquatic  environment  due  to  their  high  toxicity  and  having  several  known 
adverse  effects  on  the  natural  ecosystem  and  human  health  (Wool  et  al. 
2006;  Ullrich  et  al.  2001;  Brigham  et  al.  2009). 

As  the  scientific  and  regulatory  communities  realize  the  significant 
environmental  impacts  of  contaminants,  it  is  increasingly  imperative  to 
apply  modeling  tools  to  evaluate  and  predict  the  fate  and  transport  of  these 
contaminants  in  the  environment.  Hydrologic  and  hydraulic  (H&H)  models 
coupled  with  water  qaulity  become  indispensable  tools  used  by  environ¬ 
mental  analysts.  Over  the  last  three  decades,  a  variety  of  H&H  models  have 
been  developed  for  predictive  estimates  of  contaminant  concentration  and 
transport  pattern  (Bahadur  et  al.  2013).  However,  the  Hydrologic 
Engineering  Center- River  Analysis  System  (HEC-RAS)  (HEC  2010),  two- 
dimensional  hydrodynamics  and  water  quality  (CE-QUAL-W2)  model  (Cole 
and  Wells  2011),  and  Adaptive  Hydraulics  (AdH)  (Berger  et  al.  2012), 
continuously  developed  and  maintained  by  the  U.S.  Army  Corps  of 
Engineers  (USACE),  are  still  lacking  the  capability  of  modeling  contami¬ 
nants  in  aquatic  systems.  Recognizing  this  need,  the  research  and 
development  of  contaminant  and  mercury  simulation  modules  for  a  variety 
of  H&H  models  were  sponsored  under  the  USACE  Environmental  Quality 
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Technology  (EQT)  research  program.  These  modules  were  developed  to  add 
additional  water  quality  capabilities  for  existing  H&H  models. 

1.2  “Plug  In”  water  quality  modules 

A  set  of  “plug  in”  nutrient  simulation  modules  (NSMs)  for  existing  H&H 
models  have  been  developed  at  the  ERDC  Environmental  Laboratory 
(Zhang  and  Johnson  2016).  Following  the  same  principles  and  protocols, 
the  contaminant  simulation  module  (CSM),  mercury  simulation  module 
(HgSM)  and  general  constituent  simulation  module  (GC)  were  also 
developed  as  “plug  in”  water  quality  modules  for  existing  H&H  models. 
The  GC,  CSM,  and  HgSM  modules  were  all  written  as  dynamic  link 
libraries  (DLLs)  and  they  must  be  integrated  into  H&H  models  to  predict 
concentrations  of  water  quality  constituents  in  an  aquatic  system.  The 
H&H  models  simulate  the  physical  processes  of  advection,  dispersion,  and 
boundary  inflow,  and  are  responsible  for  computing  the  transport  of 
constituents  across  the  model  domain.  The  GC,  CSM  and  HgSM  modules 
simulate  the  speciation,  reactions,  and  transformations  of  state  variables 
and  are  responsible  for  computing  internal  source  and  sink  terms  in  each 
water  quality  cell. 

As  with  the  NSM  modules,  newly  developed  GC,  CSM,  and  HgSM  modules 
were  designed  to  be  independent  of  the  dimensionality  of  the  spatial 
domain.  The  water  quality  computation  domain  only  represents  a  well- 
mixed  cell,  including  the  water  column  and  active  sediment  layer.  The 
sediment  layer  represents  the  shallow,  biologically  active  layer  of  surficial 
sediments.  The  depth  of  the  sediment  layer  should  be  chosen  based  on  the 
problem.  The  GC,  CSM,  and  HgSM  modules  compute  kinetics  of  water 
quality  constituents  occurring  in  the  water  column  and  sediment  layer. 

The  CSM  was  designed  for  simulating  general  contaminants  such  as 
pesticides,  polychlorinated  biphenyls,  halogenated  aliphatic  hydrocarbons, 
halogenated  ethers,  monocyclic  aromatics,  phthalate  esters,  polycyclic 
aromatic  hydrocarbons,  and  nitrosamines.  These  chemicals  can  cause 
toxicity  to  aquatic  organisms,  or  bioconcentrate  through  the  food  chain, 
affecting  human  health.  The  HgSM  was  designed  specifically  for 
simulating  mercury  cycling  and  kinetics  in  aquatic  systems.  Many  of  the 
concepts  and  kinetic  algorithms  employed  in  CSM  and  HgSM  modules  are 
adopted  from  the  U.S.  EPA  developed  Water  Quality  Analysis  Simulation 
Program  (WASP)  (Wool  et  al.  2006)  and  the  Electric  Power  Research 
Institute  (EPRI)  developed  Dynamic  Mercury  Cycling  Model  (D-MCM) 
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(EPRI  2013).  Besides  modeling  solids  and  user-defined  constituents,  the 
GC  module  was  developed  for  computing  time  varying  parameters  for 
CSM  and  HgSM. 

1.3  Report  outline 

This  report  describes  the  theory  and  mathematical  equations  that  are 
implemented  in  the  GC,  CSM,  and  HgSM  modules  and  consists  of  five 
chapters: 

•  Chapter  1  introduces  the  background  and  concept  of  “plug  in”  water 
quality  modules 

•  Chapter  2  describes  the  processes  and  mathematical  equations 
included  in  the  GC  module 

•  Chapter  3  describes  the  processes  and  mathematical  equations 
included  in  the  CSM  module 

•  Chapter  4  describes  the  processes  and  mathematical  equations 
included  in  the  HgSM  module 

•  Chapter  5  provides  the  summary. 

The  appendices  provide  the  user  with  the  information  necessary  to 
understand  the  symbols  used  in  each  module: 

•  Appendix  A  summarizes  mathematical  symbols  used  in  the  GC  module 

•  Appendix  B  summarizes  mathematical  symbols  used  in  the  CSM 
module 

•  Appendix  C  summarizes  mathematical  symbols  used  in  the  HgSM 
module. 
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2  General  Constituent  Simulation  Module 
(GC) 


Many  waterborne  organic  and  inorganic  chemicals  exist  both  in  dissolved 
form  and  as  forms  adsorbed  to  solids  and  other  sorbents.  Solids  are  impor¬ 
tant  as  an  adsorbing  surface  for  contaminants  and  play  an  important  role  in 
the  transport  of  contaminants  within  aquatic  systems.  Contaminants  in  the 
water  column  may  be  lost  or  gained  due  to  the  settling  and  re-suspension  of 
sediment  solids  onto  which  contaminants  are  adsorbed.  Simulation  of  solids 
mass  balance  is  critical  step  for  predicting  contaminant  transport  and 
ultimate  fate  in  aquatic  systems.  The  mechanisms  and  modeling  the 
transport  of  sediment  solids  are  well  described  in  the  literature  (e.g., 
Winterwerp  and  ver  Kerstern  2004;  Khan  and  Wu  2013).  The  solids 
simulation  included  in  GC  was  specifically  designed  for  use  with  describing 
the  transport  of  bound  contaminants  in  the  CSM  and  HgSM  modules. 

The  GC  module  computes  internal  sources  and  sinks  of  multiple  size  solids 
in  the  water  column  and  sediment  layer.  This  module  can  compute  the 
time  varying  input  parameters  such  as  solids  settling,  re-suspension,  and 
solids  concentrations  for  CSM  and  HgSM  modules.  Any  number  of  solids 
can  be  modeled  in  GC.  Solids  may  be  modeled  as  a  single  variable  or 
alternatively  modeled  as  multiple  independent  solids  classes.  The  number 
of  the  solids  classes  is  a  user  defined  parameter.  Each  of  the  solid  classes  is 
treated  as  an  independent  variable  with  no  interactions,  such  as  blocking 
or  flocculation.  Depending  on  the  aquatic  system  being  modeled,  these 
solids  can  be  assigned  to  mineral  abiotic  solids,  detritus,  or  various  classes 
or  size  categories  of  solids.  A  simple  kinetics  for  user-defined  constituents 
is  also  included  in  the  GC  module. 

2.1  Solids  characterization  and  conceptual  model 

Solid  particles  consist  of  an  inorganic  fraction  (e.g.,  silts,  clays)  and  an 
organic  fraction  (e.g.,  algae,  zooplankton,  bacteria,  and  detritus).  The 
inorganic  portion  is  usually  considerably  higher  than  the  organic  portion  in 
the  water  column.  Both  are  important  in  creating  turbidity  and  limiting 
light  within  the  water  column.  Natural  waters  can  contain  a  mixture  of  solid 
particles  ranging  from  gravel  (2  mm  to  20  mm)  or  sand  (0.07  mm  to  2  mm) 
down  to  very  small  particles  classified  as  silt  or  clay  (smaller  than  0.07  mm). 
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The  fine  fraction  of  solids  is  characterized  as  cohesive  sediment.  Cohesive 
sediment  can  include  silt  and  clay  particles  as  well  as  particulate  organic 
matter  such  as  detritus  and  other  forms  of  organic  carbon,  diatoms  and 
other  algae.  From  a  water  quality  perspective,  cohesive  sediments  are 
usually  of  greater  importance  in  water  quality  modeling.  In  addition  to  the 
considerations  of  flow  transport  and  solids  settling  deposition,  the  bed 
processes  of  consolidation,  bioturbation  and  re-suspension  are  also 
important  in  water  quality  modeling.  The  behavior  of  fine-grained 
suspended  solids  affects  water  quality  in  several  different  ways.  First, 
turbidity  and  its  effect  on  underwater  light  is  an  important  environmental 
condition  for  algae  growth.  The  presence  of  suspended  solids  increases  the 
attenuation  of  light  in  the  water  column,  which  leads  to  an  inhibition  of 
photosynthetic  activity;  and  therefore,  a  reduction  in  primary  production. 
Second,  the  fate  of  contaminants  in  waters  is  closely  related  to  the  amount 
of  suspended  solids,  due  to  their  large  adsorbing  capacities. 

An  overview  of  the  major  processes  of  solids  modeled  in  GC  is  shown  in 
Figure  l.  The  model  domain  is  conceptualized  as  a  well-mixed  water 
column  underlain  by  an  active  sediment  layer.  The  sediment  is  modeled  as 
a  well-mixed  surface  layer.  The  specification  of  a  mixed  surface  layer  is 
included  because  this  layer  is  often  observed  at  the  surface  of  sediments 
and  interacts  with  the  overlying  water. 


Figure  1.  Conceptual  model  representation  of  water  column  and  bed 
sediment  interactions  for  solids. 


As  shown  in  Figure  l,  solids  and  attached  contaminants  in  the  water  column 
can  be  transferred  to  the  sediment  layer  by  settling.  An  opposite  reaction, 
contaminants  can  be  released  by  sediment  re-suspension,  mixing  by  benthic 
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organisms  and  diffusion  from  the  sediment  pore  water.  A  single  sediment 
layer  with  a  constant  thickness  is  simulated  because  of  its  interactions  with 
the  overlying  water.  The  depth  of  the  benthic  sediment  layer  should  be 
chosen  based  on  the  problem.  The  three  most  important  sediment  processes 
associated  with  contaminant  modeling  are  settling,  re-suspension  or 
“erosion,”  and  “burial”  of  sediment.  In  addition  to  these  processes,  the  fate 
of  solids  in  an  aquatic  system  is  also  affected  by  physical  transport. 

The  foundation  of  contaminant  simulation  models  is  a  mass  balance,  which 
accounts  for  the  production,  loss,  and  accumulation  of  the  contaminant 
within  a  specified  control  volume.  The  rate  of  accumulation  of  mass  must 
balance  the  rates  of  production  within  the  control  volume,  input  from 
outside,  loss  across  the  boundaries,  and  loss  by  reaction.  As  stated  in  the 
introduction,  H&H  models  are  responsible  for  computing  the  physical 
transport  of  materials  in  aquatic  systems.  The  “plug  in”  water  quality 
modules  only  compute  the  internal  source  and  sink  terms  associated  with 
each  state  variable.  These  rates  contained  in  each  water  quality  module  in 
this  report  are  written  as  derivatives  of  concentration  with  respect  to  time. 

2.2  Suspended  solids 

2.2.1  Internal  sources  and  sinks 

In  the  water  column,  suspended  solids  may  settle  and  deposit  on  the 
sediment  bed.  The  settling  process  moves  solids  downward  through  the 
water  column  and  deposition  results  in  their  removal  from  the  water 
column  to  the  bed.  Sediment  re-suspension  adds  solids  into  the  water 
column.  For  a  completely  mixed  water  column  cell  with  an  underlying 
sediment  layer,  the  internal  source  (+)  and  sink  (-)  equation  of  suspended 
solids  for  each  class  is  written  as 


dmn 

dt 


Jdpn 


m„ 


where 

mn  =  concentration  of  solid  “n”  in  water  (mg  L1) 
mn2  =  concentration  of  solid  “n”  in  sediment  (mg  L1) 
Vdpn  =  deposition  velocity  of  solid  “n”  (m  d1) 

Vrpn  =  re-suspension  velocity  of  solid  “n”  (m  d1) 
h  =  water  depth  (m)  (h  >  o). 


(2.1) 
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Where  the  subscript  n  denotes  the  index  number  of  solids  class  which  is 
unitless.  Multiple  solids  classes  can  be  modeled  in  GC.  Equation  2.1 
represents  the  time  rate  of  change  of  the  solids  concentration  in  the  water 
column.  In  the  sediment  layer,  mU2  is  the  sediment  bulk  density.  The  user 
input  parameters  include  the  settling  velocities  of  solids  in  the  water 
column  and  the  sediment  solid  re-suspension  velocities.  The  settling  and 
re-suspension  parameters  depend  on  the  characteristics  of  the  solids  and 
the  flow  regime  and  can  be  varied  over  time  and  space. 

2.2.2  Settling  and  deposition 

Settling  velocities  must  be  specified  for  each  of  the  solids  classes.  These 
parameters  can  be  determined  from  field  and  laboratory  experiments,  or 
estimated  using  solids  properties.  Values  of  settling  velocities  calculated 
for  a  range  of  solids  classes  and  densities  are  provided  in  Table  1. 


Table  1.  Settling  velocities  (m  d1)  for  a  range  of  solid  sizes  and  densities. 


Solid 

Diameter  (mm) 

Solid  Density  (g  cm  3) 

1.80 

2.0 

2.5 

2.7 

Fine  sand 

0.3 

300 

400 

710 

800 

0.05 

94 

120 

180 

200 

Silt 

0.02 

15 

19 

28 

32 

0.01 

3.8 

4.7 

7.1 

8.0 

0.005 

0.94 

1.2 

1.8 

2.0 

Clay 

0.002 

0.15 

0.19 

0.28 

0.32 

0.001 

0.04 

0.05 

0.07 

0.08 

As  an  alternative,  solids  settling  velocities  can  also  be  calculated  in  the  GC 
module.  No  universal  model  of  computing  settling  velocities  is  available. 
Instead,  a  variety  of  simple  equations  have  been  used  to  estimate  the 
solids  settling  velocity  in  the  water  column.  Two  types  of  formulations  are 
implemented  in  GC:  (1)  Van  Rijn’s  formula  (Van  Rijn  1989)  and  (2) 
Cheng’s  formula  (Cheng  1997).  If  the  solids  settling  velocity  is  internally 
calculated  in  GC,  then  the  solid  parameters  listed  in  Table  2  must  be 
specified  in  the  model  simulation. 
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Option  l:  Van  Rijn’s  formula: 


v.„„  = 


0.0864  A gd 


0.065  mm  <  dpn  <  0.1  mm 


18  v 


8.64x10s— fJl.O  +  O.Old3  -1.0)’  0.1  mm  < dpn  <  1  mm 
d 


pn 


&gdp 


1.1  J  -  pn  x  86400 
1000 


dpn  >  1mm 


and 


A  = 


(2.2b) 


d. 


_  dP» 

lAg] 

1000 

\v2) 

1/3 


(2.2c) 


where 

Vspn  =  settling  velocity  of  solid  “n”  (m  d_1) 
g  =  gravity  acceleration  (9.81  m  s-2) 
v  =  kinematic  viscosity  of  water  (m2  s_1) 
pPn  =  density  of  solid  “n”  in  water  (g  cm  3) 
pw  =  density  of  water  (1.0  g  cm  3) 
dpn  =  diameter  of  solid  “n”  (mm) 
d*  =  dimensionless  solids  parameter(unitless). 


The  second  equation  is  Cheng  (1979).  On  the  basis  of  experimental  data 
for  natural  particles,  Cheng  (1997)  developed  the  following  equation  that 
approximates  solids  settling  velocities 


v  =8.64) 

spn 


107 


y]25  +  1.2d,2  -5.0 


(2.3) 


Kinematic  viscosity  versus  water  temperature  can  be  estimated  from  a 
regression  relationship 


_ 1.79 -10~6 _ 

1.0  +  0.033687;,  ^  0.00022  lTj 


(2.4) 


ERDC/EL  TR-16-8 


9 


where 


Tw  =  water  temperature  (°C). 

In  the  GC  module,  the  settling  velocities  are  set  to  zero,  when  the  shear 
stress  exceeds  a  certain  critical  value,  or  when  the  water  depth  is  smaller 
than  a  certain  critical  depth  (Krone  1962). 


For  cohesive  sediments,  the  settling  behavior  may  not  be  non-interactive. 
For  example,  it  has  been  found  in  both  the  laboratory  and  field  that 
settlement  of  cohesive  solids  is  often  influenced  by  suspended  sediment 
concentrations  and  flocculation  (both  coagulation  and  the  break-up  of 
particles).  Additionally,  suspended  solids  in  natural  systems  exhibit 
variations  in  drag  characteristics  and  settling  behavior  due  to  their  non- 
uniform  shapes.  The  above  equations  may  not  adequately  capture  these 
variations.  Therefore,  model  calibration  may  be  needed  when  these 
formulae  are  applied  to  different  sites  or  conditions. 


For  solids  settling  onto  the  sediment  bed,  a  distinction  may  be  made 
between  settling  and  deposition.  A  portion  of  the  solids  which  settle 
through  the  water  column  may  remain  in  suspension  instead  of  depositing 
onto  the  surface  of  the  sediment  bed.  Whether  a  solid  will  settle  to  the  bed 
sediment  layer  depends  on  its  physical  and  chemical  properties  as  well  as 
the  hydraulic  conditions  over  the  depth  of  the  water  column.  Deposition 
onto,  and  attachment  to  the  sediment  bed,  are  usually  described  as 
probabilistic  processes  (Krone  1962).  Ariathuri  and  Krone  (1976)  used  this 
concept  to  develop  a  linear  probability  of  deposition  function  to  determine 
what  fraction  of  a  particulate  class  of  solids  is  truly  depositional 


p  — 

rdpn  ~ 


1.0 


T  <  X 


cdl 


10 _ - _ Tcdl  1  T  <T<T 

Lcdl  —  1  —  lcdu 


0 


X  >  T 


cdu 


(2.5) 


where 

r  =  bottom  shear  stress  (N  nr2) 

Tcdi  =  lower  critical  shear  stress  for  deposition  (N  nr2) 
Tcdu  =  upper  critical  shear  stress  for  deposition  (N  nr2). 
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Lower  critical  shear  stress  for  deposition  (r Cdi)  is  below  which  all  solids 
have  a  full  probability  to  deposit  on  the  bed.  Upper  critical  shear  stress  for 
deposition  (r cdu)  is  above  which  all  solids  remain  in  suspension  yielding  a 
zero  deposition  rate.  Once  the  settling  velocity  is  determined,  the 
deposition  velocity  can  be  calculated  as  the  product  of  the  settling  velocity 
and  the  probability  of  deposition 

Vdpn  =  PdpnVspn  (Z6) 


where 

Pdpn  =  deposition  probability  of  solid  “n”  (o  -  l.o). 

2.3  Sediment  layer 

When  suspended  solids  settle  to  the  bed  of  a  water  body,  they  become 
sediments.  Sediment  is  composed  of  many  materials,  including  individual 
primary  particles,  aggregates,  organic  materials,  and  associated  chemicals. 
Particles  can  be  mineral  or  organic  in  origin.  The  role  of  sediment  in  water 
quality  modeling  is  tied  both  to  the  particle  size  of  sediment,  and  to  the 
amount  of  particulate  organic  carbon  associated  with  the  sediment. 

The  mass  balance  of  solids  must  be  tracked  for  the  sediment  layer  in  the 
underlying  the  water  column.  The  GC  uses  a  mass  balance  equation  to 
calculate  solids  mass  and  concentrations  for  underlying  sediment  layer.  The 
sediment  layer  is  assumed  in  GC  as  homogeneous  (well  mixed)  and  the 
model  does  not  change  the  thickness  of  sediment  layer  during  the 
simulation.  As  the  sediment  layer  accepts  a  new  sediment  increment  on  the 
top,  an  equal  amount  is  removed  from  this  layer  for  burial.  When  sediment 
re-suspension  occurs,  it  is  assumed  that  deeper  sediments  below  this  layer 
are  entrained  into  the  upper  layer  and  then  mixed  with  the  remaining 
sediments.  The  density  of  the  sediment  layer  can  vary,  depending  on  the 
variable  sediment  composition.  The  porosity  within  the  sediment  layer  is 
assumed  constant  and  user  defined.  The  solids  simulation  module  also 
requires  a  set  of  initial  conditions  of  the  sediment  solids. 

2.3.1  Solids  mass  balance 

The  mass  balance  of  solids  in  the  sediment  layer  is  a  flux  balance  between 
the  incoming  solids  due  to  deposition  from  the  water  column,  and  loss  due 
to  sediment  re-suspension  and  burial.  Deposition  from  the  water  column 
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adds  the  source  of  solids  to  the  sediment  layer.  Re-suspension  and  deep 
burial  result  in  a  loss  of  solids  in  the  sediment  layer.  Exchange  of  solids 
between  the  water  column  and  the  sediments  are  simulated  via  the 
deposition  and  re-suspension  processes.  The  mass  balance  of  sediment 
solids  for  each  class  is  written  as 


dmn: 

dt 


J dpn 


777  —  - 


K 


-m,. 


(2.7) 


where 

vb  =  sediment  burial  velocity  (m  d1) 

h.2  =  sediment  layer  thickness  (m)  ( h2  >  o). 

Once  concentrations  of  the  sediment  layer  solids  for  each  class  are  solved 
from  equation  2.7.  The  interaction  between  solids  in  the  water  column  and 
active  sediment  layer  can  be  appropriately  quantified  in  GC. 

2.3.2  Sediment  re-suspension 

Suspended  solids  settle  through  the  water  column  to  the  surficial  sediment 
layer  and  can  be  re-introduced  back  to  overlying  water  through  re¬ 
suspension.  The  re-suspension  and  erosion  of  sediments  refers  to  the 
process  by  which  solids,  and  their  attached  contaminants,  are  transported 
from  bed  sediments  into  the  water  column.  As  with  settling  and  deposition, 
sediment  re-suspension  may  be  influenced  by  a  variety  of  physical  and 
chemical  properties  within  the  bed  and  the  overlying  water  column.  Re¬ 
suspension  of  bottom  sediments  can  be  induced  by  an  increase  of  hydro- 
dynamic  stress  or  a  weakening  of  sediment  resistance,  depending  on  both 
hydrodynamics  and  the  physicochemical  characteristics  of  the  sediment. 
There  are  two  distinct  erodibility  attributes  for  a  given  solid.  These 
attributes  are  critical  shear  stress  (i.e.,  the  critical  level  of  flow  shear  stress 
that  is  sufficient  to  dislodge  the  sediment  solids),  and  the  erosion  rate  (i.e., 
the  erosion  velocity  of  bed  sediment  into  the  water  column).  Both  attributes 
are  a  function  of  shear  stress.  The  rate  of  erosion  varies  with  the  energy 
available  for  erosion,  shear  stress,  and  the  characteristics  of  the  sediments. 

The  sediment  re-suspension  velocity  is  a  user-defined  input  parameter. 
Alternatively,  given  bottom  shear  stresses,  the  GC  can  compute  the  re¬ 
suspension  velocity  using  one  of  three  equations:  1)  Lick  et  al.  (1995),  2) 
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Parchure  and  Mehta  (1985),  and  3)  Lick  (2009).  The  first  two  equations  are 
applied  for  cohesive  sediments. 

The  first  equation  is  based  on  the  equation  given  by  Lick  et  al.  (1995)  in 
which  the  critical  shear  stress  is  the  major  parameter  controlling  the  re¬ 
suspension  mechanism. 


E  =  P 


— — 1.0 
T 


(2.8) 


where 

E  =  sediment  erosion  rate  (g  cm  2  s_1) 

P  =  surface  erosion  rate  (g  cnr2  s-1) 

Tee  =  critical  shear  stress  for  erosion  (N  nr2) 
m  =  empirically  determined  constant  (unitless). 

The  second  equation  is  based  on  the  equation  given  by  Parchure  and 
Mehta  (1985) 


E  =  E0e 


“o(t 


(2.9) 


where 

Eo  =  surface  erosion  rate  (g  cnr2  s-1) 

cio  =  empirical  constant  (unitless). 

The  empirical  constants  in  the  above  two  equations  used  for  computing 
the  re-suspension  rate  are  site  specific.  The  erodibility  coefficients  P  or  Eo, 
can  be  a  function  of  the  degree  of  consolidation  of  sediments,  their  size 
distribution  and  factors  affecting  their  stability  (such  as  armoring  due  to 
the  existence  of  mixtures  of  large  and  fine  particles).  Because  these 
properties  are  poorly  known,  this  entry  can  simply  be  considered  as  a 
calibration  parameter. 

The  critical  shear  stress  defines  the  point  at  which  the  erosion  of  bed 
sediments  is  initiated.  Erosion  or  re-suspension  of  bed  sediments  is 
predicted  to  occur  when  the  computed  shear  stress  is  greater  than  the 
critical  sheer  stress  (r  >  Tee).  No  erosion  will  occur  for  a  shear  stress  less 
than  the  critical  shear  stress  (t  <  Tee).  The  critical  shear  stress  for  erosion  is 
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a  function  of  the  degree  of  compaction  of  the  sediments.  The  critical  shear 
stress  must  either  be  measured  or  estimated  by  model  calibration.  Typical 
values  of  critical  shear  stresses  are  on  the  order  of  o.i  -  0.4  N  nr2. 

Once  the  sediment  re-suspension  rate  is  known,  re-suspension  velocity  for 
each  solids  class  is  calculated  by 

V  r  =8.64-l08—  (2.10) 

m„. 


The  third  equation  used  in  the  GC  module  is  based  on  a  general  equation 
given  by  Lick  (2009)  and  can  be  applied  for  both  cohesive  and  non-cohesive 
sediment  solids.  This  equation  directly  computes  the  re-suspension  velocity 


v.  =  10 


T  —  T 


(2.11) 


where 

Vr  =  sediment  re-suspension  velocity  of  solid  “n”  (cm  d1) 

Tm  =  critical  shear  stress  for  the  initiation  of  movement  of 
noncohensive  solids  (N  nr2). 

The  critical  shear  stress  (rcn)  for  the  initiation  of  movement  of 
noncohensive  solids  is  specified  as  an  input  parameter.  This  parameter  has 
been  determined  from  experiments  and  more  refined  theoretical  analyses. 

2.3.3  Bottom  shear  stress 

The  resuspension  formulations  described  above  all  relate  the  resuspension 
rate  constant  for  cohesive  sediments  to  the  difference  between  the  bottom 
shear  stress  (t)  and  a  critical  shear  stress  (rce).  The  bottom  stress  might  be 
induced  by  flow  velocities,  tidal  and  wind-induced  advective  flows,  and 
surface  waves.  The  bottom  shear  stress  may  be  estimated  from  wind 
speeds  and  direction  and/or  current  velocities  using  a  bottom  boundary 
layer  model,  such  as  Kajiura’s  (1968).  Otherwise,  this  parameter  can  be 
obtained  from  the  H&H  models. 
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2.3.4  Sediment  burial 

Sediment  burial  velocity  is  a  user  defined  term  used  to  convey  that  bed 
sediment  is  no  longer  available  for  re-suspension  because  it  has  been 
covered  by  much  newer  sediment  (i.e.,  it  has  been  buried.).  In  GC, 
sediment  burial  velocity  can  also  be  calculated  according  to  the  following 
mass  balance  for  the  sediment  layer  solids 

Vb  (!  -  9)  P.s 106  =  ^  Vdpnmn  -Y,(Vrpnmn2)  (2  -12) 


where 

cp  =  sediment  layer  porosity  (unitless) 
ps  =  sediment  dry  density  (g  cm  3). 

Equation  2.12  can  be  used  to  calculate  the  burial  velocity  ( Vb )  in  GC.  Within 
this  equation  suspended  and  sediment  solids  concentrations  are  determined 
respectively  from  equations  2.1  and  2.7. 

2.4  Solids  and  sediment  layer  parameters 

The  “plug  in”  water  quality  module  operates  based  on  the  schematization 
already  set  up  for  the  H&H  model  and  the  flows  computed  by  the  H&H 
model.  This  means  that  the  user  only  has  to  specify  a  limited  amount  of 
water  quality  data: 

•  definition  of  the  water  quality  cells 

•  initial  concentrations  of  the  state  variables 

•  concentrations  of  all  state  variables  on  the  inflow  boundaries 

•  forcing  functions 

•  water  quality  parameters. 

Table  2  summarizes  the  solids  and  sediment  layer  input  parameters  that 
must  be  specified  in  the  GC  module  when  all  options  are  turned  on. 

Typical  values  are  taken  from  a  limited  literature  review.  All  of  these 
parameters  are  spatially-variable  parameters  and  do  not  vary  with  time. 
They  are  treated  as  constant  values  throughout  the  model  simulation. 
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Table  2.  List  of  solids  and  user-defined  constituent  parameters. 


Symbol3 

Definition 

Typical  values 

Units 

koi(T) 

Zero-order  decay  rate  of  constituent  i 

n/ab 

mg  L_1d“1 

kn(T) 

First-order  decay  rate  of  constituent  i 

n/a 

d-1 

Vs  i 

Settling  velocity  of  constituent  i 

n/a 

m  d-1 

n 

Solids  classes 

1-5 

- 

Vs  pn 

Solids  settling  velocity 

n/a 

m  d'1 

dpn 

Solids  diameter 

n/a 

mm 

Tcdu 

Upper  critical  shear  stress  for  deposition 

0.06  -  1.1 

N  nr2 

Tcdl 

Lower  critical  shear  stress  for  deposition 

0.2d 

N  nr2 

h2 

Sediment  layer  thickness 

0.1 

m 

<t> 

Sediment  layer  porosity 

0.3  -  0.9 

- 

Ps 

Sediment  dry  density 

1.6 

g  cm3 

Vrpn 

Re-suspension  velocity 

n/a 

m  d-1 

Tee 

Critical  shear  stress  for  erosion 

0.4C 

N  nr2 

Ten 

Critical  shear  stress  for  the  initiation  of 
movement  of  noncohensive  solids 

0.414-103  d 

N  nr3 

CXO 

Empirical  constant 

0.5 

- 

p 

Surface  erosion  rate 

1.3  -  4.7-10'5 

5-10'6  c 

g  cm-2  s1 

m 

Empirically  determined  constant 

1 

- 

Eo 

Surface  erosion  rate 

n/a 

g  cm-2  s*1 

Vb 

Sediment  burial  velocity 

n/a 

m  yr1 

a.  Subscript  i  denotes  the  index  number  of  user-defined  constituents,  subscript  n  denotes  the  index 
number  of  solid  classes. 

b.  n/a  represents  not  available. 

c.  Mulder  and  Udink  (1991). 

d.  Winterwerp  et  al.  (1993). 


Suspended  solids  parameters  consist  of  a  number  of  solids  classes  and  the 
settling  velocities  for  each  class.  Inclusion  of  the  sediment  layer  introduces 
several  additional  sediment  parameters.  In  addition  to  a  layer  thickness, 
the  other  sediment  parameters  include  sediment  re-suspension  velocity, 
sediment  burial  rate,  sediment  bulk  density  (pt>),  sediment  layer  porosity 
(0),  and  the  percent  by  weight  of  each  solid  class.  Sediment  bulk  density 
( pb )  is  also  the  layer  total  solids  concentration  (i.e.,  layer  total  solids  mass 
divided  by  layer  volume);  the  total,  or  mixed  solids  density  (i.e.,  density  of 
solids),  is  the  layer  total  solids  mass  divided  by  layer  total  solids  volume; 
and  sediment  layer  porosity  ($)  is  the  ratio  of  voids  volume  (i.e.,  water 
volume  for  saturated  conditions  in  sediments)  to  total  volume  (i.e.,  water 
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plus  solids  volume).  Solids  density  averages  about  2.65  g  m  3  for  most 
sediments  (sands,  silts,  and  clays).  However,  solids  density  for  carbon- 
occluded  minerals  ranges  from  1.5  to  2.2  g  nr3  (Di  Toro  2001).  The 
porosity  of  the  upper  few  centimeters  of  the  sediment  layer  typically 
ranges  from  0.7  to  0.9  or  higher.  Sediment  porosity  is  known  to  decrease 
with  solid  particle  diameter,  where  values  range  from  approximately  0.7  to 
0.9  for  0.001  mm  to  0.3  for  1.0  mm  (Di  Toro  2001). 

2.5  Simple  kinetics  for  user-defined  constituents 

In  addition  to  solids,  the  GC  can  model  any  number  of  user-defined 
constituents.  These  constituents  are  modeled  with  simple  kinetics  and 
only  for  the  water  column.  Three  processes  are  included  in  the  GC  module 
for  any  user-defined  constituent:  zero-order  decay,  first-order  decay,  and 
net  settling  loss.  The  user  defines  a  temperature  dependent  zero  -  or  first- 
order  decay  rate,  or  a  net  settling  loss  rate  for  modeled  constituents.  The 
net  settling  rate  of  constituents  can  be  negative,  depending  on  the  degree 
of  sediment  re-suspension. 

For  any  user  defined  constituent  in  the  water  column,  the  internal  source 
(+)  and  sink  (-)  equation  is 

^  =  -koi(T)- ku ( T)Ct  -^C;  (2.15) 


where 

G  =  concentration  of  constituent  i  (mg  Lr1) 
koi(T)  =  zero-order  decay  rate  of  constituent  i  (mg  L"1  d”1) 
kii(T)  =  first-order  decay  rate  of  constituent  i  (d_1) 

Vsi  =  net  settling  loss  rate  of  constituent  i  (m  d-1). 


Where  the  subscript  i  denotes  the  index  number  of  user-defined 
constituents.  The  GC  allows  the  user  to  specify  simple  reaction  rates  for 
each  of  the  contaminants  modeled.  More  complex  reactions  and 
transformations  are  simulated  in  CSM.  User  defined  constituent 
parameters  are  listed  in  Table  2.  Table  2.  These  are  spatially-variable 
parameters,  do  not  vary  with  time,  and  they  are  treated  as  constant  values 
during  the  model  simulation. 
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2.6  Water  temperature  correction  equations 

The  reaction  coefficients  used  in  GC,  CSM,  and  HgSM  modules  are  usually 
specified  by  the  user  at  20°  C.  These  coefficients  must  be  corrected  based 
on  the  temperature  difference  between  the  actual  temperature  and  a 
reference  temperature  (usually  20°  C).  Three  temperature  correction 
formulations  are  used  in  the  GC,  CSM,  and  HgSM  modules  to  correct  the 
temperature  dependent  coefficients.  The  first  one  is  the  Arrhenius 
Equation  (AE).  The  second  one  is  Qio  equation.  The  third  one  is  the 
modified  Arrhenius  Equation  (MAE). 

2.6.1  Arrhenius  Equation  (AE) 

The  general  form  of  Arrhenius  Equation  (Fogler  2005)  can  be  written  as 

k(T)  =  k(Twr)ea]^  (2.16) 


where 

k(T)  =  value  of  the  rate  at  local  temperature  (d1) 
k(Twr)  =  value  of  the  rate  at  reference  temperature  (d-1) 

Ea  =  activation  energy  (J  mol1) 

R  =  universal  gas  constant  (8.314  J  mol1  K1) 

Twr  =  reference  water  temperature  (°K). 

For  some  chemicals,  the  Qio  equation  is  also  used  to  correct  temperature 
dependent  reaction  coefficients  (Behradek  1930) 

k(T)  =  k(Twr)Q^  (2.17) 

In  this  equation,  Ow  gives  the  change  of  reaction  coefficient  for  every  io°  C 
temperature  change. 


The  Q  10  equation  is  based  on  an  approximation  of  the  AE,  Qio  is  defined  by 
(Jonsson  and  Agerberg  2015) 


(2.18) 
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2.6.2  Modified  Arrhenius  Equation  (MAE) 

The  MAE  makes  explicit  the  temperature  dependence  of  the  exponential 
factor.  The  MAE  is  usually  of  the  form 

k(T)  =  k(TwrW-‘  T-  (2.19) 


where 


6  =  temperature  correction  coefficient. 

In  contrast  to  the  AE,  the  temperature  coefficient  6*  is  needed  in  the  MAE, 
which  usually  ranges  between  l.oi  and  l.io.  The  relationship  of  the  rate 
coefficient  and  water  temperature  defined  in  this  function  is  illustrated  in 
Figure  2. 


Figure  2.  Relationship  of  the  rate  coefficient  and  temperature  as  defined  in  modified 

Arrhenius  Equation. 
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3  Contaminant  Simulation  Module  (CSM) 

This  chapter  describes  the  theory  and  mathematical  formulations 
implemented  in  CSM.  The  CSM  is  capable  of  modeling  contaminants  in  an 
aquatic  system  as  influenced  by  the  following  processes;  ionization,  multi¬ 
phase  partitioning,  degradation,  photolysis,  hydrolysis,  volatilization, 
generalized  second-order  reaction,  and  transformations  where  one  chemical 
undergoes  a  reaction  and  is  transformed  to  a  daughter  product.  Any  process 
in  CSM  can  be  ignored  by  use  of  switches  where  such  processes  are  not 
applicable.  Two  types  of  contaminant  partitioning  options  are  included  for 
algae  and  solid  particulates;  equilibrium  and  non-equilibrium  in  which 
adsorption/desorption  can  be  affected  by  rate  limiting  processes.  The  CSM 
models  the  water  column  exchange  with  underlying  sediments  and 
exchange  with  the  atmosphere  through  a  volatilization  process.  The  CSM 
can  model  multiple  contaminants  in  one  simulation.  The  contaminants 
themselves  are  arbitrary,  in  that  the  specific  contaminant  to  be  simulated  is 
defined  through  the  specification  of  processes  and  kinetic  rates.  All 
concentrations  of  contaminants  in  CSM  are  expressed  in  terms  of 
micrograms  per  volumetric  units  of  liter  (pg  L1)- 

3.1  Overview 

The  transport  and  fate  of  contaminants  in  aquatic  systems  involve  four 
media;  atmosphere,  water,  sediment,  and  biota.  The  importance  of  each  of 
these  media  depends  on  the  chemical  properties  of  the  contaminant  and  the 
form  in  which  the  contaminant  enters  the  environment.  The  atmosphere 
can  be  an  important  source  of  contaminants  to  the  water  through  both  dry 
and  wet  deposition,  and  can  act  as  a  sink  for  more  volatile  contaminants. 
The  bottom  sediment  is  often  a  repository  for  contaminants  discharged  or 
released  into  the  aquatic  environment.  Sediment  re-suspension,  burial  and 
bioturbation,  reactions,  and  partitioning  between  the  water  and  sediment 
all  affect  the  fate  of  contaminants  in  aquatic  systems.  The  biological  and 
chemical  processes  can  alter  the  form  and  the  rate  of  the  contaminant. 

As  discussed  in  Chapter  l,  the  CSM  computes  kinetics  of  contaminants  in 
aquatic  systems.  At  a  minimum,  the  CSM  needs  to  simulate  the  water 
column  and  a  bed  layer,  and  include  both  chemical  reactions  and  sorption 
to  solids.  The  conceptual  representation  of  a  water  column  and  sediment 
contaminant  interactions  and  processes  modeled  in  CSM  is  depicted  in 
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Figure  3.  It  is  important  to  note  that  Figure  3  does  not  show  all  processes 
and  fluxes  of  the  contaminant  cycle,  but  only  the  most  relevant  for  the  water 
column  and  benthic  sediment  layer.  Each  contaminant  in  the  water  column 
is  subject  to  adsorption  and  desorption  with  dissolved  organic  carbon 
(DOC)  and  solids.  The  dissolved  phase  in  the  bulk  water  (aqueous  phase), 
the  adsorbed  phase  to  DOC  in  the  bulk  water,  and  the  adsorbed  phases  to 
organic  and  inorganic  solids  are  modeled  in  CSM.  The  chemical  species  in 
the  active  sediment  layer  is  also  partitioned  into  corresponding  forms:  the 
truly  dissolved  phase  in  the  pore  water,  the  adsorbed  phase  to  DOC  in  the 
pore  water,  and  the  adsorbed  phase  to  the  sediments.  Two  types  of 
partitioning  (equilibrium  and  non-equilibrium)  are  included  for  algae; 
particulate  organic  matter  (POM),  and  solids.  Contaminants  can  enter  the 
water  column  via  surface  runoff,  atmospheric  deposition,  and/or  direct 
discharge.  Once  contaminants  are  in  aquatic  systems,  they  may  be  degraded 
or  transformed  by  various  processes. 


Figure  3.  Conceptual  representation  of  contaminant  processes  modeled  in  CSM. 


The  contaminant  in  the  water  column  is  exchanged  with  benthic 
sediments  by  dispersive  mixing.  Sorbed  contaminant  settles  through  water 
column  and  deposits  to  or  erodes  from  benthic  sediments.  Within  the  bed, 
dissolved  contaminant  migrates  downward  or  upward  through  percolation 
and  pore  water  diffusion.  Sorbed  chemical  migrates  downward  or  upward 
through  net  sedimentation  or  erosion.  All  of  these  processes  are  simulated 
in  CSM  and  are  combined  with  contaminant  reactions  and 
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transformations  into  internal  source  and  sink  terms.  In  CSM,  the  sediment 
layer  represents  the  upper  mixed  layer  and  may  be  on  the  order  of  several 
centimeters  thick  due  to  bioturbation  and  mixing.  The  conceptualization  is 
based  on  the  observation  that  bioturbation  forms  a  well-mixed  sediment 
layer  in  a  variety  of  locations  and  environments  (Boudreau  1998).  The 
single,  well-mixed  layer  forms  the  basis  for  some  of  the  earliest  toxic 
models  (O’Connor  et  al.  1983).  For  example,  the  RECOVERY  model  (Boyer 
et  al.  1994;  Ruiz  et  al.  2000)  places  a  well-mixed  layer  at  the  sediment- 
water  interface,  above  a  succession  of  deeper  sediment  layers. 

A  variety  of  toxic  contaminant  transport  and  fate  models  are  available.  Over 
the  last  few  decades,  considerable  effort  has  been  made  to  describe  and 
model  contaminant  transport,  reactions,  transformations,  and  bio¬ 
availability  in  aquatic  ecosystems  (Wool  et  al.  2006).  Model  representations 
of  organic  contaminants  in  the  water  column  are  usually  similar  with  regard 
to  kinetic  processes  and  formulations;  however,  representations  of  the 
bottom  sediment  vary  widely.  The  biogeochemical  reactions  and 
transformations  implemented  in  CSM  build  upon  previous  contaminant 
modeling  efforts.  The  following  processes  are  implemented  in  CSM: 
ionization,  multi-phase  partitioning,  degradation,  photolysis,  hydrolysis, 
volatilization,  generalized  second-order  reaction,  and  transformations.  Each 
process  is  controlled  by  a  on/off  switch.  The  relative  importance  of  each  of 
the  processes  included  in  the  model  is  directly  governed  by  the  contaminant 
species  and  their  associated  properties.  Each  of  these  processes  is  discussed 
in  detail  below. 

3.2  Contaminant  ionization 

Ionization  is  the  dissociation  of  a  contaminant  into  multiple  charged 
chemical  species.  In  an  aquatic  environment  some  contaminants  may 
occur  only  in  their  neutral  forms  while  others  may  react  with  water 
molecules  to  form  positively  (cationic)  or  negatively  (anionic)  charged 
ions.  Ionization  can  be  important  because  of  the  different  toxicological  and 
chemical  properties  of  the  neutral  and  ionized  species.  For  example,  in 
some  cases  only  the  neutral  form  of  the  chemical  may  react  or  be 
transported  through  biotic  membranes  resulting  in  toxicity.  As  a  result,  it 
is  often  necessary  to  compute  the  distribution  of  the  contaminant  among 
the  ionic  forms  to  allow  them  to  react  or  transform  at  different  rates. 

In  CSM,  the  ionization  process  simulates  the  contaminant  for  as  many  as 
five  species:  1)  the  neutral  molecule,  2)  singly  charged  cations,  3)  doubly 
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charged  cations,  4)  singly  charged  anions,  and  5)  doubly  charged  anions. 
Each  of  the  neutral  or  ionic  species  may  also  occur  in  the  dissolved  phase  or 
adsorbed  to  DOC,  algae,  POM,  or  the  multiple  solids.  Each  chemical  species 
may  have  different  reactivity  as  reflected  by  different  transformation  rates. 
The  formation  of  the  contaminant  is  controlled  by  the  user  specified  inputs. 
Ionization  reactions  are  rapid  and  are  generally  assumed  to  be  at  (local) 
equilibrium.  At  equilibrium,  the  distribution  of  the  contaminant  between 
the  neutral  and  the  ionized  species  is  controlled  by  the  pH  and  temperature 
of  the  water  and  the  ionization  constants.  Given  the  total  concentration,  the 
pH,  and  the  equilibrium  constants,  the  fractions  of  the  contaminant 
occurring  in  each  of  the  chemical  species  k,  respectively,  are  defined  as 
(Wool  et  al.  2006) 


/■= i 

D 

(3.1a) 

2  _  Kj  10^-14 

J  D 

(3.1b) 

3  KblKb2/(l0^f 

J  D 

(3.1c) 

4  KjlO  pH 

J  D 

(3. Id) 

5  _  KalKa2/(l0  ^f 

J  D 
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-pH 


10 


-pH 


+  1 


(3.2) 


where 


f  =  fraction  of  the  chemical  occurring  in  each  of  the  chemical 
species  i  (unitless) 

Ki  =  equilibrium  constant  for  the  formation  of  the  acid  (Kai),  or 

anionic  species,  or  the  base  (Kbi)  or  cationic  species  (unitless). 
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Where  superscript  i  is  the  chemical  species  index  (z  =  l  to  5).  When  the 
ionization  is  simulated  in  CSM,  the  five  chemical  species  are  defined  for 
each  contaminant.  Different  partition  coefficients  and  reaction  rates  (e.g., 
for  degradation,  hydrolysis,  photolysis,  volatilization,  etc.)  may  be 
specified  for  each  ionic  species  of  the  contaminant  modeled  in  CSM. 

3.3  Contaminant  partitioning  and  distribution 

The  partitioning  process  is  used  to  describe  the  distribution  of  a 
contaminant  between  pairs  of  media.  Partitioning,  also  called  sorption,  is 
the  bonding  of  dissolved  chemicals  onto  solid  phases,  such  as  suspended 
solids,  biological  material,  and  sometimes  dissolved  or  colloidal  organic 
material.  The  distribution  between  dissolved  and  adsorbed  forms  impacts 
the  transport  and  fate  of  the  contaminants  in  aquatic  systems.  Partitioning 
of  the  contaminant  is  often  associated  with  DOC  and  the  organic  carbon 
content  of  the  solids  (Karickhoff  1984;  Di  Toro  et  al.  1991).  The  sorption  of 
xenobiotic  organic  compounds  to  DOC  is  a  driving  force  in  determining  the 
bioavailability  (Traina  et  al.  1996).  Therefore,  multiple  phase  partitioning 
(dissolved  in  water,  adsorbed  to  DOC,  and  adsorbed  to  solids)  is  modeled  in 
CSM.  Organic  solids  are  derived  from  plant  materials,  dead  bacterial  or 
algal  cells,  and  decaying  aquatic  organisms.  Depending  on  the  aquatic 
system  being  modeled,  solids  can  be  assigned  to  mineral  abiotic  solids, 
detrital,  or  various  classes  or  size  categories  of  solids.  Sorption  or  uptake  of 
the  contaminant  associated  with  algae  in  the  water  column  is  also  included 
in  CSM.  Inside  the  algae  contaminant  kinetics  and  transformations  are  not 
modeled  in  CSM. 

Slow  diffusion  in  solid  matter  has  been  acknowledged  to  take  place  after 
fast  equilibrium  adsorption  or  prior  to  fast  equilibrium  desorption  (Runkel 
et  al.  1999;  Zheng  and  Bennett  2002).  An  alternative  to  equilibrium 
partitioning  is  to  simulate  sorption  using  kinetic  rate  equations.  In  CSM, 
partitioning  can  be  simulated  with  two  options:  1)  equilibrium  approaches, 
and  2)  non-equilibrium  adsorption/desorption  kinetics.  Eiither  nonlinear 
sorption  isotherms  or  non-equilibrium  sorption  and  desorption  can  be 
simulated  in  CSM  for  contaminants.  Some  limitations  should  be  kept  in 
mind  when  applying  the  equilibrium  partitioning.  First,  contaminant 
concentrations  should  be  near  trace  levels,  (i.e.,  below  half  the  solubility). 
At  higher  concentrations,  the  assumptions  of  linear  partitioning  and 
transformation  begin  to  break  down. 
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3.3.1  Equilibrium  partitioning 

3.3. 1.1  Partitioning  isotherms 

When  sorption  reactions  are  fast  relative  to  other  environmental 
processes,  then  an  equilibrium  partitioning  isotherm  is  often  used  in  the 
modeling  contaminant  transport  and  fate  in  aquatic  systems.  Under 
equilibrium  partitioning,  adsorption  and  desorption  can  be  simulated  as 
either  a  linear  or  non-linear  function  of  organic  matter  and  solids. 
Therefore,  three  commonly  reported  sorption  isotherms  are  implemented 
in  CSM:  l)  Linear,  2)  Langmuir,  and  3)  Freundlich.  Figure  4  shows  the 
adsorbed  concentration  as  a  function  of  aqueous  concentration  defined  by 
three  isotherms.  The  input  coefficients  and  requirements  for  the  linear, 
Freundlich,  and  Langmuir  isotherm  equations  are  different. 

Figure  4.  Three  equilibrium  sorption  isotherms. 


The  sorption  isotherm  indicates  how  the  contaminant  distributes  between 
the  aqueous  phase  and  the  solid  phase  when  the  sorption  process  reaches  an 
equilibrium  state.  Adsorption  isotherms  at  very  low  solute  concentrations 
are  often  linear. 

•  Linear  isotherm:  The  linear  partition  coefficient  is  simply  a  ratio  of  the 
adsorbed  phase  concentration  (expressed  in  mg  per  kg  adsorbing 
material)  to  the  dissolved  phase  concentration  (expressed  in  pg  per  L  of 
solution)  at  equilibrium.  The  linear  isotherm  is  defined  as 
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Cps  =  KdCd 


10 


(3.3) 


where 

CPs  =  mass  of  solids  adsorbed  phase  in  water  (pg  g-1) 

Kd  =  sorption  or  partition  coefficient  (L  kg1) 

Cd  =  concentration  of  dissolved  phase  in  water  (pg  L1). 

The  Kd  value  is  a  direct  measure  of  the  partitioning  of  a  contaminant 
between  the  solid  and  aqueous  phases.  The  unit  of  Kd  is  therefore  L  Kg1. 
Values  for  Kd  not  only  vary  greatly  between  contaminants,  but  also  vary  as 
a  function  of  aqueous  and  solid  phase  chemistry  (Delegard  and  Barney 
1983;  Kaplan  and  Serne  1995).  Kd  for  partitioning  between  suspended 
solids  and  surface  water,  and  Kd  for  partitioning  between  sediment  solids 
and  its  pore  water  can  be  defined  as  different  values  in  CSM.  Linear 
adsorption  best  describes  the  partitioning  at  low  concentrations. 

•  Freundlich  isotherm:  Contaminant  adsorption  on  solids  can  deviate 
from  the  above  linear  relationship.  The  Freundlich  isotherm 
(Freundlich  1926)  has  been  widely  used  to  describe  solute  adsorption 
by  solids.  The  Freundlich  isotherm  is  defined  as 

Cp,  =  KJ(C,f  (3.4) 


where 


Kf  =  Freundlich  adsorption  constant  ([pg  g1]  [pg  L1]^) 
b  =  Freundlich  exponent  (unitless). 

The  Freundlich  equation  is  sometimes  written  with  the  exponent  in 
equation  3.4  being  l/b  instead  of  b.  The  linear  isotherm  is  a  special  case  of 
the  Freundlich  isotherm  where  the  Freundlich  exponent  b  is  equal  to  1. 
When  b  >  1,  the  adsorption  rate  increases  with  increasing  solution 
concentration.  When  b  <  1,  the  adsorption  rate  decreases  with  solution 
concentration,  as  the  low  energy  sites  are  occupied.  The  Freundlich  model 
does  not  account  for  finite  adsorption  capacity  at  high  concentrations  of 
solute.  However,  when  considering  trace  constituent  adsorption,  ignoring 
such  physical  constraints  is  usually  not  critical. 
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•  Langmuir  isotherm:  The  Langmuir  model  was  originally  proposed  to 
describe  the  adsorption  of  gas  molecules  onto  homogeneous  solid 
surfaces  (crystalline  materials)  that  exhibit  one  type  of  adsorption  site 
(Langmuir  1918).  The  Langmuir  adsorption  model  has  been  extended 
to  describe  adsorption  of  solution  species  onto  solid  adsorbents 
including  heterogeneous  solids.  The  Langmuir  isotherm  assumes  that  a 
saturation  point,  adsorption  capacity,  is  reached  at  which  no  further 
adsorption  occurs.  The  lower  portion  of  the  Langmuir  isotherm  is  also 
linear.  The  Langmuir  isotherm  is  defined  as 


_  Q cKjC, 

1  +  K,Cd 


(3.5) 


where 

Ki  =  Langmuir  adsorption  constant  (L  pg1) 

qc  =  maximum  amount  adsorbed  by  the  solid  (pg  g1). 

The  Langmuir  isotherm  was  based  on  different  assumptions,  one  of  which 
is  that  a  solid  surface  possesses  a  finite  number  of  adsorption  sites.  Once 
filled,  the  surface  will  no  longer  adsorb  solute  from  solution.  As  a  result,  a 
major  advantage  of  the  Langmuir  isotherm  over  linear  and  Freundlich 
isotherms  is  that  a  maximum  adsorption  capacity  is  incorporated  into  the 
formulation  of  the  model. 

3.3. 1.2  Computation  of  linear  equilibrium  partition  coefficients 

For  a  particular  contaminant,  the  partition  coefficient  in  water  depends  on 
the  nature  of  suspended  or  bed  sediment  solids,  and  key  geochemical 
parameters  of  the  water  and  sediment  layer.  In  CSM,  partition  coefficients 
are  user  specified  parameters  for  each  chemical  species.  Differences  in  the 
characteristics  of  water  column  and  sediment  layer  may  also  result  in 
different  partition  coefficients  of  the  contaminant.  Contaminant 
partitioning  for  POM  can  be  simulated  separately  from  inorganic  solids, 
rather  than  as  a  fraction  of  the  solids.  Values  for  the  partition  coefficients 
can  be  obtained  from  laboratory  experiments. 

For  organic  chemicals,  laboratory  studies  have  shown  that  the  partition 
coefficient  is  related  to  the  hydrophobicity  of  the  chemical  and  the  organic 
matter  content  of  the  sediment.  The  adsorption  of  hydrophobic  organic 
compounds  in  general,  is  considered  to  be  primarily  to  the  organic  matter  in 
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solids.  Studies  related  to  a  variety  of  organic  contaminants  have  generally 
shown  a  linear  log -log  correlation  between  Kow  and  Kd.  This  has  been 
demonstrated  by  several  researchers  (e.g.,  Karickhoff  et  al.  1979;  Seth  et  al. 
1999).  As  a  alternative,  the  equilibrium  partition  coefficients  for  DOC,  algae, 
POM,  and  solids  can  be  computed  in  CSM  from  a  user-defined  Kow 


Kdoo  ~adoc 


(3.6a) 


Kap=aap-Kow  (3-6b) 

K pom  =  apom '  Kow  (3.6c) 

Kpn  =  apn-Kow  (3.6d) 


where 


Kow  =  octanol-water  partition  coefficient  (unitless) 

(2 doc,  Clap,  Clpom,  Clpn  =  partition  correlation  coefficient  with  DOC,  algae,  POM, 
solids. 

The  user  may  input  Kow  rather  than  Kd.  In  general,  log  Kow  values  can  be 
found  from  the  literature.  The  above  equations  indicate  that  partition 
coefficients  and  the  rates  of  sorption  are  determined  not  only  by  the 
amount  of  organic  matter  but  also  by  the  amount  of  each  type  of  organic 
matter,  adoc,  aap,  ap0m,  a,m  vary  widely.  LaGrega  (1994)  reports  a  value  of 
(aP  om  —  0.63)  as  a  commonly  used  value  while  Seth  et  al.  (1999)  calculate  a 
value  of  (ap0m  =  0.35)  with  a  variation  in  a  by  a  factor  of  2.5  in  either 
direction.  To  use  this  relationship,  the  contaminant  must  be  nonionic, 
because  sorption  of  ionic  contaminants  are  affected  by  pH. 

3. 3. 1.3  Linear  equilibrium  partitioning 

Under  a  linear  equilibrium  partitioning,  the  distribution  among  the  phases 
is  controlled  by  the  partition  coefficients  Kd.  The  total  mass  of  contaminant 
in  each  phase  is  controlled  by  Kd  and  the  concentration  of  sorbent  present. 
For  a  contaminant  partitioning  to  DOC,  algae,  POM  and  multiple  solids 
phases,  the  fractions  of  each  phase  are  functions  of  the  partition  coefficients 
and  the  concentrations  of  DOC,  algae,  POM  and  solids.  Solids  should  be 
representative  of  fine-grained,  cohesive  solids  (e.g.,  clay,  silt,  organic 
matter).  Formulations  of  computing  the  linear  equilibrium  partitioning 
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fractions  of  the  contaminant  on  a  total  volume  basis  in  the  water  column  are 
written  as 


and 


J  d  ~  ^ 


f  _  C doc 

Jdi 


’  doc 


c. 


f  =^L-_ 

J  ap 


_  ^ pom 

pom 


f  =  £». 

J  pn 


106 

(3.7a) 

R 

KdocDOC 

(3.7b) 

R 

*§ 

•> 

a- 

(3.7c) 

R 

.  KpomPOM 

(3.7d) 

R 

_  Kpnmn 

(3.7e) 

R 

N 

d  +  KpomPOM  +  YJKpnmn 

(3.7  f) 

In  CSM,  equilibrium  partitioning  of  contaminants  is  handled  in  the  same 
manner  in  the  sediment  layer  as  in  the  water  column.  The  description  of 
water  column  partitioning  provided  above  applies  for  partitioning  in  the 
sediment  layer.  The  fractions  associated  with  dissolved  phase  in  pore 
water,  DOC,  and  sediment  solids  adsorbed  phases,  are  calculated  by 
considering  the  water  content  and  porosity.  Formulations  of  computing 
the  linear  equilibrium  partitioning  fractions  of  the  contaminant  in  the 
sediment  layer  are  written  as 


f  _£d2 
Jd2  ~  ^ 

'-"TO 


fa, 


C, 


c., 


A 


c 


pom.2 


pom  2 


c. 


c 

f  =  Pn2 

J  pn2 


106cp 

(3.8a) 

r2 

Kdoe2cpDOC2 

(3.8b) 

r2 

Rpom2  '  POM 2 

(3.8c) 

r2 

_  Kpn2mn2 

(3.8d) 

R 2 
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and 


R2  =  10 >  +  Kdoc2(p  ■  DOC2  +  Kpom2  ■  POM2  +  J2 Kpn2mn2  (3.8e) 

n—1 

The  symbols  and  definitions  used  in  equations  3.7a  -  3.8e  are  given  in 
Table  3. 


Table  3.  Symbols  and  definitions  used  in  equations  3.7a  -  3.8e. 


Symbol 

Definition 

Units 

Water  column 

fa 

Fraction  of  dissolved  phase  in  water 

- 

fdoc 

Fraction  of  DOC  adsorbed  phase  in  water 

- 

fap 

Fraction  of  algae  adsorbed  phase  in  water 

- 

fpom 

Fraction  of  POM  adsorbed  phase  in  water 

- 

fpn 

Fraction  of  solid  “n”  adsorbed  phase  in  water 

- 

Kdoc 

Equilibrium  partition  coefficient  for  DOC  in  water 

L  kg 1 

Kap 

Equilibrium  partition  coefficient  for  algae  in  water 

L  kg 1 

Kpom 

Equilibrium  partition  coefficient  for  POM  in  water 

L  kg 1 

Kpn 

Equilibrium  partition  coefficient  for  solid  “n”  in  water 

L  kg 1 

DOC 

Dissolved  organic  carbon  in  water 

mg  L1 

Apd 

Algal  biomass  (dry  weight)  in  water 

mg  L1 

POM 

Particulate  organic  matter  in  water 

mg  L1 

rrin 

Concentration  of  solid  “n”  in  water 

mg  L1 

Ct 

Total  concentration  in  water 

pgL1 

Sediment  layer 

fd2 

Fraction  of  dissolved  phase  in  sediment 

- 

fdoc2 

Fraction  of  DOC  adsorbed  phase  in  sediment 

- 

fpom2 

Fraction  of  POM  adsorbed  phase  in  sediment 

- 

fpn2 

Fraction  of  solid  “n”  adsorbed  phase  in  sediment 

- 

Kdoc2 

Equilibrium  partition  coefficient  for  sediment  DOC 

L  kg 1 

Kpom2 

Equilibrium  partition  coefficient  for  sediment  POM 

L  kg 1 

Kpn2 

Equilibrium  partition  coefficient  for  sediment  solid  “n” 

L  kg 1 

D0C2 

Sediment  dissolved  organic  carbon  in  pore  water 

mg  L1 

POM2 

Sediment  particulate  organic  matter 

mg  L1 

rrin2 

Concentration  of  sediment  solid  “n” 

mg  L1 

Ct2 

Total  concentration  in  sediment 

pgL1 
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DOC,  POM,  POM2,  APd  used  in  partitioning  equations  are  user  specified 
parameters  and  can  be  computed  from  the  NSM  modules  (Zhang  and 
Johnson  2016).  In  CSM,  equilibrium  partitioning  is  simulated  based  on 
the  total  concentration  of  the  contaminant.  At  equilibrium,  the 
concentration  of  chemical  in  any  phase  can  be  calculated  from  the  total 
concentration.  Therefore,  only  a  single  state  variable  representing  total 
concentration  is  used  for  each  contaminant.  Contaminant  concentrations 
for  each  phase  can  be  calculated  by 


Cd= £c'=£/'c‘ 

i  i 

(3.9a) 

cd0c=EcL=E-/14 

i  i 

(3.9b) 

c  =  E  c!  =  E  f  cl 

ap  /_ _ j  ap  /  j  J  ap  T 

i  i 

(3.9c) 

c  =E  c1  =Ef  cl 

pom  /  j  pom  /  j  J  pom  T 

i  i 

(3.9d) 

c,.=Ect.=EE/i4 

i  i  n= 1 

(3.9e) 

EE  /;c' 

Cp,.=10s^=^ - 

EEm. 

(3.90 

i  n=l 


Where  the  subscript  i  denotes  the  ith  chemical  species  if  the  ionization  is 
modeled  in  CSM  (i  =  1  to  5),  the  subscript  n  denotes  the  nth  solids  class  (n  = 
ito  N). 

In  the  sediment  layer,  concentrations  of  dissolved,  DOC  and  adsorbed 
phases  can  also  be  calculated  from  the  total  concentration  with  respect  to  a 
unit  volume  of  total  sediments  by 

cd2=Ec^=E/«Q2  (3-10a) 

i  i 
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Qoo2  =  EcL2-E/1242 

i  i 

(3.10b) 

cpom2=EcPom2  =  E/pom242 

i  i 

(3.10c) 

Cr,2=T.CU=T,T.fUCT2 

i  i  n= 1 

(3.10d) 

N 

EE/;2c'2 

c  —  1 03  !  n=1 

^pts2  ±U  N 

i  n—1 

(3.10e) 

The  above  Cd2  and  Cdoc2  are  computed  based  on  the  mass  of  dissolved 
phase  and  adsorbed  on  DOC  on  the  sediment  total  volume.  Their  pore 
water  concentrations  of  the  contaminant  are  expressed  as 

Cdp2  =  Cd2/ty  (3.11a) 

Cdxr2  =  CdxJi p  (3.11b) 

The  symbols  of  the  contaminant  concentrations  used  in  equations  3.9a- 
3.11b  are  defined  in  Table  4. 


Table  4.  Symbols  and  definitions  of  the  contaminant  concentrations  in  equations  3.9a  -  3.11b. 


Symbol 

Definition 

Units 

Water  column 

Co 

Concentration  of  dissolved  phase  in  water 

pg  L1 

Cdoc 

Concentration  of  DOC  adsorbed  phase  in  water 

Pg  L 1 

Cap 

Concentration  of  algae  adsorbed  phase  in  water 

Mg  L 1 

Concentration  of  POM  adsorbed  phase  in  water 

pg  L1 

Concentration  of  solid  “n”  adsorbed  phase  in  water 

pg  L1 

Cpt 

Total  concentration  of  solids  adsorbed  phase  in  water 

pg  L1 

Cpts 

Total  concentration  of  solids  adsorbed  on  solids  in  water 

Mgg1 

Sediment  layer 

Cd2 

Concentration  of  dissolved  phase  in  sediment 

pg  L1 

Cdoc2 

Concentration  of  DOC  adsorbed  phase  in  sediment 

pg  L1 

Cpom2 

Concentration  of  POM  adsorbed  phase  in  sediment 

Pg  L 1 
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Symbol 

Definition 

Units 

Cpn2 

Concentration  of  solid  “n”  adsorbed  phase  in  sediment 

pgL1 

Cpt2 

Total  concentration  of  solids  adsorbed  phase  in  sediment 

pg  L1 

Cpts2 

Total  concentration  of  solids  adsorbed  on  sediment  solids 

hgg1 

Cdp2 

Concentration  of  dissolved  phase  in  pore  water 

pg  L1 

Cdocp2 

Concentration  of  DOC  adsorbed  phase  in  pore  water 

pg  L1 

3.3. 1.4  Non-linear  equilibrium  partitioning 


When  the  concentration  of  the  contaminant  increases  beyond  a  certain 
level,  not  enough  sites  might  be  available  for  adsorption.  Therefore, 
competition  effects  will  lead  to  an  increase  in  non-linear  partitioning  on  the 
concentration  of  the  contaminant.  However,  there  is  no  evidence  that 
competition  effects  the  adsorption  of  the  substance  onto  the  DOC  in  the 
range  of  environmentally  relevant  DOC  concentrations  in  water  (Krop  et  al. 
2001).  It  is  assumed  that  the  partitioning  distribution  between  DOC  and 
water  always  follows  a  linear  isotherm  and  this  is  true  in  HgSM  as  well.  In 
subsequent  sections,  Kdoc  is  defined  with  the  assumption  of  linearity  of 
contaminant  partitioning  isotherms.  Following  the  Freundlich  equilibrium 
sorption  isotherm,  total  concentrations  of  the  contaminant  in  the  water 
column  and  sediment  layer  can  be  determined  by  the  following  mass 
balances. 


Water  column: 


C„  +10-6  KdocCdDOC 


+  1CT3  KfapCbdapApd  +  1CT3  K fpomCbfom  POM  +  1CT3  ^  KJi>nCjpnmn  =  CT 


(3.12) 


Sediment  layer: 


Cd2  +  10-6KHn^C^DOC. 


docT-'dl* 


+  10-3^om2  (Cd2/cpf om2  POM2  + 10-35>^2  (Cd2/cp)opn2  mn2  =  C. 


bpn  2 


(3.13) 


With  the  Langmuir  equilibrium  partitioning,  similar  mass  balance 
equations  of  total  concentrations  are  developed. 


Water  column: 
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Cd+10BKdocCdDOC 

|  10 "KlapCdqcapApd  |  10 -3KlpomCdqcpomPOM  |  »  KlpnCdqcnmn  g  (3-14) 
103+^/npC,  103+^pomC,  10  3+KIpnCd 

Sediment  layer: 

Ci2+^Kdoc2Ci2DOC2 

|  10 "Klpom2(Cjy)qcpom2POM2  |  ^10~3  Klpn2(Ci2/q>)qcn2mn2  ^  (3-15) 

io3+^pom2(cd2/cp)  ti  103+Klpn2{Cd2M 

where 

Kjpom  -  Freundlich  adsorption  constant  for  POM  in  water  (pg  g1)  (pg  L'1)43 
Kjpom?  -  Freundlich  adsorption  constant  for  sediment  POM  (pg  g1)  (pg  L'1)43 
Kfap  =  Freundlich  adsorption  constant  for  algae  in  water  (pg  g-1)  (pg  L1)4’ 
Kjpn  =  Freundlich  adsorption  constant  for  solid  “n”  in  water  (pg  g1) 

(pg  L -1)43 

Kfpnp  -  Freundlich  adsorption  constant  for  sediment  solid  “n”  (pg  g1) 

(pg  L-i)-b 

bpom  =  Freundlich  exponent  for  POM  in  water  (uniteless) 
bpom2  =  Freundlich  exponent  for  sediment  POM  (uniteless) 
bap  =  Freundlich  exponent  for  algae  in  water  (uniteless) 
bPn  =  Freundlich  exponent  for  solid  “n”  in  water  (uniteless) 
bPn2  =  Freundlich  exponent  for  sediment  solid  “n”  (uniteless) 

KiaP  =  Langmuir  adsorption  constant  for  algae  in  water  (L  pg1) 

Kipom  =  Langmuir  adsorption  constant  for  POM  in  water  (L  pg1) 

Kipoim  =  Langmuir  adsorption  constant  for  sediment  POM  (L  pg1) 

Kipn  -  Langmuir  adsorption  constant  for  solid  “n”  in  water  (L  pg1) 

Kipnp  =  Langmuir  adsorption  constant  for  sediment  solid  “n”  (L  pg1) 
qcap  -  maximum  amount  adsorbed  by  algae  in  water  (pg  g1) 
qcPom  -  maximum  amount  adsorbed  by  POM  in  water  (pg  g1) 
qcpcmp  -  maximum  amount  adsorbed  by  sediment  POM  (pg  g1) 
qcn  -  maximum  amount  adsorbed  by  solid  “n”  in  water  (pg  g-1) 
qcnp  -  maximum  amount  adsorbed  by  sediment  solid  “n”  (pg  g1). 

The  above  non-linear  equations  3.12  to  3.15  contain  unknown 
concentrations  (Cd  and  Cda)  and  are  solved  numerically.  Two  numerical 
methods  are  implemented  in  CSM:  1)  Newton-Raphson,  and  2)  Bisection. 

Both  methods  are  described  in  Chapra  and  Canale  (2006). 
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3.3.2  Non-equilibrium  partitioning 

The  equilibrium  partitioning  methods  discussed  above  assumes  that  all 
sorption  sites  have  equal  energy,  an  infinite  number  of  adsorption  sites, 
and  exposure  time.  The  problem  with  equilibrium  adsorption  is  that  these 
assumptions  are  generally  false.  Often,  there  is  an  absence  of  equilibrium 
due  to  insufficient  exposure  time,  metabolic  biotransformation,  and 
nonlinear  relationships  for  very  large  and/or  super  hydrophobic 
compounds  (Bertelsen  et  al.  1998).  For  example,  polychlorinated 
biphenyls  in  Lake  Ontario  exhibit  a  25-fold  non-equilibrium  (Cook  and 
Burkhard  1998).  If  the  adsorption  process  is  slow  or  irreversible,  dissolved 
and  adsorbed  concentrations  cannot  be  directly  calculated  from  the  total 
concentration  of  the  contaminant.  The  state  variables  for  a  contaminant 
include  the  concentrations  of  potential  phases. 

Under  non-equilibrium  partitioning,  the  effects  of  the  adsorption  and 
desorption  of  the  contaminant  between  water  and  solids  must  be 
computed  separately  in  the  transport  equations.  The  rate  of  adsorption  is  a 
function  of  the  adsorptive  capacity  of  the  solid  and  the  dissolved 
concentration,  while  desorption  rate  depends  only  on  the  adsorbed 
concentration.  Net  adsorption  and  desorption  rates  of  the  contaminant 
associated  with  algae,  POM,  and  solids  can  be  written  as 

^r  =  Kdap{lO-3qcapApd -Cap)cd -kdaapCap  (3.16a) 


dc 


dt 


POm  =  Kdpom  3  QcpomPOM  -  Cpom)Cd  —  kdapomCpom 


(3.16b) 


DC, 


dt 


Kdn  {10~3(lCnmn-Cpn)Cd  ~  Karfi p 


(3.16c) 


where 

kadn  =  adsorption  coefficient  for  solid  “n”  in  water  (L  pg1  d1) 
kadap  =  adsorption  coefficient  for  algae  in  water  (L  pg-1  d1) 
ICadpom  —  adsorption  coefficient  for  POM  in  water  (L  pg1  d1) 
kdan  =  desorption  rate  for  solid  “n”  in  water  (d1) 
kdaap  =  desorption  rate  for  algae  in  water  (d1) 

ICdapom  —  desorption  rate  for  POM  in  water  (d1). 
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As  a  result,  two  parameters,  kadn,  and  kdan  are  required  to  account  for  non¬ 
equilibrium  partitioning.  The  kinetic  rate  coefficients  are  related  to  the 
equilibrium  constants. 

The  net  adsorption  and  desorption  rates  of  the  contaminant  associated 
with  sediment  POM  and  solids  can  be  written  as 


~  kadpom2  (l0~3  qcpom2POM2  —  Cpom2)— —  kdapom2Cpom2 


8C 


pn2 


qj.  adn2\^-®  (?Cn2^n2  ^  pn2j  ^  ^dan2^ pn2 


(3.17b) 


where 

kadm  =  adsorption  coefficient  for  sediment  solid  “n”  (L  pg1  d1) 
kadPom2  =  adsorption  coefficient  for  sediment  POM  (L  pg1  d1) 
kdan2  —  desorption  rate  for  sediment  solid  “n”  (d1) 
lcdapom.2  =  desorption  rate  for  sediment  POM  (d1)- 

In  the  above  non-equilibrium  equations,  the  concentrations  of  dissolved 
and  adsorbed  phases  are  dependent  upon  each  other,  and  each  is  a  function 
of  time.  In  this  case  the  contaminant  is  simulated  with  more  than  one  state 
variable,  which  may  include  dissolved,  POM,  algae,  and  solids  adsorbed 
phases.  Kinetic  source  and  sink  terms  will  be  computed  separately  for  each 
phase  and  their  transport  equations  must  be  solved  simultaneously.  For 
simplicity,  non-equilibrium  partitioning  only  applies  for  non-ionic  chemical 
species.  In  other  words,  the  ionization  and  non-equilibrium  partitioning 
processes  cannot  be  activated  at  the  same  time  in  CSM. 

The  adsorption  coefficients  are  related  to  the  equilibrium  constants 
(equation  3.3).  The  desorption  rates  are  additional  parameters  for  the 
non-equilibrium  option  and  are  water  temperature  dependent.  The 
desorption  rates  are  corrected  using  the  MAE.  Activation  energy  values 
( Ea )  are  defined  for  each  solid  class. 

3.4  Contaminant  reactions  and  transformations 

In  addition  to  species  and  phase  partitioning,  contaminants  participate  in 
chemical  or  biological  reactions  in  the  aquatic  environment.  These 
processes  can  influence  the  ultimate  fate  of  contaminants.  Contaminant 
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degradation,  photolysis,  hydrolysis,  volatilization,  generalized  second-order 
reaction,  and  transformations  are  implemented  within  CSM.  The 
importance  of  each  of  these  processes  is  highly  dependent  on  both  the 
contaminant  and  its  environment.  Some  contaminants  undergo  a  complex 
set  of  reactions,  while  others  behave  in  a  more  simplified  manner.  The  CSM 
can  be  used  as  a  simpler  model  to  conduct  simulations  of  dye  tracers,  or 
first  order  degradation.  More  complex  simulations  of  the  contaminant  may 
employ  all  processes  included  in  CSM. 

3.4.1  Degradation 

Contaminants  may  degrade,  either  chemically  or  biochemically.  The  rate 
of  contaminant  degradation  in  the  environment  is  a  complex  function 
dependent  on  the  numerous  local  conditions.  Often  the  individual 
degradation  processes  are  not  well  known  or  cannot  be  quantified  exactly 
for  a  given  contaminant.  The  degradation  kinetics  of  contaminants  is  then 
described  as  the  n  order  kinetics  without  coupling  with  the  limiting 
substrate  (Schnoor  1996). 


The  degradation  process  is  simulated  for  the  water  column  and  underlying 
sediment  layer.  Different  degradation  rate  coefficients  {kid,  kiP )  can  be 
specified  for  the  water  column  and  sediment  layer.  This  degradation 
process  of  contaminants,  whether  they  are  dissolved  or  adsorbed  to  solids, 
is  simulated  using  the  same  equations.  In  CSM,  the  overall  degradation 
rate  of  the  contaminant  is  calculated  from  each  partitioning  phase  and 
each  chemical  species 


sdec=T ; 


K/T)fd  + 

doc  )  f  doc  +  Kap 


(T)faP  +KPojT)f‘c 


+XXn<w;„ 


Cl  (3.18) 


where 


w 

KdociT) 

KaP(T) 

KPoJT) 

Kpn{T) 


contaminant  degradation  (pg  L1  d1) 
degradation  rate  for  dissolved  species  in  water  (d1) 

degradation  rate  for  DOC  adsorbed  species  in  water  (d1) 
degradation  rate  for  algae  adsorbed  species  in  water  (d1) 
degradation  rate  for  POM  adsorbed  species  in  water  (d1) 
degradation  rate  for  solid  “n”  adsorbed  species  in  water  (d1). 
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First-order  degradation  rates  appear  more  practical  and  can  be  obtained 
from  experiments  under  field  conditions.  Since  the  chemical  natures  of 
aqueous  solutions  and  solid  environments  differ  greatly,  contaminant 
degradation  for  each  form  may  occur  at  very  different  rates.  The  CSM 
allows  for  specifying  different  degradation  rates  for  each  form.  The 
degradation  rates  are  temperature  dependent  parameters  and  are 
corrected  using  a  MAE. 

3.4.2  Hydrolysis 

Hydrolysis  is  a  reaction  in  which  cleavage  of  a  molecular  bond  of  the 
contaminant  and  formation  of  a  new  bond  with  either  the  hydrogen  or  the 
hydroxyl  component  of  a  water  molecule  occurs.  Contaminants  in  water 
may  react  with  positively  charged  hydronium  ions  [H+],  negatively 
charged  hydroxide  ions  [OH  ],  or  neutral  water  molecules.  Therefore  pH 
can  affect  water  ionization  and  hydrolysis  reactions. 

Hydrolysis  may  affect  both  the  dissolved  and  particulate  phases.  The 
WASP  model  allows  for  specifying  hydrolysis  for  any  partitioning  phases. 
However,  in  many  model  applications,  hydrolysis  and  photolysis  reactions 
are  considered  to  affect  only  the  dissolved  phase  (Schnoor  1996). 
Therefore,  the  hydrolysis  process  is  simulated  in  CSM  for  only  the 
dissolved  phase  in  the  water  column  and  underlying  sediment  layer. 
Different  hydrolysis  rate  coefficients  ( kha ,  khn,  khb )  can  be  specified  for  the 
water  column  and  sediment  layer.  In  CSM,  the  hydrolysis  rate  of  the 
contaminant  is  calculated  based  on  the  first  order  reaction  for  the  neutral, 
and  second  order  reaction  for  the  acidic  or  basic  ions. 

sw=E 


=  contaminant  hydrolysis  (pg  L1  d1) 

=  acid  hydrolysis  rate  constant  for  dissolved  species  in  water 
(L  mol1  d1) 

=  base  hydrolysis  rate  constant  for  dissolved  species  in  water 
(L  mol1  d1) 

=  neutral  hydrolysis  rate  for  dissolved  species  in  water  (d1) 


where 

Shyd 

w 

w 


[kla(T)[H+}  +  kihn(T)  +  kihb(T)[OH-])f^  + 

(Kadoc(T)[H+}  +  Kndoc(T)  +  Kmoc  (T)[OH})fic 
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k\adoc ( T )  =  acid  hydrolysis  rate  constant  for  DOC  adsorbed  species  in  water 
(L  mob1  d1) 

k'h hdoc(T)  =  base  hydrolysis  rate  constant  for  DOC  adsorbed  species  in  water 
(L  mob1  d1) 

Kndoc(T )  =  neutral  hydrolysis  rate  for  DOC  adsorbed  chemical  species  in 
water  (d1) 

[H+],  [OH  ]  =  molar  concentrations  of  hydrogen  and  hydroxide  ions  (mol  L1), 
respectively. 

The  pH  is  a  water  quality  parameter.  H+  and  OH"  can  be  determined  from 
the  pH  value 


[H+]  =  10pH  (3.20a) 

[OH  ]  =  10-pOH  =  io(14-pH)  (3.20b) 

Hydrolysis  is  one  of  the  most  important  mechanisms  in  the  environment 
for  the  breakdown  of  a  parent  chemical  and  these  rates  are  affected  by 
water  temperature.  In  CSM,  the  hydrolysis  rates  are  corrected  using  the 
AE.  An  activation  energy  value  (Ea)  of  75  KJ  mol1  (a  mid-range  value  for 
organic  chemicals)  is  used  as  the  default  value. 

3.4.3  Photolysis 

Contaminant  photolysis  is  assumed  to  occur  due  to  the  water  surface  being 
exposed  to  sun  light.  However,  the  relative  significance  of  photolysis  in 
relation  to  other  processes  may  vary.  Aquatic  photolysis  is  the  transforma¬ 
tion  or  degradation  of  a  contaminant  that  results  directly  from  the 
adsorption  of  light  energy.  The  amount  of  contaminant  loss  due  to 
photolysis  is  a  function  of  the  quantity  and  wavelength  distribution  of 
incident  light,  the  light  absorption  characteristics  of  the  contaminant,  and 
the  efficiency  at  which  absorbed  light  yields  a  chemical  reaction.  Only 
certain  chemical  species  of  the  contaminant  can  be  photolytically  altered 
and  CSM  allows  for  specifying  photolyzed  species.  The  photolysis  process  is 
simulated  in  CSM  for  only  the  dissolved  phase  of  the  contaminant  in  the 
water  column. 

A  first-order  rate  equation  has  been  used  to  estimate  aquatic  photolysis 
rate  (Thomann  and  Mueller  1987;  Chapra  1997).  This  equation  relates  the 
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observed  surface  photolysis  rate  under  laboratory  conditions,  and  such 
values  are  then  corrected  to  the  field  conditions  by  using  a  light 
attenuation  coefficient 


I  l  — e 

Spht  ~  1-33 

opht 


-Am 


(1-0.56  +^oo/ioc)4  (3.21) 


where 

Spht  =  contaminant  photolysis  (pg  L1  d1) 
h  =  water  depth  ( h  >  o)  (m) 

=  surface  photolysis  rate  for  dissolved  species  (d1) 

k'phtdac  =  surface  photolysis  rate  for  DOC  adsorbed  species  (d1) 

Io  =  solar  radiation  at  the  water  surface  (W  m~2) 

Iopht  =  light  intensity  when  kPht  is  measured  (W  nr2) 

Amax  =  maximium  light  extinction  coefficient  (nr1) 

and 

XIMB  =  aI.X  (3.22) 

where 

cct  =  correction  factor  for  light  attenuation  (l.i  to  1.6) 

A  =  light  extinction  coefficient  (nr1). 

The  light  extinction  coefficient  (A)  is  computed  in  NSM  modules  (Zhang 
and  Johnson  2016).  In  NSMs,  it  is  assumed  that  suspended  solids,  algae, 
and  organic  matter  affect  the  light  extinction.  The  surface  photolysis  rate 
of  a  contaminant  is  not  a  temperature  dependent. 

3.4.4  Volatilization 

Volatilization  is  the  movement  of  a  contaminant  from  the  bulk  water  phase 
of  a  water  body  across  the  water-air  interface  and  into  the  air.  This  process 
depends  on  the  contaminant  and  the  physical  properties  of  the  water  body 
and  the  atmosphere.  The  CSM  assumes  that  only  dissolved  contaminants 
can  be  transferred  across  the  water-air  interface,  and  adsorption  to  solids 
or  organic  carbon  reduces  volatilization.  The  overall  volatilization  rate  is 
determined  by  the  well-known  two-film  theory  (Whitman  1923).  The  two- 
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film  method  assumes  that  two  “stagnant  films”  are  bounded  on  either  side 
by  well  mixed  compartments.  The  transfer  coefficients  Kl  and  Kg  are  used 
to  quantify  the  exchange  of  contaminants  between  water  and  atmosphere. 
The  volatilization  transfer  velocity  is  then  derived  as  the  reciprocal  of  the 
total  resistance 


Vv(T) 


kl  kgkhI 


RT 


wk  ) 


(3.23) 


where 

Vv(T)  =  volatilization  velocity  across  the  air-water  interface  (m  d  ‘J 
Kl  =  mass  transfer  velocity  from  the  liquid  film  (m  d1) 

Kg  =  mass  transfer  velocity  from  the  gaseous  film  (m  d1) 

R  =  universal  gas  constant  (8.314  Pa  m3  mol-1  °K_1) 

Twk  =  water  temperature  (°K) 

Kh  =  Henry’s  Law  constant  (Pa  m3  mol1). 

The  value  of  vv(T),  depends  on  the  intensity  of  turbulence  in  a  water  body 
and  in  the  overlying  atmosphere.  Mackay  and  Leinonen  (1975)  have 
discussed  conditions  under  which  the  value  of  Vv(T)  is  primarily  determined 
by  the  intensity  of  turbulence  in  the  water.  The  Kh  relates  the  concentration 
of  a  contaminant  in  the  air  phase  to  its  concentration  in  the  water  phase, 
strongly  affecting  the  air-phase  resistance.  Depending  on  the  Kh  value,  the 
water  phase,  the  air  phase,  or  both  may  control  the  volatilization  rate.  A 
common  way  to  define  Henry’s  Law  constant  is  by  dividing  the  partial 
pressure  by  the  aqueous-phase  concentration  as  in  equation  3.23.  The 
Henry  volatility  can  also  be  expressed  as  the  dimensionless  ratio  between 
the  gas-phase  concentration  of  a  species  and  its  aqueous-phase 
concentration. 

As  the  Henry's  Law  coefficient  increases,  the  transfer  tends  to  be 
increasingly  influenced  by  the  intensity  of  turbulence  in  water.  As  the  Kh 
decreases,  the  air-water  transfer  coefficient  tends  to  be  increasingly 
influenced  by  the  intensity  of  atmospheric  turbulence.  For  small  values  of 
Kh,  transfer  velocity  through  the  air  side  of  the  interface  will  control  the 
overall  volatilization  rate  of  the  contaminant.  At  very  low  values  of  Kh 
(<  10-7),  the  volatilization  velocity  is  so  slow  and  therefore  can  be 
considered  unimportant  as  an  inter-media  transfer  mechanism.  The 


ERDC/EL  TR-16-8 


41 


volatilization  velocity  is  usually  of  relatively  less  magnitude  in  lakes  and 
reservoirs  than  in  rivers  and  streams. 


Two  transfer  coefficients  need  to  be  determined,  Kl  for  the  liquid  film  and 
Kg  for  the  gas  film  bordering  the  interface  between  water  and  atmosphere. 
These  coefficients  are  in  fact  mass  transfer  velocities.  Numerous  empirical 
relations  exist  that  describe  the  transfer  coefficients  as  functions  of  the 
wind  speed  and/or  the  water  flow  velocity  (Lyman  et  al.  1990).  The 
following  formulations  are  implemented  in  CSM  for  computing  the 
transfer  coefficients  of  Kl  and  Kg  (Mill  et  al.  1982). 


kl  =  K 


32 


MW 


0.25 


(3.24a) 


Kg  = 168 uw 


18 


MW 


0.25 


(3.24b) 


where 

ka  =  oxygen  reaeration  rate  (m  d1) 
MW  =  molecular  weight  (g  mol-1) 

Uw  =  wind  speed  (m  s-1)- 


Oxygen  reaeration  rate  ka  in  riverine  systems  can  be  calculated  using  the 
formulations  implemented  in  NSM  modules  (Zhang  and  Johnson  2016). 
Concentration  differences  serve  as  the  driving  force  for  the  water  layer 
diffusion.  In  CSM,  the  volatilization  rate  of  the  contaminant  is  calculated  as 


Svlt=-VJT) 


fjCT- 


KJR-T 


wk 


(3.25) 


where 


Svit  =  contaminant  volatilization  (pg  L1  d1) 

Co  =  air  concentration  (gaseous)  of  contaminant  (pg  L1). 

Note  that  only  the  un-ionized  form  (neutral  molecule)  of  the  dissolved 
contaminant  is  subject  to  volatilization  if  the  ionization  process  is 
simulated  in  CSM.  The  volatilization  velocity,  which  is  specified  by  the 
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user  or  calculated  from  equation  3.23,  is  corrected  with  the  local  water 
temperature  through  the  MAE. 


3.4.5  Generalized  second-order  reaction 


A  generalized  second-order  reaction  is  included  to  allow  the  user  to 
simulate  the  effect  of  processes  that  are  not  considered  in  the  above 
mechanisms.  The  reaction  depends  upon  a  rate  constant  and  an 
environmental  factor  which  may  be  taken  to  represent,  for  example,  some 
reducing  or  oxidizing  agent.  In  CSM,  the  generalized  second-order 
reaction  is  linear  with  respect  to  the  concentration  of  a  contaminant  and 
environmental  property 


s  =  V 

exa  /  j 


k'(T)f'+k'  (T)f‘  +k‘  (T)f  +k‘  (T)f  +Yk‘  (T)f' 

ed v  y  J  d  edoc  v  y  J  doc  eap v  y  J  ap  epom  v  y  J  pom  /  j  epn  v  y  y  pi 


CilE]  (3.26) 


where 


Sexa 

[ E 1 

ked(T) 

kedoc(T) 

keap(T) 

kepom(T) 

kepn(T) 


contaminant  generalized  second-order  reaction  (pg  L1  d1) 
concentration  of  environmental  property  driving  generalized 
second-order  reaction  in  water  (mg  L1) 
second  order  rate  for  dissolved  species  in  water  ([mg  L1]1  d1) 
second  order  rate  for  DOC  adsorbed  species  in  water  ([mg  L1]1  d1) 
second  order  rate  for  algae  adsorbed  species  in  water  ([mg  L1]-1  d1) 
second  order  rate  for  POM  adsorbed  species  in  water  ([mg  L-1]-1  d1) 
second  order  rate  for  solids  adsorbed  species  in  water  ([mg  L1]1  d1) 


In  CSM,  the  second-order  rates  are  temperature  dependent  and  are 
corrected  using  the  AE.  Activation  energy  values  ( Ea )  are  needed  for  each 
chemical  species. 

3.4.6  Transformations 

In  an  aquatic  environment,  chemical  reactions  can  result  in  the  removal  of 
primary  contaminants,  while  sometimes  also  creating  secondary 
contaminants,  (e.g.,  a  parent  chemical-daughter  product  sequence).  In  CSM 
the  contaminants  may  be  independent  or  they  are  allowed  to  transform  into 
other  species.  This  process  is  useful  to  estimate  the  persistence  of 
contaminants,  including  their  degradation  products.  Due  to  a  lack  of 
detailed  process  information,  all  transformation  processes  are  grouped 
together  in  one  model  formulation.  This  is  represented  as  a  lumped 
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transformation  yield  operating  on  the  pool  of  the  reactants  with  rate 
constants  that  can  vary  spatially.  When  two  or  more  contaminants  are 
modeled,  linked  transformations  that  convert  one  state  variable  into 
another  can  be  calculated  by  specifying  a  pathway  and  transformation 
efficiency  coefficient  for  each  reaction  process. 

j 


where 

Strm  =  total  contaminant  transformation  (pg  L1  d1) 

j  =  species  index  that  can  be  transformed  into  daughter  products 
(unitless) 

kj~>i(T)  =  transformation  reaction  rate  (d1) 

Yj->i  =  transformation  yield  coefficient  from  species  /  to  i  (g-i/g-j). 

In  a  simulation,  the  contaminant  can  undergo  several  chemical 
transformations.  Any  combinations  of  contaminant  transformations 
between  two  state  variables  can  be  modeled  in  CSM  for  the  water  column 
and  sediment  layer.  Transformation  processes  include  degradation, 
hydrolysis,  photolysis  and  generalized  second-order  reaction.  Yield  rate 
coefficients  can  be  applied  at  varying  strengths. 

3.5  Contaminant  air  deposition,  settling  and  sedimentation 

3.5.1  Air  deposition 

Contaminants  may  be  added  to  a  water  body  by  both  wet  and  dry 
atmospheric  deposition.  A  user  defined  load  is  used  to  describe  the 
atmospheric  deposition  in  CSM. 

satm  =103  AsL0  (3.28) 


where 

Satm  =  contaminant  air  deposition  (pg  L1  d1) 

As  =  surface  water  area  (nr2) 

Lo  =  areal  deposition  of  contaminant  (g  nr2  d1). 

Contaminant  concentrations  in  air  and  rain  can  be  obtained  from 
monitoring  and  then  converted  into  loads.  The  wet  deposition  load  is 
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calculated  by  multiplying  the  precipitation  amount  by  the  contaminant 
concentration  in  the  rain.  The  dry  deposition  load  is  the  product  of  the 
deposition  velocity  and  the  air  concentration  of  the  contaminant. 


3.5.2  Settling 

The  settling  and  erosion  of  sediments  and  associated  contaminants  are 
two  important  processes  in  water  quality  modeling.  When  adsorbed  on 
suspended  solids,  contaminants  can  be  carried  to  the  sediments  by  solids 
settling  from  the  water  column,  possibly  becoming  a  part  of  the  bed 
sediment.  The  settling  of  contaminants  is  coupled  to  the  settling  of 
suspended  solids.  In  CSM,  the  settling  rate  of  the  contaminant  (g  nr3  d1)  is 
calculated  as  the  production  of  the  settling  velocity  and  the  adsorbed 
phase  concentration  in  the  water  column 


s-=|e 


v  f1  +v  f 

/  j  snnJ  nn  '  sapJ  c 


sapJ  ap  T  ^ somJ  pom 
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nJ  no 


cl 


(3.29) 


where 

Sset  =  contaminant  settling  (pg  L1  d1) 

Vspn  =  solids  settling  velocity  (m  d1) 

Vsap  =  algae  settling  velocity  (m  d1) 

Vsom  =  POM  settling  velocity  (m  d1)- 

Settling  of  suspended  solids  results  in  adsorbed  contaminants  being 
subtracted  from  the  water  column.  Although  a  constant  sediment  layer  is 
assumed  in  CSM,  its  position  relative  to  the  sediment-water  interface 
changes  as  additional  material  is  deposited  on  the  bed.  Thus,  deposition  of 
solids  from  the  water  column  results  in  older  sediments  moving  further 
from  the  sediment -water  interface. 

3.5.3  Re-suspension 

The  sediment  erosion  and  re-suspension  processes  bring  solids  and 
attached  contaminants  from  the  bed  sediment  to  the  water  column.  The  re¬ 
suspension  results  in  older  sediments  moving  toward  the  sediment-water 
interface.  Because  erosion  and  re-suspension  rates  are  highly  variable  in 
space  and  time,  contaminant  fluxes  due  to  the  erosion  and  re-suspension  of 
solids  with  their  attached  contaminants  are  also  highly  variable  in  space  and 
time.  In  CSM,  the  chemical  composition  of  the  re-suspending  sediment  is 
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assumed  to  be  the  same  as  that  of  the  bed  sediment.  The  re-suspension  rate 
of  the  contaminant  (g  nr  3  d1)  is  calculated  based  on  the  rate  at  which 
sediment  particles  are  eroded  (resuspended) 


v-  ^ rpn  J~pn  2 


f" 

'-'T2 


(3.30) 


where 

Sres  =  contaminant  re-suspension  (pg  L1  d1) 
h2  =  thickness  of  active  sediment  layer  (m) 

Vrpn  =  re-suspension  velocity  (m  d1)- 

The  settling  and  re-suspension  processes  of  contaminants  are  assumed  to 
operate  on  all  particulate  fractions  of  contaminants.  Solids  adsorbed 
contaminant  mass  migrates  according  to  the  migration  of  the  solids 
classes.  The  amount  of  chemical  mass  associated  with  each  solids  class  is 
computed  from  partitioning  relationships,  and  the  mass  of  each  adsorbed 
chemical.  This  is  distributed  in  the  same  manner  as  the  solids,  (e.g., 
deposition,  re-suspension,  and  sediment  layer  redistribution). 


3.5.4  Sediment  burial 

Since  a  constant  sediment  layer  is  simulated  in  CSM,  the  burial  of 
contaminants  adsorbed  on  sediment  solids  results  from  net  settling  at  the 
sediment-water  interface.  The  burial  rate  of  sediment  contaminants  is 
calculated  as  the  production  of  the  burial  velocity  and  the  adsorbed  phase 
concentration  in  the  sediment  layer 


S‘=V?' 


£/pn  2  +  /p, 


pom2 


V  n=l 


cl 

'-'T  2 


(3.31) 


where 

Sb  =  contaminant  sediment  burial  (pg  L1  d1) 

Vb  =  sediment  burial  velocity(m  d1). 

3.5.5  Sediment-water  transfer 

Dissolved  contaminants  may  be  exchanged  between  sediment  and 
overlying  water  by  means  of  several  processes.  The  sediment-water 


ERDC/EL  TR-16-8 


46 


transfer  of  contaminants  is  primarily  due  to  sediment  erosion,  deposition, 
diffusion  and  bioturbation.  In  general,  they  occur  more  or  less 
simultaneously  and  there  are  several  interactions  among  them.  In  CSM, 
the  sediment-water  interaction  and  transfer  of  contaminants  are  modeled 
using  with  a  simplified  approach.  The  overall  mass  transfer  coefficient  is 
applied  to  calculate  a  diffusion  flux  proportional  to  a  concentration 
gradient  across  the  sediment -water  interface. 

^Et^AP-C’)  (3.32) 

n  2  i 


where 


Ssw  =  contaminant  sediment-water  transfer  (pg  L1  d1) 

Vm  =  sediment-water  transfer  velocity  (m  d1)- 

The  sediment-water  transfer  velocity  is  either  specified  as  a  constant  input 
parameter  or  computed  using  empirical  formulations.  A  theoretical  model 
based  on  the  well-known  resistance  in-series  concepts  for  chemical 
transport  between  phases  was  proposed  by  several  researchers  (e.g., 
Thibodeaux  et  al.  2001;  Erickson  et  al.  2005).  The  in-bed  process  on  the 
sediment  side  is  predominantly  particle  bio  diffusion.  Whereas,  that  on  the 
water  side  is  based  on  the  well-known  benthic  boundary  layer  mass- 
transport  coefficient  concept.  Desorption  occurs  from  a  particle-bound 
state  into  solution  at  the  interface.  Individually,  these  processes  have  a 
sound  theoretical  basis,  with  numerous  observational  data  sets  both  in 
laboratory  studies  and  in  the  field.  Figure  5.  illustrates  the  chemical 
release  pathway  from  the  sediment. 

Four  equations  are  implemented  in  CSM  for  computing  the  sediment- 
water  transfer  of  the  contaminant:  1)  Thibodeaux  et  al.  (2001),  2)  Di  Toro 
et  al.  (1981),  3)  Boyer  et  al.  (1994),  and  4)  Schink  and  Guinasso  (1977). 


The  first  equation  is  Thibodeaux  et  al.  (2001)  who  gives  the  following 
equation  to  calculate  the  mass  transfer  velocity  from  bed  sediment 


1 

5  +A 

^bPb^d2  Vms 


(3.33) 


where 
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Db  =  biodiffusion  coefficient  representing  particle  diffusivity  (m2  d1) 
pb  =  sediment  bulk  density  (g  cm  3) 

Vms  =  user-defined  sediment-water  transfer  velocity  (m  d1) 
z2  =  sediment  bioturbation  layer  thickness  (m). 


Figure  5.  Schematic  representation  for  soluble  chemical  release  from  bed 

sediment. 


In  this  equation,  the  concentration  gradient  across  the  sediment-water 
interface  is  affected  by  the  sediment  adsorption.  The  second  equation  is 
based  on  the  one  given  by  Di  Toro  et  al.  (1981).  The  sediment-water  mass- 
transfer  velocity  is  defined  as  a  function  of  the  molecular  weight  of  the 
constituent  and  the  sediment  porosity. 

vm  =  0.19cp  (MW  )2/3  (3.34) 

The  third  equation  is  determined  from  a  sediment  diffusion  coefficient 
where  zi  is  the  thickness  of  the  pore  water  diffusion  layer  (Boyer  et  al. 

1994) 


V 


m 


cp3Z) 


zi 


(3.35) 


where 
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Dm  =  molecular  diffusion  coefficient  (m2  d1) 
zi  =  pore  water  diffusion  layer  thickness  (=  o.oi  m). 

The  fourth  equation  is  given  by  Schink  and  Guinasso  (1977) 


v.„  =8.64x10' 


MDm/v 

24 


2/3 


(3.36) 


where 


u*  =  flow  shear  velocity  along  the  bed,  which  is  approximately  10 
percent  of  the  mean  velocity  of  flow  (m  s_1) 
v  =  kinematic  viscosity  of  water  (m2  d1)- 

3.6  Environmental  factors 

It  is  known  that  environmental  properties  such  as  pH,  water  temperature, 
microbial  population,  light  attenuation,  and  organic  and  inorganic  solids 
influence  the  the  reactions  and  transformations  of  the  contaminants 
modeled  in  the  aquatic  environment.  For  example,  water  temperature 
affects  most  reaction  coefficients  and  kinetic  rates.  Also,  pH  values  are 
important  to  ionization  and  hydrolysis  reactions,  their  chemical  speciation 
is  directly  affected  by  pH.  Adorption  is  associated  with  organic  and 
inorganic  solids.  These  environmental  properties  are  important  forcing 
functions  for  computing  the  transport  and  fate  of  contaminants.  These 
properties  often  vary  significantly  with  space  and  time.  Table  5  provides  a 
list  of  environmental  variables  and  affected  contaminant  processes  in  CSM. 


Table  5.  List  of  environmental  dependent  variables  and  affected  processes  in  CSM. 


Symbol 

Definition 

Units 

Affected  Process 

Tw 

Water  temperature 

°C 

All  in  water  column 

Tsed 

Sediment  temperature 

°C 

All  in  sediment 

lo 

Solar  radiation  at  the  water  surface 

W  nr2 

Photolysis 

A 

Light  extinction  coefficient 

nr1 

Photolysis 

pH 

pH 

- 

Ionization,  Hydrolysis 

DOC 

Dissolved  organic  carbon 

mg-C  L1 

Partitioing 

DOC2 

Sediment  dissolved  organic  carbon 

mg-C  L1 

Partitioing 

Apd 

Algal  biomass 

Partitioing 

POM 

Particulate  organic  matter 

Partitioing 

POM2 

Sediment  particulate  organic  matter 

Partitioing 
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Symbol 

Definition 

Units 

Affected  Process 

Lo 

Deposition  load 

g  m  2  yr1 

- 

Vspn 

Solids  settling  velocity 

m  d'1 

Settling 

Vrpn 

Solids  re-suspension  velocity 

m  d'1 

Re-suspension 

Vb 

Sediment  burial  velocity 

m  d'1 

Burial 

The  dependent  variables  listed  above  must  be  specified  either  from 
observed  data  or  generated  from  other  modules  before  the  CSM  can  be 
employed.  Most  of  them  can  be  computed  from  the  GC  module  as 
described  in  this  report  and  NSM  modules  as  described  in  Zhang  and 
Johnson  (2016).  For  example,  NSMI  (nutrient  simulation  module  I)  can 
be  used  to  compute  Apd,  POM,  POM2,  DOC,  A,  and  pH  listed  above.  The 
NSMI  module  does  not  explicitly  compute  pore  water  concentrations  of 
DOC  in  the  sediment  layer,  which  must  be  defined. 

The  three  most  important  sediment  processes  associated  with  water 
quality  modeling  are  settling,  re-suspension  or  “erosion,”  and  “burial”  of 
sediment.  Adsorbed  contaminants  can  settle  to  the  bottom  with  suspended 
solids  and  algae.  Contaminated  sediments  may  be  re-suspended  into  the 
water  column.  Their  fates  are  controlled  by  hydraulic  factors  (e.g.,  flow, 
sediment  transport),  and  depend  strongly  on  the  sediment  size 
fractionation  (e.g.,  clay,  silt,  sand  and  gravel).  Relatively  simple  sediment 
algorithms  for  calculating  settling,  re-suspension,  and  burial  of  sediment 
particles  have  been  incorporated  directly  into  the  GC  module  discussed  in 
the  previous  chapter.  The  GC  module  computes  water  column  and 
sediment  bed  concentrations  of  cohesive  and  non-cohesive  solids,  water 
column  settling  rates,  deposition  and  re-suspension  of  solids  between  the 
water  column  and  sediment  bed. 

3.7  Water  column  contaminant  source/sink  equations 

The  pathways  modeled  for  each  contaminant  are  depicted  in  Figure  2. 

Major  physical  processes  simulated  in  CSM  are  particulate  settling,  re¬ 
suspension,  and  sediment-water  diffusion  of  contaminants.  In  addition, 
CSM  has  options  for  simulating  the  following  kinetic  processes:  multi -phase 
partitioning,  ionization,  sorption,  degradation,  hydrolysis,  photolysis, 
volatilization,  and  transformations.  They  can  all  be  included  into  the  source 
and  sink  term  in  the  transport  and  mass  balance  equation.  For  each 
simulation  time  step,  the  CSM  computes  the  rate  changes  of  contaminant 
concentrations  due  to  internal  sources  and  sinks. 


ERDC/EL  TR-16-8 


50 


A  well-mixed  water  column  with  an  underlying  sediment  layer  is  assumed 
for  the  source  and  sink  equations.  Concentrations  in  these  equations, 
including  those  for  sediment  contaminants,  are  expressed  in  terms  of  mass 
per  unit  volume  of  water  plus  solids  (pg  L1).  The  method  employed  here  is 
numerically  equivalent  and  simpler  because  all  the  state  variables  are  in 
the  same  units,  and  no  correction  term  is  needed  in  the  transport 
equations.  Furthermore,  this  method  is  more  amenable  to  integration  with 
H&H  transport  models,  for  which  the  common  concentration  units  are 
mass  per  unit  volume. 

3.7.1  Equilibrium  partitioning 

In  CSM,  the  source  (+)  and  sink  (-)  terms  of  water  column  contaminant 
under  equilibrium  partitioning  implementation  are  computed  for  the  total 
concentration  and  are  listed  in  Table  6. 

3.7.2  Non-equilibrium  partitioning 

Under  non-equilibrium  partitioning  implementation,  the  concentrations 
of  dissolved  and  adsorbed  phases  of  the  contaminant  are  dependent  of 
each  other.  The  state  variables  for  a  contaminant  must  include 
concentrations  of  potential  phases.  The  source  (+)  and  sink  (-)  terms  of 
the  contaminant  are  computed  in  CSM  for  each  contaminant  dissolved  in 
water,  adsorbed  to  DOC,  algae,  POM  and  solids  and  are  listed  in  Table  7. 
The  first  step  in  modeling  a  contaminant  with  non-equilibrium 
partitioning  is  to  determine  the  concentration  of  the  contaminant  in  each 
phase.  The  transport  equation  describing  the  mass  balance  of  a 
contaminant  with  non-equilibrium  pardoning  must  also  be  developed  for 
each  individual  phase  and  solved  simultaneously. 


Table  6.  Water  column  total  contaminant  source/sink  terms  and  pathways  under  equilibrium  partitioning. 


Source/Sink  term* 

Pathway 

II 

E-h 

to 

tO 

+  io34A> 

Air  deposition 

--Y 

hr 

(  N  \ 

v  f  +v  r  +y v  r  ci 

sapJ  ap  somJ  pom  spnJ  pn  T 

K  n= 1  ) 

Settling 

Re-suspension 

+  V;Z((/«  +  fdoc2  It  ~  (fj  +  fL  h  ) 

h  i 

Water-sediment  transfer 
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Source/Sink  term* 

Pathway 

-I 

i 

rt=l  J 

Degradation 

i  \ 

[  ]+ (7’) + D/j + ^Woc  (r) 

l  {+Kbdoc(T)[OH\  ) 

/iW 

Hydrolysis 

-1.33-^l1  e  "T*  (1  - 0.56C, )X {k‘pht +  k‘phldjT)f-oc )c; 

^  Opht  Anax  '  ^  i 

Photolysis 

i  r  c  ^ 

--v„(7)  fdCT-—±f— 
h  ^  KH!RTk) 

Volatilization 

-I 

i 

Ka  mi,  +  Kaoc  C T)fL  +  Kap  C T)fap  +  Kpom  C T)fpom  +  £  k Jepn  (T)f‘n 

n= 1 

c;[i] 

Generalized  second-order 
reaction 

+Z^(r)K^cv 

j 

Total  transformations 

*  SSris  the  sum  of  internal  source  and  sink  terms  of  total  contaminant  in  water  (|jg  L1  d'1). 


Table  7.  Water  column  contaminant  partitioning  phase  source/sink  terms  and  pathways  under 

non-equilibrium  partitioning. 


Source/Sink  term* 

Pathway 

Dissolved  phase 

ss„  = 

+  io34A> 

Air  deposition 

+j^in(cd2/<p-cd) 

Water-sediment  transfer 

~kld(T)Cd 

Degradation 

~{kjT)\ 

H+\+khn(T)  +  khb{T)[0H-\fd 

Hydrolysis 

J  1  -^max  'h 

i  1  e  fi  n  Vr  fr'ir 

Photolysis 

1-5-5  x  *  7  V1  ^  L  )K phA1  d 

-*0pht  Anax  *  ^ 

h 

(  c  ^ 

c  0 

Volatilization 

-ked(T)Cd[E] 

Generalized  second-order  reaction 

•  >.('/'»;  ,x;; 

j 

Total  transformations 
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Source/Sink  term* 

Pathway 

- kadap P)Cd (l O'3 qcapApd  - Cap)+  kdaap ( T)Cap 
-  Kipom ( T)Cd  (l 0“3  qcpomPOM  -  Cpom )+  kiapom ( T)Cpom 

~  Z  Vadn ( T)Cd  (l 0“3 qcnmn  -  Cpn )-  kdan ( T)Cpn } 

n= 1 

Adsorption-desorption 

DOC  adsorbed  phase 

ss^  = 

+  Vm  if' doc  2  !<t>~  C doc  ) 

Water-sediment  transfer 

—  KdofT)Cdoc 

Degradation 

-ikhadjT){H 

+  J +kndOC  if)  +  khbdoc  (t)[0H- J) Cdoc 

Hydrolysis 

T  1  —  ^hiax'^ 

Photolysis 

J  v1  u  30L. L  )K  phtdoc\L  doc 

-*0pht  Anax*" 

~ ked0c(T)Cdoc  e\ 

Generalized  second-order  reaction 

Algae  adsorbed  phase 

II 

§• 

co 

co 

1  r- 

- v  C 

^  sa/?  ap 

Settling 

-K„(T)Cap 

Degradation 

-K,rmc,r[i 

;] 

Generalized  second-order  reaction 

+  fdap(T)Cd(\ 

0  QcapApd  —  Cap  )-  kdaap  t T)Cap 

Adsorption-desorption 

POM  adsorbed  phase 

sspom  — 

1  r’ 

- v  C 

som  pom 

Settling 

^1  /?om  (^0^-' pom 

Degradation 

-kPo„AT)cPo„ 

,[£] 

Generalized  second-order  reaction 

+  ^0jr)C,(lO-3VmPOM-C;;oJ-^o;„(r)C;,om 

Adsorption-desorption 

Solids  adsorbed  phase 

^  = 

1 

- v  C 

^  5/?«  pn 

Settling 

+  -v  C  , 

h  rn  pn  1 

Re-suspension 
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Source/Sink  term* 

Pathway 

-kipn(T)Cpn 

Degradation 

-KJt)C„[e} 

Generalized  second-order  reaction 

+  Kin (T)Cd (l  (T3 qcnmn  -  Cpn )-  kdan (T)Cpn 

Adsorption-desorption 

*  SSd  is  the  sum  of  internal  source  and  sink  terms  of  dissolved  phase  in  water  (pg  L 1  d1). 


SSdoc  is  the  sum  of  internal  source  and  sink  terms  of  DOC  adsorbed  phase  in  water  (pg  L1  d  r). 

SSap  is  the  sum  of  internal  source  and  sink  terms  of  algae  adsorbed  phase  in  water  (pg  L 1  d1). 

SSpom  is  the  sum  of  internal  source  and  sink  terms  of  POM  adsorbed  phase  in  water  (pg  L 1  d1). 

SSpn  is  the  sum  of  internal  source  and  sink  terms  of  solids  adsorbed  phase  in  water  (pg  L 1  d1). 

3.8  Sediment  contaminant  mass  balance  equations 

The  sediment  layer  is  envisioned  as  a  single,  well-mixed  layer  (Figure  l). 
The  sediment  layer  is  assumed  to  have  constant  thickness,  volume,  and 
porosity.  The  sediment  layer  gains  mass  through  deposition,  but  loses  it 
through  re-suspension  as  well  as  sediment  burial.  As  the  active  sediment 
layer  accepts  a  new  sediment  increment  on  the  top,  an  equal  amount  is 
removed  from  the  bottom  for  burial.  When  sediment  re-suspension 
occurs,  it  is  assumed  that  deeper  sediments  below  the  active  layer  are 
entrained  into  the  active  layer  and  mixed  with  the  remaining  sediments. 

Similar  to  the  water  column,  contaminant  kinetics  in  the  sediment  layer 
include  partitioning,  chemical  reactions,  and  transformations. 
Contaminants  in  the  sediment  layer  are  partitioned  into  those  that  are 
freely  dissolved  in  pore  waters  and  those  that  are  adsorbed  to  DOC, 
organic  and  inorganic  solids.  The  one  exception  is  the  assumption  that 
algae  do  not  exist.  The  chemical  reactions  and  transformations  in  the 
sediment  layer  are  largely  based  on  the  same  processes  applied  in  the 
water  column.  However,  contaminant  volatilization  and  aquatic  photolysis 
are  ignored  for  the  sediment  layer.  The  description  of  most  water  column 
algorithms  provided  above  applies  for  the  sediment  layer  as  well.  The 
differences  come  into  play  via  the  parameters  used  for  the  reactions. 

In  CSM,  the  transport  of  contaminants  in  the  sediment  layer  is  simulated 
basically  as  two  separate  procedures.  One  procedure  accounts  for  changes 
in  the  mass  of  contaminants  within  the  bed  as  affected  by  deposition,  re¬ 
suspension,  and  deep  burial.  This  procedure  involves  mostly  particulate 
chemicals,  which  are  added  through  settling,  and  lost  through  re¬ 
suspension  and  burial.  The  other  procedure  solves  for  pore  water  transport 
of  the  dissolved  contaminants  (water  and  DOC  adsorbed  phases).  The 
effects  of  both  procedures  are  combined  with  the  biochemical  reactions  to 
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determine  the  changes  of  total  concentration  in  the  sediment  layer.  The 
coupling  between  the  water  column  and  sediment  layer  is  performed 
explicitly  using  internal  source  and  sink  terms.  The  equation  presented  here 
assumes  a  well-mixed  sediment  layer  with  an  overlying  water  column. 

3.8.1  Equilibrium  partitioning 

Under  equilibrium  partitioning  implementation,  the  source  (+)  and  sink 
(-)  terms  of  the  sediment  contaminant  are  computed  in  CSM  for  the  total 
concentration  (per  unit  volume  of  sediment)  and  are  listed  in  Table  8. 


Table  8.  Sediment  total  contaminant  source/sink  terms  and  pathways  under  equilibrium  partitioning. 


Source/Sink  terms* 

Pathway 

CO 

Co 

to 

II 

+  -YI 
h2  y1 

(  N  \ 

v  f  +v  r  +Yv  f  ci 

sapJ  ap  somJ  pom  spnJ  pn  T 

V  n= 1  y 

Deposition 

— — Y 
h2r 

f  N  \ 

2  CT2 

V  «=i  J 

Re-suspension 

-T -'LV»(fU+itfU)CtT2 

"2  i  V  n= 1  ) 

Deep  burial 

-v-fL  (to + fLi  te + fL  Yt  ) 

^2  i 

Water-sediment  transfer 

-z 

i 

Kd2(T)f'd2  +  k[doc2(T)f’oc2  +  Kpom2{T)fpom2  +±kipn2(T)fU 

n= 1 

c* 

'-T  2 

Degradation 

-z- 

i 

f  {kL(T)[H+]+  kUr )  +  kUr)[0H-])fj2  +  1. 

\{kLoAT)[H+]+  k‘hndoc2(T)  +  Kbdoc2(T)[oH-y'doc2\  T1 

Hydrolysis 

-I 

i 

k‘ed2(T)fd2  +kiedoe2(T)floc2  +kiepom2(T)f‘pom2  +  fjkiepn2(T)f'pn2 

n= 1 

ct2  \e2  ] 

Generalized  second-order 
reaction 

+  Tk2j^(T)Y2J^Ci2 

j 

Total  transformations 

*  SSt2  is  the  sum  of  source  and  sink  terms  of  total  contaminant  concentration  in  sediment  (pg  L1  d'1). 


3.8.2  Non-equilibrium  partitioning 

With  non-equilibrium  partitioning  implementation,  the  source  (+)  and 
sink  (-)  terms  of  the  sediment  contaminant  are  computed  in  CSM  for  each 
contaminant  dissolved  in  pore  water,  adsorbed  to  DOC,  POM  and  solids 
(again  per  unit  volume  of  sediment)  and  are  listed  in  Table  9. 
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Table  9.  Sediment  contaminant  partitioning  phase  source/sink  terms  and  pathways  under  non¬ 
equilibrium  partitioning. 


Source/Sink  term* 

Pathway 

Dissolved  phase 

II 

CO 

Co 

V,n 

K 

f  c  N 

,  d 
\  <P  J 

Water-sediment  transfer 

kdl  (T)Cd2 

Degradation 

~  [hal  (T)[h+  ]+  khn2  ( T )  +  khb2  (T)[OH-]fdl 

Hydrolysis 

~  ke(n  (T)Cd2  \E2  ] 

Generalized  second-order 
reaction 

+Z^(r)^,.ci2 

j 

Total  transformations 

k adp 

N 

n= 1 

torn  2  V  /  , 

(l  ^  qcpom2POM2  Cpom2)+  kdapom  2  (T)C pom2 

f{^-3qcpn2mn2-Cpn2)-kdan2(.T)Cpn2 

Adsorption-desorption 

DOC  adsorbed  phase 

ssdoc2  = 

-v-t(cdoc2l<l>-cdoc) 

h2 

Water-sediment  transfer 

—  doc2 

Degradation 

-{khad0jT)[H^\ 

+  khndoc 2  (^)  +  ^7iMoc2  (r)[(9// 

Hydrolysis 

kedocl  (T)Cdoc 2  [^2  ] 

Generalized  second-order 
reaction 

POM  adsorbed  phase 

sspom2  = 

+i' 

f  c  ,  c  ) 

^ sap ^ ap  ^ som^ pom  ) 

Deposition 

1  r< 

h2b  pom2 

Deep  burial 

^1  poml  p0m2 

Degradation 

^ epoml  p0m2 

fe] 

Generalized  second-order 
reaction 

+kdpom2(T)<^L( 

10  qcpom2k>OM2  —  Cpgm2)—  kdapom2(T)Cpgm2 

Adsorption-desorption 
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Source/Sink  term* 

Pathway 

Solids  adsorbed  phase 

SSpn  2  ~ 

£  1  V  C 

spn  pn 

Deposition 

_  1  c 

^rpn^pnl 

Re-suspension 

1  r 

h2  Vb  pn2 

Deep  burial 

~  k\pnl(?)C pn2 

Degradation 

-k,AT)C„M 

Generalized  second-order 
reaction 

^adn2  (T)  — — —  (l  0  qcn2mn2  —  Cp„ 2  )_  & dan2  (P)^pn2 

Adsorption-desorption 

*  SSd2  is  the  sum  of  source  and  sink  terms  of  dissolved  phase  in  sediment  (pg  L1  d-1). 

SSdoc2  is  the  sum  of  source  and  sink  terms  of  DOC  adsorbed  phase  in  sediment  (pg  L 1  d1). 

SSp0m2  is  the  sum  of  source  and  sink  terms  of  POM  adsorbed  phase  in  sediment  (pg  L 1  d1). 

SSpn2  is  the  sum  of  source  and  sink  terms  of  solids  adsorbed  phase  in  sediment  (pg  L 1  d1). 

3.9  CSM  parameters 

The  CSM  operates  based  on  the  schematization  already  set  up  for  the  H&H 
model  and  the  flows  computed  by  the  H&H  model.  Table  to  provides  a  list 
of  the  contaminant-specific  input  parameters  needed  to  run  the  CSM  model 
simulation  when  all  options  are  turned  on.  Each  contaminant  may  exist  as  a 
neutral  compound  and  up  to  four  ionic  species.  The  neutral  and  ionic 
species  can  exist  in  multiple  phases.  CSM  input  parameters  can  be 
determined  from  field  and  laboratory  experiments,  or  estimated  using 
chemical  properties  and  from  the  literature.  Knowledge  of  a  contaminant’s 
chemical  form  and  of  the  chemical’s  ability  to  be  adsorbed,  to  be  degraded, 
or  reacted  is  necessary  if  the  fate  of  the  chemical  is  to  be  modeled 
accurately.  Literature  values  were  taken  from  a  compilation  of  several 
studies.  It  gives  a  range  of  values  for  each  of  the  parameters  (if  available).  It 
should  be  noted  that  these  values  are  based  on  a  limited  literature  review. 
Most  of  them  are  model  calibration  parameters  in  real  world  applications. 

This  table  will  be  repeated  for  each  contaminant  and  each  water  quality 
region,  allowing  the  user  to  define  the  different  values  for  input 
parameters.  In  CSM,  several  simulation  options  are  available  and  the 
required  input  parameters  depend  upon  their  selections.  The  complete 
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mass  balance  approach  involves  adding  a  sediment  layer  to  the 
computational  grid.  For  the  same  contaminants  and  processes,  different 
parameter  values  are  allowed  for  the  water  column  and  sediment  layer.  All 
contaminant  parameters  can  be  specified  as  either  a  uniform  constant,  a 
single  value  for  all  water  quality  cells,  or  a  spatially-varying  value  by  user- 
defined  water  quality  regions. 


Table  10.  List  of  CSM  input  parameters. 


Symbol 

Definition 

Typical 

values3 

Units 

Temperature 

Correction11 

Global 

MW 

Molecular  weight 

10  - 103 

g  mol1 

Csd 

Solubility 

n/a 

pg  L1 

Dm 

Molecular  diffusivity 

n/a 

msd-1 

Vm 

Sediment-water  mass  transfer  velocity 

n/a 

m  d-1 

h2 

Sediment  layer  thickness 

0.1  -  0.15 

m 

12 

Sediment  bioturbation  layer  thickness 

0.05  -  0.1 

m 

Vss 

Solids  settling  velocity 

- 

n/a 

m  d'1 

Vs  om 

Organic  matter  settling  velocity 

- 

n/a 

m  d'1 

Water 

column 

Sediment 

layer 

Ionization 

Kal 

Kal2 

Ionization  constant  for  anionic  species 

n/a 

unitless 

Ea 

n/a 

Ka2 

Ka22 

Ionization  constant  for  anionic  species 

n/a 

unitless 

Ea 

n/a 

Kbl 

Kbl2 

Ionization  constant  for  cationic  species 

n/a 

unitless 

Ea 

n/a 

Kb2 

Kb22 

Ionization  constant  for  cationic  species 

n/a 

unitless 

Ea 

n/a 

Equilibrium  partitioning 

Kdoc 

Kdoc2 

Equilibrium  partition  coefficient  for 

DOC 

n/a 

L  kg 1 

Kpn 

Kpn2 

Equilibrium  partition  coefficient  for 
solid  “n” 

10-1- 105 

L  kg 1 

foe 

foc2 

Fraction  of  organic  carbon  in  solids 

0.005  - 
0.5 

unitless 

b 

b2 

Freundlich  exponent 

n/a 

unitless 

Kf 

Kf2 

Freundlich  adsorption  constant 

n/a 

(Mg  g1)  (Mg 

l-iyb 

Ki 

K/2 

Langmuir  adsorption  constant 

n/a 

L  mg1 

Qcn 

Qcn2 

Adsorption  capacity  for  solid  “n” 

n/a 

Mg  g1 

Kow 

- 

Octanol-water  partitioning  coefficient 

n/a 

unitless 

Of  doc 

- 

Partition  correlation  coefficient  for  DOC 

0.4 

unitless 
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Symbol 

Definition 

Typical 

values3 

Units 

Temperature 

Correction*1 

(Xap 

- 

Partition  correlation  coefficient  for 
algae 

n/a 

unitless 

CXpom 

- 

Partition  correlation  coefficient  for 

POM 

n/a 

unitless 

Ofpn 

- 

Partition  correlation  coefficient  for 
solids 

n/a 

unitless 

Non-equilibrium  partitioning 

kadn 

kadn2 

Adsorption  coefficient 

n/a 

L  pg 1  d1 

kdan 

kdan2 

Desorption  rate 

n/a 

d'1 

Qcn 

Cfcn2 

Adsorption  capacity  for  solid  “n” 

n/a 

Mgg1 

Volatilization 

Vv(T) 

- 

Volatilization  velocity 

0.6-25 

m  d-1 

0 

n/a 

Kh 

- 

Henry’s  constant 

1- 107 

Pa  m3  moL1 

Cog 

Air  (gaseous)  concentration 

0  -  1.0 

pgL1 

Photolysis 

kpht 

- 

Aquatic  photolysis  rate  for  dissolved 
phase 

0-10 

d-1 

kphtdoc 

- 

Aquatic  photolysis  rate  for  DOC 
adsorbed  phase 

0-10 

d'1 

lopht 

- 

Light  intensity  when  kpht  is  measured 

n/a 

W  nr2 

0(1 

- 

1.2  -  1.6 

unitless 

Degradation 

n 

Degradation  order 

1.0 

unitless 

kid(T) 

kid2(T) 

Degradation  rate  for  dissolved  phase 

0-0.5 

d-1 

0 

n/a 

kidoc(T) 

kldoc2(T) 

Degradation  rate  for  DOC  adsorbed 
phase 

0-0.5 

d'1 

0 

n/a 

kiap(T) 

Degradation  rate  for  algae  adsorbed 
phase 

0-0.5 

d-1 

0 

n/a 

klpom(T) 

klpom2(T) 

Degradation  rate  for  POM  adsorbed 
phase 

0-0.5 

d'1 

0 

n/a 

kiP(T) 

kip2(T) 

Degradation  rate  for  solids  adsorbed 
phase 

n/a 

d-1 

0 

n/a 

Hydrolysis 

khb(T) 

khb2(T) 

Alkaline  hydrolysis  rate  of  dissolved 
phase 

h- 

O 

1 

o 

m3  mol1  d'1 

Ea 

75 
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Symbol 

Definition 

Typical 

values3 

Units 

Temperature 

Correction*1 

kfm(T) 

khnl(T) 

Neutral  hydrolysis  rate  of  dissolved 
phase 

0- 102 

d-1 

Ea 

75 

kha(T) 

kha2(T) 

Acid  hydrolysis  rate  of  dissolved  phase 

0- 107 

m3  mol 1  d'1 

Ea 

75 

khbdoc(T) 

khbdoc2(T) 

Alkaline  hydrolysis  rate  of  DOC 
adsorbed  phase 

0- 107 

m3  mot1  d'1 

Ea 

75 

khndoc(T) 

khndoc2(T) 

Neutral  hydrolysis  rate  of  DOC 
adsorbed  phase 

0- 102 

d-1 

Ea 

75 

khadoc(T) 

khadoc2(T) 

Acid  hydrolysis  rate  of  DOC  adsorbed 
phase 

r-~ 

O 

1 

O 

m3  mot1  d'1 

Ea 

75 

Generalized  second-order  reaction 

ked(T) 

ked2(T) 

Second-order  rate  for  dissolved  phase 

n/a 

(mg  L1)-1  d-1 

Ea 

n/a 

kedoc(T) 

kedoc2(T) 

Second-order  rate  for  DOC  adsorbed 
phase 

n/a 

(mg  L1)-1  d-1 

Ea 

n/a 

keap(T) 

Second-order  rate  for  algae  adsorbed 
phase 

n/a 

(mg  L1)-1  d-1 

Ea 

n/a 

kepom(T) 

kepom2(T) 

Second-order  rate  for  POM  adsorbed 
phase 

n/a 

(mg  L-1)-1  d-1 

Ea 

n/a 

keP(T) 

keP2(T) 

Second-order  rate  for  solids  adsorbed 
phase 

n/a 

(mg  L*1)'1  d'1 

Ea 

n/a 

Transformations 

Yji 

Transformation  yield  coefficient  from 
chemical  j  to  i 

1.0 

unitless 

a.  Wool  etal.  (2006). 

b.  Ea  units  is  KJ  mol1. 


3.10  CSM  outputs 

This  section  summarizes  the  CSM  outputs  from  a  model  simulation.  The 
fundamental  output  produced  by  CSM  consists  of  predicted  concentrations 
and  predicted  mass  fluxes  at  specified  locations  within  the  system.  Concen¬ 
trations  of  the  contaminant  computed  for  each  physical  compartment 
(water  column  and  sediment  layer)  include  total,  dissolved,  DOC  adsorbed, 
algae,  POM  and  solids  adsorbed  phases.  Table  n  lists  the  symbols  and 
definitions  of  concentrations  of  the  contaminant  computed  in  CSM. 

In  addition,  contaminant  pathway  fluxes  of  the  contaminant  are  internally 
computed  in  CSM  and  can  be  reported  in  the  model  outputs.  Table  12  lists 
the  symbols  and  definitions  of  the  pathway  fluxes  computed  in  CSM.  The 
CSM  allows  the  user  to  selectively  turn  on  and  off  each  pathway  flux  in  the 
model  outputs. 


ERDC/EL  TR-16-8 


61 


Table  11.  List  of  contaminant  concentrations  computed  in  CSM. 


Symbol 

Definition 

Units 

Water  column 

Ct 

Total  concentration  in  water 

PgL1 

Cion 

Ionic  concentration  in  water 

PgL1 

cd 

Concentration  of  dissolved  phase  in  water 

pgL1 

Cdoc 

Concentration  of  DOC  adsorbed  phase  in  water 

pgL1 

Cap 

Concentration  of  algae  adsorbed  phase  in  water 

pgL1 

Cpom 

Concentration  of  POM  adsorbed  phase  in  water 

pgL1 

Cpt 

Total  concentration  of  solids  adsorbed  phase  in  water 

pgL1 

Cpts 

Total  concentration  of  solids  adsorbed  phase  in  water 

Pgg1 

Sediment  layer 

CT2 

Total  concentration  in  sediment 

pgL1 

Cion2 

Ionic  concentration  in  water 

Cd2 

Concentration  of  dissolved  phase  in  pore  water 

Cdoc2 

Concentration  of  DOC  adsorbed  phase  in  pore  water 

Cpom2 

Concentration  of  POM  adsorbed  phase  in  sediment 

Cpt2 

Total  concentration  of  solids  adsorbed  phase  in  sediment 

pgL1 

Cpts2 

Total  concentration  of  solids  adsorbed  phase  in  sediment 
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Table  12.  List  of  contaminant  pathway  fluxes  computed  in  CSM. 


Symbol 

Definition 

Units 

Water  column 

Sdec 

Contaminant  degradation  in  water 

pgL-id-i 

Shyd 

Contaminant  hydrolysis  in  water 

pg  L1 d1 

Spht 

Contaminant  photolysis  in  water 

pgL-id-i 

Svlt 

Contaminant  volatilization  in  water 

pg  L1  d'1 

Sexa 

Contaminant  generalized  second-order  reaction  in  water 

pg  L1  d'1 

Strm 

Total  contaminant  transformations  in  water 

pg  L1  d'1 

Satm 

Contaminant  air  deposition 

Pg  L1 d1 

Sset 

Contaminant  settling 

Pg  L1 d1 

Sres 

Contaminant  re-suspension 

pg  L-1  d'1 

Ssw 

Contaminant  sediment-water  transfer 

pg  L-1  d'1 

Sediment  layer 

Sdec 

Sediment  contaminant  degradation 

pg  L1  d1 

Shyd 

Sediment  contaminant  hydrolysis 

pg  L1 d1 

Sexa 

Sediment  contaminant  generalized  second-order  reaction 

pg  L-1  d'1 

Strm 

Total  sediment  contaminant  transformations 

pg  L-1  d'1 

Sset 

Contaminant  deposition 

pg  L-1  d'1 

Sres 

Sediment  contaminant  re-suspension 

pg  L-1  d'1 

Ssw 

Contaminant  sediment-water  transfer 

Pg  L1 d1 

sb 

Sediment  contaminant  burial 

Pg  L1 d1 
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4  Mercury  Simulation  Module  (HgSM) 

This  chapter  describes  the  theory  and  mathematical  formulations 
implemented  in  the  Mercury  Simulation  Module  (HgSM).  The  HgSM  is 
capable  of  modeling  three  major  mercury  (Hg)  species  in  an  aquatic 
system;  elemental  mercury  (abbreviated  as  Hgo),  divalent  mercury 
(abbreviated  as  Hgll),  and  methylmercury  (referred  to  as  MeHg ),  and 
their  interactions  with  DOC,  algae,  organic,  and  inorganic  solids.  The 
HgSM  computes  internal  source  and  sink  terms  for  Hgo,  Hgll,  and  MeHg 
in  the  water  column  and  Hgll  and  MeHg  in  the  sediment  layer.  Due  to  the 
small  concentration  of  mercury  species  compared  to  the  major  water 
quality  constituents,  all  concentrations  of  mercury  species  in  HgSM  are 
expressed  in  terms  of  nanograms  per  volumetric  units  of  liter  (ng  Iv1)- 

4.1  Overview 

Mercury  is  introduced  into  the  environment  by  a  variety  of  natural 
processes  and  anthropogenic  activities.  Mercury  exists  in  several  different 
chemical  species  with  highly  differing  behaviors  and  toxicities.  Hgo  is  the 
most  common  form  of  mercury  found  in  the  atmosphere  (USEPA 1997). 
Besides  Hgo,  most  of  the  mercury  encountered  in  water,  soil,  sediments, 
and  biota  (all  environmental  media  except  in  the  atmosphere)  is  in  the  form 
of  inorganic  mercury  as  Hgll  and  organic  mercury,  particularly  MeHg 
(Mason  et  al.  1993, 1994).  Hgll  and  MeHg  can  exist  in  the  dissolved  phase 
as  well  as  in  the  colloidal  and  suspended  phases  within  aquatic  systems. 
These  mercury  species  often  undergo  complex  cycling  (Krabbenhoft  and 
Rickert  2003;  Harris  et  al.  2007;  EPRI  2006;  EPRI  2013).  One  of  the  major 
concerns  of  mercury  and  its  impact  to  the  environment  and  human  health  is 
its  transformation  to  MeHg,  a  highly  toxic  form  of  mercury.  A  crucial 
determinant  of  the  fate  of  MeHg  in  aquatic  systems  is  the  cycling  of  various 
species.  Cycling  consists  of  inputs  to,  and  outputs  from,  a  system,  as  well  as 
transport  and  transformation  of  the  mercury  species  within  the  system. 

The  HgSM  was  designed  to  simulate  three  species  (Hgo,  Hgll,  and  MeHg) 
as  model  state  variables.  Among  the  mercury  state  variables,  Hgll,  and 
MeHg  exist  in  both  the  water  column  and  the  sediment  layer  as  dissolved 
phase  and  associated  with  the  DOC,  POM  and  inorganic  solids.  Hgo  exists 
in  the  system  only  in  the  dissolved  phase.  Hgll  is  used  to  represent 
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inorganic  mercury  that  can  exist  as  Hg2+  complexed  to  chloride,  hydroxide 
and  sulfide  ions  (i.e.,  HgCl2,  Hg(OH)2,  HgS,  etc).  Mercury  species  can  be 
adsorbed  by  DOC,  algae,  POM  and  other  inorganic  solids.  Physical  and 
biochemical  processes  modeled  in  HgSM  include:  l)  adsorption  and 
desorption  of  mercury,  2)  volatilization,  3)  atmospheric  deposition, 

4)  diffusive  exchange  between  the  water  column  and  sediment  layer, 

5)  deposition  and  re-suspension,  6)  sediment  burial  of  adsorbed  mercury, 
and  7)  biogeochemical  transformations  among  three  species.  Similar  to  the 
CSM,  the  HgSM  only  computes  the  kinetics  of  mercury  state  variables  in 
aquatic  systems.  Figure  6.  provides  an  overview  of  the  model  representation 
of  the  mercury  species  and  cycling  processes  modeled  in  HgSM  (with  the 
exception  of  transport). 


Figure  6.  Conceptual  representation  of  mercury  speciation  and  processes  modeled  in  HgSM. 


Mercury  cycling  in  aquatic  systems  has  been  extensively  studied,  but  it  still 
represents  the  major  uncertainty  in  modeling  transport  and  fate  of  mercury 
species.  In  HgSM,  the  transformation  processes  of  Hgo,  Hgll,  and  MeHg 
and  their  reaction  rates  are  described  using  simple  kinetic  equations 
adopted  in  part  from  the  WASP  mercury  module  (Wool  et  al.  2006)  and 
D-MCM  (EPRI  2013).  The  WASP  mercury  module  was  developed  as  a 
companion  to  the  other  two  WASP  water  quality  submodels:  EUTRO  and 
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TOXI.  The  WASP  model  simulates  the  transport  and  fate  of  Hgo,  Hgll,  and 
MeHg  in  surface  water.  The  D-MCM  is  a  time  dependent,  mechanistic  mass 
balance  model  for  simulating  mercury  cycling  and  bioaccumulation.  The 
D-MCM  simulates  Hgo,  Hgll,  and  MeHg,  and  the  simplified  food  web  that 
includes  phytoplankton,  zooplankton,  benthos,  and  three  fish  species. 

4.2  Mercury  speciation  and  partitioning 

4.2.1  Mercury  speciation 

Mercury  speciation  plays  a  major  role  in  determining  the  cycling  and 
therefore,  the  effects  of  transport  and  fate  of  mercury  in  the  environment 
(Turner  1987;  Stein  et  al.  1996;  Jackson  1997;  Morel  et  al.  1998).  The  HgSM 
models  three  primary  mercury  species  of  environmental  concern:  Hgo, 
Hgll,  and  MeHg.  Gaseous  Hgo  is  the  dominant  form  in  the  atmosphere 
(Nater  and  Grigal  1992).  Hgo  has  an  approximate  atmospheric  residence 
time  of  6  months  to  1  year  (Fitzgerald  1989).  The  mercury  species  differ 
greatly  in  their  solubility.  Hgo  is  slightly  water  soluble  and  has  a  high 
Henry’s  Law  constant  (Schroeder  and  Munthe  1998).  Hgo  constitutes  very 
little  of  the  total  mercury  in  the  surface  water  but  may  provide  a  significant 
pathway  for  the  volatilization  of  mercury  from  surface  waters.  In  HgSM, 
Hgo  is  assumed  to  exist  in  the  dissolved  phase  and  only  in  the  water 
column. 

The  main  dissolved  mercury  species  in  aquatic  systems  is  Hg2+complexes 
with  various  inorganic  and  organic  ligands.  MeHg  is  the  methylated  form  of 
mercury  that  can  be  synthesized  in  aqueous  systems.  Hg2+  and  MeHg  in 
aquatic  environments  are  generally  not  free  ions,  but  complex  to  various 
inorganic  or  organic  ligands,  including  hydroxide,  chloride,  sulfides,  and 
dissolved  organic  matter  (Morel  et  al.  1998).  In  oxidizing  conditions,  the 
dominant  forms  are  Hg2+and  Hg22+  as  chloride  and  hydroxide  complexes. 
While  in  reducing  conditions,  the  dominant  forms  are  sulphur-mercury 
compounds  such  as  mono-  and  bi-sulphide  complexes  (HgS  (s),  HgS2H2, 
HgS2H-  and  HgS22)  (Puk  and  Weber  1994;  Morel  et  al.  1998).  In 
intermediate  conditions,  the  most  common  forms  are  alkyl  mercury 
compounds  (Rundgren  et  al.  1992).  The  HgSM  combines  all  Hg2+  species  of 
mercury  into  a  state  variable  identified  by  Hgll.  Hgll  is  modeled  for  the 
water  column  and  sediment  layer. 

Hgo  and  Hgll  are  the  major  chemical  forms  of  mercury  input  into  the 
environment  from  anthropogenic  or  natural  sources  (Stein  1996).  During 
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the  biogeochemical  cycling  of  mercury,  organic  mercury  species  can  be 
produced.  In  freshwater  systems  the  predominant  form  of  organic 
mercury  compound  is  methylmercury  (CH3Hg+),  which  is  represented  by 
MeHg  Although  MeHg  only  accounts  for  a  small  fraction  of  mercury  in  the 
environment,  it  is  the  species  of  most  concern  to  humans  due  to  its  high 
toxicity,  prevalent  existence,  and  capability  of  being  accumulated  and 
amplified  along  the  food  chain.  The  bioaccumulation  of  MeHg  in  the  food 
chain,  in  particular  in  aquatic  ecosystems,  has  caused  exposure  of  both 
humans  and  wildlife  to  MeHg,  posing  severe  health  risks  (Zillioux  et  al. 
1993;  Carroll  et  al.  2000;  Wang  et  al.  2004;  Evers  2005).  Therefore,  MeHg 
is  modeled  for  the  water  column  and  sediment  layer. 

4.2.2  Mercury  partitioning 

Mercury  species  may  be  present  under  different  phases  in  aquatic  systems: 

1)  elemental  Hgo  as  a  non-aqueous  liquid  phase,  2)  dissolved  in  the  aqueous 
phase  as  a  free  ion  (Hg2+)  or  complexed  with  inorganic  and/or  organic 
ligands,  3)  sorbed  on  inorganic  solids  and  organic  matter,  and  4)  in  solid 
(precipitated)  phase.  The  importance  of  adsorption/desorption  on  mercury 
bioavailability  has  been  confirmed  by  both  laboratory  and  field  studies.  It 
has  been  well  documented  that  DOC  and  solids  are  two  principal  carriers  of 
MeHg  and  Hgll  (Lyon  et  al.  1997).  DOC  is  an  effective  complex  for  ligand 
for  many  trace  metals,  including  mercury.  Hgll  and  MeHg  may  exist  in 
various  dissolved  and  particulate  forms  in  aquatic  systems.  In  HgSM,  all 
Hgo  is  assumed  in  the  dissolved  phase,  Hgll  and  MeHg  are  allowed  to  be 
partitioned  among  the  dissolved  phase  and  DOC,  algae,  and  solids  adsorbed 
phases.  Similar  to  the  CSM,  solids  can  be  assigned  to  mineral  abiotic  solids, 
detrital,  or  various  classes  or  size  categories  of  solids.  Adsorption  or  uptake 
of  Hgll  and  MeHg  with  algae  in  the  water  column  is  simulated  in  a  similar 
manner  as  is  done  for  contaminants  in  CSM.  Inside  the  algae  mercury 
kinetics  and  transformations  are  not  considered  in  HgSM.  Sediment  Hgll 
and  MeHg  are  allowed  to  be  partitioned  between  dissolved,  DOC,  algae, 
POM  and  solids  adsorbed  phases.  Under  an  equilibrium  partitioning,  all 
Hgll  and  MeHg  are  assumed  to  be  instantly  exchangeable  between  the 
dissolved  and  adsorbed  phases. 

It  is  recognized  that  adsorbed  Hgll  may  be  not  at  true  "equilibrium"  with 
the  dissolved  phase.  Alternatively,  sorption  of  Hgll  on  algae,  organic  and 
inorganic  solids  can  be  simulated  according  to  non-equilibrium 
partitioning  kinetics  described  in  chapter  3. 
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4.2.2. 1  Linear  equilibrium  partitioning 

In  HgSM,  partitioning  factions  of  Hgll  and  MeHg  in  the  water  column  and 
sediment  layer  are  identically  calculated  based  on  the  formulations 
developed  for  contaminants  (Chapter  3.3).  Table  13  summarizes  formula¬ 
tions  for  calculating  the  linear  equilibrium  partitioning  fractions  of  Hgll 
and  MeHg.  The  symbols  in  these  equations  are  defined  in  Table  14. 

Partition  coefficients  of  Hgll  and  MeHg  have  been  calculated  for  sediment 
or  soil  over  water.  Values  of  the  partition  coefficient  Kd  on  the  order  of  10- 
los  L  kg1  for  soil,  ios  L  kg1  for  sediment  and  105+  L  kg1  for  suspended 
material  are  typically  found  for  Hgll  and  MeHg  (USEPA 1997),  indicating 
a  strong  preference  for  Hgll  and  MeHg  to  remain  bound  to  soil,  bottom 
sediment,  or  suspended  matter  (increasing  affinity  in  that  order).  Harris  et 
al.  (1996)  reported  that  the  partitioning  of  Hgll  to  total  suspended  solids  is 
20-25%  higher  than  partitioning  to  abiotic  solids,  while  the  partitioning  of 
MeHg  to  total  suspended  solids  is  2.5  times  higher  than  partitioning  to 
abiotic  solids. 

For  the  benthic  sediment,  the  values  of  the  partition  coefficient  Kd  ranging 
from  1.6-104  to  9.9-105  L  kg1  with  median  values  between  5.4-104  and 
7.9-104  L  kg1  for  Hgll  and  6.5-102  to  l.i-ios  between  6.1-ios  and  9.0-103  L 
kg-1  for  MeHg  are  reported  by  Lyon  et  al.  (1997).  Harris  et  al.  (1996) 
reported  that  the  values  of  the  partition  coefficient  of  Hgll  to  benthic 
solids  in  four  lakes  range  from  3.0-104  to  105  L  kg1,  while  the  values  of  the 
partition  coefficient  of  MeHg  range  from  2.2-103  to  7.8-103  L  kg1. 
Generally,  the  distribution  coefficient  for  water  column  partitioning  is 
higher  than  for  sediment  pore  water.  Similarly,  the  distribution 
coefficients  of  MeHg  are  typically  several  folds  smaller  than  the 
corresponding  values  for  Hgll. 


Table  13.  Formulations  of  calculating  Hgll  and  MeHg  partitioning  fractions  in  HgSM. 


Inorganic  Mercury  (Hgll) 

Methylmercury  (MeHg) 

Water  column 

to6 

J d-Hgll  ~  „ 

RHgII 

to6 

J d— MeHg  j-, 

^ MeHg 

,  _  Kdoc_HgIIDOC 

J  doc-Hgll 

KHgII 

r  _  K  doc-MeHgDOC 

J  doc— MeHg 

^ MeHg 
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Inorganic  Mercury  (Hgll) 

Methylmercury  (MeHg) 

K  A 

r  _  ap—Hglr^pd 

J  ap-Hgll  ~  D 

RHgII 

K  A 

r  V  ap-MeHg  pd 

J  ap—MeHg  q 

MeHg 

,  _  Kpom_HgIIPOM 

J  pom-Hgll  tj 

RHgII 

r  _  K  pom-MeHgP OM 

J  pom-MeHg  q 

^ MeHg 

r  _  Kpn-Hglimn 

J  pn-Hgll  ~  D 

KHgII 

r  K  pn-MeHg^n 

J  pn—MeHg  r, 

^ MeHg 

RHgn  =  106  +Kdoc_HglIDOC  +Kap_HgIIApd  +K  pom_HglIPOM  +  j^K pn_HgIImn 

n= 1 

N 

RMeHg  =  10  doc_MeHgDOC  +K ap  MeHgApd  +K pom_MeHgPOM  +  ^ 

n= 1 

Sediment  layer 

io  V 

J d-HgII2  R 

K  Hgll  2 

_  ioV 

J  d-MeHg2  R 

^  MeHg  2 

r  K  doc-HgII2t/>DOC2 

J  doc-Hglll  ~  r, 

^  Hgll  2 

r  _  K doc-MeHgl^OC \ 

J  doc—MeHg2  R 

^  MeHg  2 

,  _  Kpom_Hgii2POM2 

r  _  K pom-MeHgd^OM 2 

J  pom—HgH2  q 

^  Hgll  2 

J  pom— MeHg  2  ^ 

^  MeHg  2 

r  ^  pn-Hgll  2mn2 

J  pn-Hgll  2  -  R 

^  Hgll  2 

r  _  Kpn-MeHg2mn2 

J  pn-MeHg2  R 

^ MeHg  2 

R  Hgll  2  ~  10  (j)+K  -  Hgll  2  </>DOC2  +Kpom-Hgll2POM2  +  pn-Hgll  2mn2 

n= 1 

N 

R-MeHg2  =  10  V  +K doc- MeHg2(!)DOC 2  pom-MeHg2^>^^  2  ^ ^ pn - MeHg 2mn2 

71  =  1 

Table  14.  Symbols  and  definitions  used  in  Table  13. 


Name 

Definition 

Units 

Water  column 

fd-Hgll 

Fraction  of  dissolved  Hgll  in  water 

- 

fdoc-Hgll 

Fraction  of  DOC  adsorbed  Hgll  in  water 

- 

fap-Hgll 

Fraction  of  algae  adsorbed  Hgll  in  water 

- 

f pom-Hgll 

Fraction  of  POM  adsorbed  Hgll  in  water 

- 

fpn-Hgll 

Fraction  of  solids  adsorbed  Hgll  in  water 

- 

Kdoc-Hgll 

Hgll  equilibrium  partition  coefficient  for  DOC  in  water 

L  kg 1 

Kap-Hgll 

Hgll  equilibrium  partition  coefficient  for  algae  in  water 

L  kg 1 
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Name 

Definition 

Units 

Kpom-Hgll 

Hgll  equilibrium  partition  coefficient  for  POM  in  water 

L  kg 1 

Kpn-Hgll 

Hgll  equilibrium  partition  coefficient  for  solids  in  water 

L  kg 1 

fd-MeHg 

Fraction  of  dissolved  MeHg  in  water 

- 

fdoc-MeHg 

Fraction  of  DOC  adsorbed  MeHg  in  water 

- 

fap-MeHg 

Fraction  of  algae  adsorbed  MeHg  in  water 

- 

fpom-MeHg 

Fraction  of  POM  adsorbed  MeHg  in  water 

- 

fpn-MeHg 

Fraction  of  solids  adsorbed  MeHg  in  water 

- 

Kdoc-MeHg 

MeHg  equilibrium  partition  coefficient  for  DOC  in  water 

L  kg 1 

Kap-  MeHg 

MeHg  equilibrium  partition  coefficient  for  algae  in  water 

L  kg 1 

Kpom-MeHg 

MeHg  equilibrium  partition  coefficient  for  POM  in  water 

L  kg 1 

Kpn-MeHg 

MeHg  equilibrium  partition  coefficient  for  solids  in  water 

L  kg 1 

Sediment  layer 

fd-Hgll2 

Fraction  of  dissolved  Hgll  in  sediment 

- 

fdoc-Hgl!2 

Fraction  of  DOC  adsorbed  Hgll  in  sediment 

- 

fpom-Hgll2 

Fraction  of  POM  adsorbed  Hgll  in  sediment 

- 

fpn-Hgll2 

Fraction  of  solids  adsorbed  Hgll  in  sediment 

- 

Kdoc-Hgll2 

Hgll  equilibrium  partition  coefficient  for  sediment  DOC 

L  kg 1 

Kpom-Hgll2 

Hgll  equilibrium  partition  coefficient  for  sediment  POM 

L  kg 1 

Kpn-Hgll2 

Hgll  equilibrium  partition  coefficient  for  sediment  solids 

L  kg 1 

fd-MeHg2 

Fraction  of  dissolved  MeHg  in  sediment 

- 

fdoc-MeHg2 

Fraction  of  DOC  adsorbed  MeHg  in  sediment 

- 

fpom-MeHg2 

Fraction  of  POM  adsorbed  MeHg  in  sediment 

- 

fpn-MeHg2 

Fraction  of  solids  adsorbed  MeHg  in  sediment 

- 

Kdoc-MeHg2 

MeHg  equilibrium  partition  coefficient  for  sediment  DOC 

L  kg 1 

Kpom-MeHg2 

MeHg  equilibrium  partition  coefficient  for  sediment  POM 

L  kg 1 

Kpn-MeHg2 

MeHg  equilibrium  partition  coefficient  for  sediment  solids 

L  kg 1 

Concentrations  of  Hgll  and  MeHg  in  each  phase  in  the  water  column  and 
sediment  layer  can  be  calculated  from  the  partitioning  factions  and  are 
listed  in  Table  15.  Dissolved  concentrations  of  Hgll  and  MeHg  in  the 
sediment  layer  are  the  mass  of  dissolved  and  DOC  adsorbed  phases 
relative  to  the  total  volume  of  solids  and  water.  The  porosity  corrected 
concentrations  in  pore  water  are  also  listed  in  Table  15.  The  symbols  in 
Table  15  are  defined  in  Table  16. 
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Table  15.  Formulations  of  calculating  Hgll  and  MeHg  partitioning  concentrations  in  HgSM. 


Inorganic  Mercury  (Hgll) 

Methylmercury  (MeHg) 

Water  column 

HgHd  =fd.miHgII 

MeHgd=fd_MeHgMeHg 

Hgndoc  =  fdoc-HgIIHgII 

MeHSdcc  =  fdoc-MeHgMeH8 

HgIIpt=fJfpn.HgnHgII 

n  =  1 

MeHSP,  =itfP„-M.HgMeHg 

n —\ 

N 

y.  f  pn-HgIlHgH 

—  1  O3  n= 1 

N 

Y^fpn-MeHgMeHg 

Mr  Her  —  1  O3  "=l 

pts  ^  N 

Z'«„ 

n= 1 

mengpls  -  iu  N 

77=1 

Sediment  layer 

VgHd2=fd_HgII2HgII2 

MeHgd2  =  fd_MeHg2MeHg2 

HgHdoc2  =  fdoc-HglllHglli 

MeHgdoc2  =  fdoc-MeHglMeHgl 

Hgnpt2=fJfPn-HgInHgn2 

n= 1 

N 

MeHgpl2  =Tjfpn-MeHg2MeHg2 

77=1 

N 

y  f pn- Hgll 2^ gll 2 

H^TT  —  1  03  "=1 

N 

y  f  pn- Me Hg 2 Hg 2 

MpHrr  —  1  f)3  n= 1 

11 811  pts  2  1U  Af 

77=1 

iviengpts2  iu  # 

i«,2 

77=1 

Pore  water 

tq 

t 

II 

£q 

=MeH§J2l(l) 

MeH8docp2=  MeHgdoci!  <!> 

Table  16.  Symbols  and  definitions  of  Hgll  and  MeHg  concentrations  used  in  Table  16. 


Symbol 

Definition 

Units 

Water  column 

Hgll 

Concentration  of  total  Hgll  in  water 

ng  L-1 

Hglld 

Concentration  of  dissolved  Hgll  in  water 

ng  L-1 

Hgll  doc 

Concentration  of  DOC  adsorbed  Hgll  in  water 

ng  L-1 

Hgllpt 

Total  concentration  of  solids  adsorbed  Hgll  in  water 

ng  L-1 

Hgllpts 

Total  concentration  of  solids  adsorbed  Hgll  in  water 

ngg1 

MeHg 

Concentration  of  total  MeHg  in  water 

ng  L-1 

MeHgd 

Concentration  of  dissolved  MeHg  in  water 

ng  L-1 

MQHgdoc 

Concentration  of  DOC  adsorbed  MeHg  in  water 

ng  L-1 

MeHgpt 

Total  concentration  of  solids  adsorbed  MeHg  in  water 

ng  L-1 
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Symbol 

Definition 

Units 

MeHgpts 

Total  concentration  of  solids  adsorbed  MeHg  in  water 

ngg1 

Sediment  layer 

Hgll  2 

Concentration  of  total  Hgll  in  sediment 

ng  L-1 

Hglldp2 

Concentration  of  dissolved  Hgll  in  pore  water 

ng  L-1 

Hglldocp2 

Concentration  of  DOC  adsorbed  Hgll  in  pore  water 

ng  L-1 

Hglipt2 

Total  concentration  of  solids  adsorbed  Hgll  in  sediment 

ng  L1 

Hgllpts2 

Total  concentration  of  solids  adsorbed  Hgll  in  sediment 

ngg1 

MeHg2 

Concentration  of  total  MeHg  on  in  sediment 

ng  L-1 

MeHgdp2 

Concentration  of  dissolved  MeHg  in  pore  water 

ng  L-1 

MeHgdocp2 

Concentration  of  DOC  adsorbed  MeHg  in  pore  water 

ng  L1 

MeHgPt2 

Total  concentration  of  solids  adsorbed  MeHg  in  sediment 

ng  L-1 

MeHgPts2 

Total  concentration  of  solids  adsorbed  MeHg  in  sediment 

ngg1 

4.2.2.2  Non-linear  equilibrium  partitioning 

In  HgSM,  two  of  the  most  popular  adsorption  isotherms  including  the 
Langmuir  and  Freundlich  isotherms  are  used  for  Hgll  and  MeHg.  As  with 
non-linear  equilibrium  partitioning  implemented  in  CSM,  concentrations 
of  Hgll  and  MeHg  for  each  phase  must  be  solved  numerically  from  the 
total  concentration  of  Hgll  and  MeHg.  With  the  Freundlich  equilibrium 
partitioning  implemented,  concentrations  of  total  Hgll  and  MeHg  in  the 
water  column  and  sediment  layer  are  determined  with  the  following  mass 
balances. 

Water  column: 

HgIId  + 10-6  Kdoc  HgII HgIId DOC  + 10  3bpom  Kfpom  Hgfr HgIIdbpom POM 

n  (4.1a) 

+  10  3hapKfop  HgllHgIl/apApd  +  10  3bp"J2KfPn  ^Hgll'rm,,  =  Hgll 


MeHgd  +10  eKt 


doc- MeHg 


MeHgdDOC  + 10  3bpom  K 


MeHg  dbpom  POM 


+  10 


-3  bap 


K 


fap—MeHg 


MeHgdbapApd  + 10 


-3  bpn 


fpom — MeHg  H  d 

n  (4.1b) 

bpn 


fpn—MeHg 


MeHgdmn  =  MeHg 


n= 1 


Sediment  layer: 


HgHd2  + 10  "Kdnc  Hgll2HgIId2L)OC2  +10  —K, 

N 

+  ^q-3 bpa^^R  .  (  tt„tt  l^\bPn2. 


>  {HgIId2/q>)P  POM 


(4.2a) 


fpn-Hgll  2 


(HgIId2/(p) P"  mn2  =  HgII2 
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MeHgrl2  + 10  6  K 


doc- MeHg  2 


MeHgd2DOC2  + 10 


-3bpom2 


K 


fpom—MeHg2 


( MeHgdJ<p ) 


bpom2 


POM „ 


(4.2b) 


+ 10 -3bpn2J2K^n_MeHg2  (MeHgd2/q>)  P"  mn2  =  MeHg 


where 

J^-fpom-HgII  =  Hgll  Freundlich  adsorption  constant  for  POM  in  water  (pg  g1) 
(pg  L-1)-15 

J£.fpom-HgIl2  =  Hgll  Freundlich  adsorption  constant  for  sediment  POM  (pg  g- 

9  (pg  L'9_b 

I^/ap-Hgii  =  Hgll  Freundlich  adsorption  constant  for  algae  in  water  (pg  g1) 
(pg  L-i)'b 

t^Jpn-Hgll  Hgll  Freundlich  adsorption  constant  for  solid  “n”  in  water  (pg 
g-1)  (pg  L-1)-13 

I£fpn-Hgll2  —  Hgll  Freundlich  adsorption  constant  for  sediment  solid  “n” 

(pg  g-1)  (pg  L-i)'b 

J^-fpom-MeHg  =  MeHg  Freundlich  adsorption  constant  for  POM  in  water  (pg  g- 

9  (pg  L-j)-b 

J^-fpom-MeHg2  =  MeHg  Freundlich  adsorption  constant  for  sediment  POM  (pg 
g-1)  (pg  L-i)-b 

J^fap-MeHg  =  MeHg  Freundlich  adsorption  constant  for  algae  in  water  (pg  g- 

9  (pg  L-i)-b 

I^fpn-MeHg  =  MeHg  Freundlich  adsorption  constant  for  solid  “n”  in  water 
(pg  g-1)  (pg  L-i)'b 

I^fpn-MeHg2  —  MeHg  Freundlich  adsorption  constant  for  sediment  solid  “n” 
(pg  g-1)  (pg  L-i)'b 

bpom  =  Freundlich  exponent  for  POM  in  water  (uniteless) 
bpom2  =  Freundlich  exponent  for  sediment  POM  (uniteless) 
bap  =  Freundlich  exponent  for  algae  in  water  (uniteless) 
bPn  =  Freundlich  exponent  for  solid  “n”  in  water  (uniteless) 
bPn2  =  Freundlich  exponent  for  sediment  solid  “n”  (uniteless). 


With  the  Langmuir  equilibrium  partitioning  implemented,  concentrations 
of  total  Hgll  and  MeHg  in  the  water  column  and  sediment  layer  are 
determined  with  the  following  mass  balance  equations 
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t^lap-Hgll  =  Hgll  Langmuir  adsorption  constant  for  algae  in  water  (L  pg-1) 
I^lpom-Hgll  =  Hgll  Langmuir  adsorption  constant  for  POM  in  water  (L  pg_1) 

Klpo  m  -Hgll2  =  Hgll  Langmuir  adsorption  constant  for  sediment  POM  (L  pg_1) 
■K-lpn-Hgll  =  Hgll  Langmuir  adsorption  constant  for  solid  “n”  in  water  (L  pg-1) 
Klpn-Hgll2  =  Hgll  Langmuir  adsorption  constant  for  sediment  solid  “n”  (L  pg_1) 
qcaP-Hgii  =  maximum  amount  of  Hgll  adsorbed  by  algae  in  water  (pg  g1) 
q cpom-Hgii  -  maximum  amount  of  Hgll  adsorbed  by  POM  in  water  (pg  g1) 
q cp0m-Hgii2  =  maximum  amount  of  Hgll  adsorbed  by  sediment  POM  (pg  g1) 
q cn-Hgii  =  maximum  amount  of  Hgll  adsorbed  by  solid  “n”  in  water  (pg  g1) 
qcn -Hgnz  =  maximum  amount  of  Hgll  adsorbed  by  sediment  solid  “n”  (pg  g1) 
I^lap-MeHg  =  MeHg  Langmuir  adsorption  constant  for  algae  in  water  (L  pg-1) 
I^lpom-MeHg  =  MeHg  Langmuir  adsorption  constant  for  POM  in  water  (L  pg-1) 


I^lpom-MeHg2  =  MeHg  Langmuir  adsorption  constant  for  sediment  POM  (L  pg1) 
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-K-lpn-MeHg  =  MeHg  Langmuir  adsorption  constant  for  solid  “n”  in  water  (L  pg-1) 
■K-lpn-MeHg2  =  MeHg  Langmuir  adsorption  constant  for  sediment  solid  “n”  (L  pg1) 
q cap-MeHg  =  maximum  amount  of  MeHg  adsorbed  by  algae  in  water  (pg  g1) 
q cpom-MeHg =  maximum  amount  of  MeHg  adsorbed  by  POM  in  water  (pg  g1) 
q cpom-MeH92  —  maximum  amount  of  MeHg  adsorbed  by  sediment  POM  (pg  g1) 
qcn-MeHg  =  maximum  amount  of  MeHg  adsorbed  by  solid  “n”  in  water  (pg  g1) 
qm-M,Hg2  =  maximum  amount  of  MeHg  adsorbed  by  sediment  solid  “n”  (pg  g1)- 


4. 2. 2. 3  Non-equilibrium  partitioning 

As  discussed  by  several  researchers  (e.g.  Hudson  et  al.  1994;  Yin  et  al. 
1997),  Hgll  sorption  onto  algae  and  solids  is  partially  instantaneous  and 
partially  rate  limited  (slow).  Hgll  species  are  strongly  adsorbed  by  soils 
and  sediments  and  are  desorbed  slowly.  Clay  minerals  adsorb  mercury 
maximally  at  pH  6.  In  acid  soils,  most  mercury  is  adsorbed  by  organic 
matter.  When  organic  matter  is  not  present,  mercury  becomes  relatively 
more  mobile  in  acid  soils.  Therfore  a  non-equilibrium  sorption  of  Hgll 
onto  algae,  POM  and  inorganic  solids  is  included  in  HgSM.  Rate  limited 
adsorption  and  desorption  equations  of  Hgll  with  algae,  POM  and 
suspended  solids  can  be  written  as 


dHgIIai 

dt 

dHgll 

pom 

dt 

dHgll 


=  10- -3kadap  (o  _h-A  pd  -  Hgll  )  Hgll  d  -  kdHgIIap  (4.5a) 


dt 


adap  \Ucap-HgII^pd  ap)±±y±±d  daap *  ap 

=  10  -3kadpom  ( q„gIIPOM  -  Hgll  pom)  Hgll  d  -  kdapomHgII pom  (4.5b) 
-  =  10 -3kadn  (q,.„  Hgllmn  -  HgIIpn  )HgIId  -  kdanHgIIpn  (4.5c) 


where 


kadn  =  Hgll  adsorption  coefficient  for  solid  “n”  in  water  (L  pg1  d1) 
kadap  =  Hgll  adsorption  coefficient  for  algae  in  water  (L  pg1  d1) 
kadpom  —  Hgll  adsorption  coefficient  for  POM  in  water  (L  pg1  d1) 
kdan  =  Hgll  desorption  rate  for  solid  “n”  in  water  (d1) 
kdaap  —  Hgll  desorption  rate  for  algae  in  water  (d1) 
kdapom  —  Hgll  desorption  rate  for  POM  in  water  (d1). 


Rate  limited  adsorption  and  desorption  equations  of  Hgll  with  POM  and 
solids  in  the  sediment  layer  can  be  written  as 
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dH9lIp'lm7  =10  3 k  (q  POM  —  Hgll  ^HgII‘ 

adpom2  \  ±cpom-HqII2  2  ; 


dt 


—  k.  Hgll  ,  (4.6a) 

aapom2  pom  2 


adpom.2  \~Lcpom-HgII2  2  z)  pom  2  )  dapom2  zf  pom2 

<P 


kadn2  \Qcn-HgII2mn2  Hgll pn 2  j  ^ 


-kdan2HgIIpn2  (4.6b) 


where 

kadn2  =  Hgll  adsorption  coefficient  for  sediment  solid  “n”  (L  pg1  d1) 
1Cadpom2  =  Hgll  adsorption  coefficient  for  sediment  POM  (L  pg_1  d_1) 
kdan2  =  Hgll  desorption  rate  for  sediment  solid  “n”  (d1) 
lCdapom2  =  Hgll  desorption  rate  for  sediment  POM  (d1)- 

4.3  Mercury  transformations 

The  oxidation,  reduction,  methylation  and  demethylation  reactions  of 
mercury  species  are  assumed  to  be  widespread  in  the  environment.  Some 
of  the  reaction  processes  have  not  been  well  understood  and  defined. 
Therefore  simplified  transformation  processes  among  the  mercury  state 
variables  are  modeled  in  HgSM.  Table  17  summarizes  the  key  reaction  and 
transformation  processes  modeled  in  HgSM. 


Table  17.  Mercury  species  and  their  transformations  modeled  in  HgSM. 


Species 

Water  column 

Sediment  layer 

HgO 

Hgll 

MeHg 

Hgll 

MeHg 

HgO 

Oxidation  (HgO- 
>Hgll) 

Hgll 

Photoreduction 
(Hgll— >HgO) 

Methylation 

(Hgll~>MeHg) 

Methylation 

(Hgll~>MeHg) 

MeHg 

Photoreduction 

(MeHg->HgO) 

Demethylation 

(MeHg~>Hgll) 

Demethylation 

(MeHg->Hgll) 

In  natural  waters,  the  transformations  of  dissolved  Hgll  and  MeHg 
include  mercury  complexes  with  DOC  (Nriagu  1979;  Bloom  et  al.  1991). 
Therefore,  the  user  can  assign  the  DOC  adsorbed  phase  that  is  available 
for  the  transformation  in  HgSM. 

4.3.1  Elemental  mercury  oxidation 

Hgo  can  be  relatively  quickly  oxidized  by  different  reagents,  mainly  oxygen 
but  also  nitrates,  nitrites,  iron  hydroxides,  iron  phosphates,  sulphates, 
sulphur,  and  carbon  dioxide  (Stein  et  al.  1996).  In  HgSM,  oxidation  from 
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Hgo  to  Hgll  is  modeled  as  a  first-order  process.  The  first-order  oxidation 
rate  is  calculated  as  a  function  of  Hgo  concentration  in  the  water  column 

Hgo  — >  Hgll  —  k12(T)-  HgO  (4.7) 


where 

Hgo—>HgII  =  Hgo  oxidation  yield  into  Hgll  in  water  (ng  L1  d1) 
ki2(T)  =  Hgo  oxidation  rate  in  water  (d1). 

The  Hgo  oxidation  rate  (ki2)  are  temperature  dependent  and  corrected 
using  the  AE.  The  oxidation  yield,  Hgo—>HgII,  is  subtracted  from  the 
Hgo  as  a  sink  (-)  and  added  to  the  Hgll  as  a  source  (+)  with  a  specified 
yield  coefficient  (Fa)  of  i.o. 

4.3.2  Inorganic  mercury  reduction 

Hgo  may  be  formed  through  the  reduction  of  Hgll  or  the  demethylation  of 
MeHg  (Allard  and  Arsenie  1991).  On  the  basis  of  experimental  data  some 
authors  indicate  the  biological  processes  as  the  most  important,  while 
others  suggest  that  photoreduction  reactions  play  the  dominant  role 
(Costa  and  Liss  2000).  Photoreduction  of  Hgll  to  Hgo  is  light  dependent; 
thus,  reduction  of  Hgll  is  expected  to  peak  at  midday.  In  HgSM,  Hgll 
reduction  is  assumed  to  be  driven  by  sunlight,  and  the  surface  reduction 
rate  is  attenuated  through  the  water  column  using  a  specified  light 
extinction  coefficient.  Photoreduction  from  Hgll  to  Hgo  in  the  water 
column  is  modeled  as  a  photolysis  process  as  with  the  photoreduction 
kinetics  used  in  WASP 

J  A  ~Xmaxh 

Hgll  -  Hgo  =  1.33  0  -  (1-0.56 C L)(kd21fd_HgII  +kdoc2Jioc_HgII)HgII  (4.8) 

1  , .  A  tl 

opht  max 


where 

HgII—>Hgo  =  Hgll  photoreduction  yield  into  Hgo  in  water  (ng  L1  d1) 
kd2i  =  dissolved  Hgll  photoreduction  rate  in  water  (d1) 
kdoc2i  =  DOC  adsorbed  Hgll  photoreduction  rate  in  water  (d1)- 

The  surface  photoreduction  rates  ( kd2i ,  kdoc2i)  are  not  temperature 
dependent.  The  photoreduction  yield,  HgII—>Hgo,  is  subtracted  from  the 
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Hgll  as  a  sink  (-)  and  added  to  the  Hgo  as  a  source  (+)  with  a  specified 
yield  coefficient  ( Y21 )  of  1.0. 

4.3.3  Mercury  methylation 

Several  studies  indicated  that  Hgll  methylation  is  the  dominant  source  of 
MeHg  in  aquatic  systems  (USEPA 1997;  EPRI  2013).  Hgll  can  be 
methylated  to  MeHg  in  water  phase  through  biotic  (microbially  mediated) 
pathways  or  abiotic  pathways  (photo-mediated  or  non  photo-mediated 
chemical  methylation)  (Stein  1996;  Ullrich  2001).  Methylation  takes  place 
mainly  in  surface  layers  of  sediments  and  in  soils,  but  it  can  also  take  place 
in  water  column  although  more  weakly  (Regnell  et  al.  1996).  The 
mechanisms  of  Hgll  methylation  in  water  column  have  not  been  fully 
recognized  yet,  but  it  is  assumed  that  they  involve  microorganisms, 
similarly  as  in  bed  sediments  (Morel  et  al.  1998). 

Hgll  must  be  in  solution  or  easily  transferrable  form  for  methylation 
(Benoit  et  al.  2003;  Marvin-DiPasquale  et  al.  2009).  Distribution  of 
mercury  species  between  the  solid  and  aqueous  phase  significantly  affects 
the  bioavailability  of  the  Hg  species.  DOC  interacts  strongly  with  dissolved 
Hgll  and  likely  affects  methylation  rate  and  MeHg  production.  Gilmour 
and  Henry  (1991)  reported  that  increased  DOC  in  the  water  column  may 
increase  ligand  formation  between  DOC  and  dissolved  Hgll,  making  it 
unavailable  for  microbial  methylation.  In  HgSM,  methylation  from  Hgll  to 
MeHg  in  the  water  column  is  modeled  as  a  first-order  kinetic  process 
based  on  the  algorithm  used  in  WASP 

Hgll  MeHg  (kd23(T)fd_HgII  +  kdoc23(T)fdoc_HgII)HgII  (4.9) 


where 

HgII—>MeHg  =  Hgll  methylation  yield  in  water  (ng  L1  d1) 

kd23(T)  =  dissolved  Hgll  methylation  rate  in  water  (d1) 
kdoc23(T)  =  DOC  sorbed  Hgll  methylation  rate  in  water  (d1)- 

The  water  column  Hgll  methylation  rates  ( kd23 ,  kdoc23)  are  temperature 
dependent  and  corrected  using  the  AE. 

In  the  sediment  layer  sulfate  reducing  bacteria  (SRB)  are  the  most 
important  for  the  biological  methylation  of  Hgll.  Sulfur,  organic  carbon, 
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sediment  structure,  and  composition  in  the  sediment  layer  all  affect  MeHg 
production  by  changing  the  amount  of  bioavailable  Hgll  and  by 
stimulating  SRB  activity  (Langer  et  al.  2001;  Mason  et  al.  2005;  Lin  et  al. 
2012). 

The  Hgll  methylation  rate  is  linked  to  sulfate  reduction.  Strong 
relationships  were  observed  between  sulfate  reducing  rates  and  Hgll 
methylation  rates  in  sediments  (Benoit  et  al.  1999;  King  et  al.  1999;  King  et 
al.  2001).  Langer  et  al.  (2001)  shows  the  relationship  between  Hgll 
methylation  and  sulfate  in  sediments  (Figure  7. ). 

As  shown  in  Figure  7,  the  Hgll  methylation  rate  in  high-sulfate  estuarine 
sediments  is  far  lower  than  the  rate  in  low-sulfate  freshwater  sediments. 

In  fresh  water  sediments,  methylation  has  been  observed  to  be  optimal 
when  sulfate  concentrations  range  between  2  and  too  pM.  Above  this 
favorable  range  of  sulfate  concentrations,  the  production  of  sulfide 
through  sulfate  reduction  would  inhibit  mercury  methylation.  Below  this 
sulfate  concentration  range,  sulfate  availability  would  limit  sulfate 
reduction  and  mercury  methylation.  The  amount  of  MeHg  increased  in 
proportion  to  the  concentration  of  sulfate,  then  the  MeHg  decreased. 
Sulfate  reduction  removes  the  sulfate  and  ultimately  produces  sulfide.  The 
overall  reaction  equation  is  SO;2  +  2 CH20  ->  S  2  +  2 C02  +  2 H20  ■ 

Figure  7.  Relationship  between  Hgll  methylation  and  sulfate  in 
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The  reaction  of  Hgll  with  sulfide  to  produce  insoluble  HgS  decreases  the 
availability  of  Hgll  for  methylation  (Gilmour  et  al.  1998).  In  HgSM,  the 
sediment  Hgll  methylation  rate  is  simply  calculated  as  a  function  of 
sulfate  and  dissolved  Hgll  concentration  based  on  the  algorithm  used  in 
WARMF  (EPRI  2006) 


{Hgll  — v  MeHg)2  =  kso42(T) 


SO  4, 


K  „ 


SO  4„ 


$04 2  TmsoAf Hg n  2 II g  II 2 


(4.10) 


where 

(HgII—>MeHg)2  =  Hgll  methylation  yield  in  sediment  (ng  L1  d1) 
kso42(T)  =  sediment  sulfate  reduction  rate  (d1) 

Vmso4  =  ratio  of  sediment  methylation  rate  and  sulfate  reduction  rate  (L 
mg-1) 

Kso4  =  half-saturation  constant  for  the  effect  of  sulfate  on  methylation 
(mg-02  L-1) 

SO 42  =  sediment  pore  water  sulfate  concentration  (mg-02  L1). 

The  above  equation  relates  methylation  to  the  sediment  sulfate  reduction 
rate  but  requires  that  the  sulfate  must  be  modeled  in  the  system.  A  half¬ 
saturation  constant  is  used  to  adjust  the  sulfate  effect  on  the  methylation 
rate.  This  relationship  results  in  methylation  rate  dependencies  which  are 
nearly  linear  when  sulfate  concentrations  are  much  less  than  KSo4. 
Conversely,  the  methylation  rate  does  not  vary  with  sulfate  if  the  sulfate  is 
much  greater  than  KSo4.  However,  the  above  algorithm  does  not  adequately 
capture  the  complex  interactions  between  MeHg  production  and  the  sulfur 
cycle.  The  inhibitory  influence  of  H2S  on  Hgll  methylation  is  not 
otherwise  captured  in  the  above  equation  4.10. 

The  sediment  sulfate  reduction  rate  is  temperature  dependent.  A  Q10 
equation  is  used  to  correct  sulfate  reduction  rate  according  to  local 
temperature  (Wool  et  al.  2006;  EPRI  2013).  The  methylation  yield,  Hgll— 
>MeHg,  is  subtracted  from  the  Hgll  as  a  sink  (-)  and  added  to  the  MeHg 
as  a  source  (+)  with  a  specified  yield  coefficient  (F23)  of  1.07. 

4.3.4  Methylmercury  reduction 

Photolysis  of  MeHg  has  been  shown  to  occur  in  water  (Callahan  et  al. 

1979).  Photoreduction  from  MeHg  to  Hgo  is  simulated  as  a  light 
dependent  process  in  HgSM.  This  process  occurs  primarily  in  the  water 
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surface  due  to  light  attenuation.  Photoreduction  of  MeHg  can  also  be 
affected  by  the  speciation  of  MeHg  in  the  water  column.  A  previous  study 
reported  that  sunlight  induced  MeHg  photoreduction  could  not  occur 
when  the  MeHgCl,  MeHgOH  or  MeHg  ion  was  the  dominant  species  of 
MeHg  in  water,  and  while  phenyl  and  sulfur  bonded  MeHg  species  could 
be  decomposed  (Baughman  et  al.  1973). 

In  HgSM,  the  MeHg  reduction  rate  is  calculated  based  on  the  light 
attenuation  and  available  MeHg  in  the  water  column 

J  A  -Xmaxh 

MeHg  -*■  Hgo  =  1.33  0  ~ 6  (l  - 0.56 CL)(kd31fd_MeHg  +  kdoc31fdoc_MeHg)MeHg  (4.12) 

1  , .  A  •  /I 

opht  max 


where 

MeHg—>Hgo  =  MeHg  photoreduction  yield  into  Hgo  in  water  (ng  LA  d1) 
kd3i  =  dissolved  MeHg  photoreduction  rate  in  water  (d1) 
kdoc3i  =  DOC  adsorbed  MeHg  photoreduction  rate  in  water  (d1)- 

The  surface  photoreduction  rates  ( kd3i ,  kdoc3i )  are  not  temperature 
dependent.  The  photoreduction  yield  is  subtracted  from  the  MeHg  as  a 
sink  (-)  and  added  to  the  Hgo  as  a  source  (+)  with  a  specified  yield 
coefficient  (F3i)  of  0.93. 

4.3.5  Mercury  demethylation 

Methylation  and  demethylation  are  concurrent  reactions.  Demethylation 
is  a  reverse  process  of  Hgll  methylation.  As  not  all  of  mercury  complexes 
are  available  for  methylation  and  demethylation,  speciation  of  Hgll  and 
MeHg  in  water  and  pore  water  is  an  important  factor  determining 
bioavailability  of  mercury  species.  Demethylation  of  MeHg  is  also 
influenced  by  biotic  and  abiotic  factors  (photodemethylation  or  non 
photo-mediated  demethylation  (Hobman  et  al.  2000).  Demethylation  by 
photolysis  is  assumed  to  occur  in  the  water  column.  In  HgSM,  the  water 
column  demethylation  from  MeHg  to  Hgll  is  modeled  as  a  hydrolysis 
process 


MeHg 


Hgll  =  1.33 


A  1  —  e“ 


I  M  X 

opht  n 


■h 


(l  0.56Cl  )(A:£J32  fd_MeHg  +kfoc32  fd 


doc 32  J  doc-MeHg 


)Maf/gr(4.13) 


where 
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MeHg—>HgII  =  MeHg  demethylation  yield  into  Hgll  in  water  (ng  L1  d1) 
kd.32  =  dissolved  MeHg  demethylation  rate  in  water  (d1) 
kdoc32  =  DOC  adsorbed  MeHg  demethylation  rate  in  water  (d1). 

There  are  several  possible  pathways  for  MeHg  demethylation  that  are 
discussed  in  the  literature  (Marvin-DiPasquale  et  al.  2000).  The  biotic 
process  was  suggested  to  be  the  dominant  pathway  of  MeHg 
demethylation  in  the  sediment  layer.  The  biological  demethylation  of 
MeHg  is  a  slow  process  in  contrast  to  methylation  and  is  most  effective  in 
aerobic  conditions  (Gilmour  and  Henry  1991).  The  turnover  time  of  MeHg 
demethylation  in  sediments  is  on  the  order  of  days  to  weeks.  In  HgSM,  the 
sediment  demethylation  from  MeHg  to  Hgll  is  modeled  as  a  first-order 
process  acting  on  the  dissolved  MeHg  based  on  based  on  the  algorithm 
used  in  WASP 


{MeHg  -  Hgll\  =  kd32_2(T)fd_MeHg2MeHg2  (4.14) 


where 

(MeHg—>HgII)2  =  MeHg  demethylation  yield  in  sediment  (ng  L1  d1) 

kd32-2(T ')  =  dissolved  MeHg  demethylation  rate  in  sediment  (d1). 

The  water  column  demethylation  rates  (kd32,  kdoc.32)  are  not  temperature 
dependent.  The  sediment  demethylation  rate  is  temperature  dependent 
and  corrected  using  the  AE.  An  activation  energy  value  (£),)  of  10  Real 
mol1  (41.84  KJ  mol1)  is  used  as  the  default  value.  The  demethylation 
yield,  MeHg—>HgII,  is  subtracted  from  the  MeHg  as  a  sink  (-)  and  added 
to  the  Hgll  as  a  source  (+)  with  a  specified  yield  coefficient  (  Y32)  of  0.93. 

4.4  Mercury  volatilization 

Volatilization  of  both  Hgo  and  MeHg  are  modeled  in  HgSM  at  the  surface 
of  the  water  column.  Volatilization  of  Hgll  is  ignored  in  the  model  since 
Hgll  has  much  lower  Henry’s  Law  Constants  (USEPA 1997).  The 
volatilization  rate  is  calculated  using  a  two-phase  film  resistance  theory  as 
discussed  in  CSM.  In  HgSM,  volatilization  rates  can  also  be  defined  by  the 
user,  which  is  referenced  to  20°  C. 
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4.4.1  Elemental  mercury 

Hgo  is  a  liquid  in  water  at  ambient  temperatures  and  can  emit  to  the  air 
via  volatilization.  The  volatilization  of  Hgo  represents  the  most  important 
part  of  the  exchange  between  water  and  the  air.  In  HgSM,  the  Hgo 
volatilization  across  the  air-water  interface  is  calculated  as 

Atm  — >  HgO  =  vv_Hg0  (T)j- 


HgO  — 


HgO0 


K„  (R-T 


(4.15) 


where 

Atm—>Hgo  =  Hgo  volatilization  rate  (ng  Lr1  d1) 

Vv-Hgo(T)  =  volatilization  velocity  of  Hgo  (m  d1) 

HgOo  =  air  concentration  of  Hgo  (gaseous)  (ng  L1) 

KH-Hgo  =  Henry’s  Law  constant  of  Hgo  (Pa  m3  mob1). 

Air-surface  exchange  of  Hgo  can  occur  bi-directionally,  allowing  transfer 
of  mercury  from  the  atmosphere.  Ambient  air  concentrations  of  Hgo  are 
reported  to  range  from  about  2  to  io  ngm-3,  with  the  higher  end  of  this 
range  reflecting  contributions  from  specific  local  sources  (ATSDR  2005). 

4.4.2  Methylmercury 

Similar  to  Hgo,  the  MeHg  volatilization  across  the  air-water  interface  is 
calculated  as 


Atm  MeHg  =  v 


v-MeHg 


d-MeHg 


MeHg 


MeHg0 

KH/{R-Twk) 


(4.16) 


where 


Atm—>MeHg=  MeHg  volatilization  rate  (ng  L1  d1) 
vv-MeHg(T)=  volatilization  velocity  of  MeHg  (m  d1) 

MeHgo  =  air  concentration  of  MeHg  (gaseous)  (ng  L1) 

KH-MeHg  =  Henry’s  Law  constant  of  MeHg  (Pa  m3  mol1)- 

4.5  Mercury  air  deposition,  settling  and  sedimentation 

Mercury  may  be  added  to  a  water  body  through  both  wet  and  dry 
atmospheric  deposition.  Atmospheric  deposition  of  Hgll  and  MeHg  is 
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included  in  HgSM.  In  addition,  mercury  attached  to  suspended  solids  can 
settle  onto  the  sediments  where  it  can  diffuse  into  the  water  column,  be  re¬ 
suspended,  or  be  buried  to  deep  sediments.  Mercury  settling,  sediment  re¬ 
suspension  and  the  burial  processes  are  assumed  to  operate  on  all 
particulate  fractions  of  Hgll  and  MeHg  in  HgSM. 

4.5.1  Air  deposition 

Atmospheric  deposition  of  Hgll  and  MeHg  represent  the  sum  of  wet  and 
dry  mercury  deposition.  In  HgSM,  deposition  rates  of  Hgll  and  MeHg  are 
calculated  by 


Atm  -y  Hgll  =  AsLHgII  (4.17a) 

Atm  — >  MeHg  =  AsLMeHg  (4.17b) 


where 

Atm—>HgII  =  Hgll  air  deposition  (ng  Lr1  d1) 

Atm—>MeHg  =  MeHg  air  deposition  (ng  L1  d_1) 

As  =  surface  water  area  (nr2) 

Ltigii  =  Hgll  areal  air  deposition  rate  (pg  nr2  d1) 
LMeHg  =  MeHg  areal  air  deposition  rate  (pg  nr2  d1)- 


4.5.2  Settling 

Mercury  may  enter  the  aquatic  system  in  the  insoluble  form  adsorbed  onto 
suspended  solids.  Suspended  solids  can  then  settle  either  to  the  sediment 
layer  or  to  an  underlying  water  compartment.  Settling  of  suspended  solids 
will  remove  adsorbed  mercury  species  in  proportion  to  their  bulk 
concentration.  In  HgSM,  settling  rates  of  solids  adsorbed  Hgll  and  MeHg 
from  the  water  column  are  calculated  by 


Hgll  Bed  =  -- 


\  v  f  -\-v  f  -) -v  f  Hpii  (4.18a) 

/  /  usrmJ  nn  Hull  1  xnnJ  nn-HnU  1  ^ snmJ  nnm-HnTT  l.) 1 1  '  ' 


MeHg  -►  Bed  =  -— 


^  ^  VSpn  jpn  MeHg  ^  sap  f  ap- MeHg  ^somf  pom— MeHg 


MeHg  (4.18b) 


where 


HgII—>Bed  =  Hgll  settling  rate  in  water  (ng  Lr1  d1) 
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MeHg—>Bed  =  MeHg  settling  rate  in  water  (ng  L 1  d1)- 


4.5.3  Re-suspension 

The  re-suspension/erosion  of  sediments  refers  to  the  process  by  which 
solids,  and  their  adsorbed  contaminants,  are  transported  from  the 
sediment  layer  into  the  water  column.  In  HgSM,  the  re-suspension  rates  of 
Hgll  and  MeHg  are  calculated  respectively  as  the  product  of  sediment 
adsorbed  concentration  and  the  re-suspension  velocity 


Bed  Hgll 


_1 

h. 


N 

^  rnn  J  pn—HgII2 
71  —  1 


f 

nJ  pi 


Hgll  2 


(4.19a) 


Bed  — >  MeHg 


1 


N 

y  '  ^ rnri .f  i pn—MeHg2 
71  =  1 


MeHg2 


(4.19b) 


where 


Bed—  > Hgll  =  sediment  Hgll  re-suspension  rate  (ng  L1  d1) 
Bed—>MeHg  =  sediment  MeHg  re-suspension  rate  (ng  Lr1  d1). 


4.5.4  Sediment  burial 

As  new  sediments  are  deposited,  an  equal  mass  of  sediment  is  lost  to 
burial  due  to  the  constant  sediment  layer  assumption.  In  HgSM,  the  fluxes 
of  Hgll  and  MeHg  out  of  the  active  sediment  layer  and  into  underlying 
sediments  are  calculated  respectively  as  the  product  of  active  sediment 
adsorbed  mercury  concentration  and  burial  velocity 


HgII2  burial  —  — u, 
h2 


N 

E/p.  —Hgll 2  dporn  Hgll  2 

V  71=1 


HgII2 


(4.20a) 


MeHg2  burial  =  — v, 
h2 


(  N 

E/, 


/ 
J  Dl 


pn—MeHg2  '  J  pom— MeHg 2 


MeHg2  (4.20b) 


where 

Hgll 2  burial  =  sediment  Hgll  burial  rate  (ng  L1  d1) 
MeHg 2  burial  =  sediment  MeHg  burial  rate  (ng  L_1  d1)- 
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4.5.5  Sediment-water  transfer 

Dissolved  mercury  species  in  the  sediment  layer  can  transfer  across  the 
sediment-water  interface,  or  vice  versa.  Pore  water  mass  transfer  is 
assumed  to  operate  on  all  dissolved  fractions  of  Hgll  and  MeHg.  In 
HgSM,  sediment-water  transfer  fluxes  of  Hgll  and  MeHg  are  calculated  by 

Hgll  ~  HgII2  =  ^\(fd_HgIT2  +  fdoc_  HgII2 ) HgII2 /cp  - (fd_HgII  +  fdoc_HgU ) Hgll]  (4.21a) 
MeHg  ~  MeHg2  =  ^-[(fd  MeHg2  +  fdoc  MeHg2)MeHg2/<p-(fd  MeHg+fd0C  MeHg)MeHg\  (4.21b) 


where 

HgII<—>HgII2  =  Hgll  sediment-water  transfer  (ng  L1  d1) 

MeHg<—>MeHg2  =  MeHg  sediment-water  transfer  (ng  L1  d1)- 

A  positive  flux  results  in  the  transport  of  mercury  from  the  sediment  to  the 
overlying  water,  a  negative  flux  means  the  transport  of  mercury  from  the 
overlying  water  to  the  sediment.  A  high  concentration  of  MeHg  at  the 
interface  between  the  oxic  and  anoxic  layers  could  be  a  result  of  either 
production  or  accumulation  of  settling  particulate  matter  (Cossa  et  al. 
1994).  The  sediment -water  mass-transfer  velocity  of  Hgll  and  MeHg  is  a 
user  specified  parameter  or  internally  computed  in  the  model.  Four 
alternative  sediment-water  transfer  equations  implemented  in  HgSM  are 
based  on  thoses  included  in  CSM:  1)  Thibodeaux  et  al.  (2001),  2)  Di  Toro 
et  al.  (1981),  3)  Boyer  et  al.  (1994),  and  4)  Schink  and  Guinasso  (1977). 

4.6  Water  column  mercury  source/sink  equations 

The  pathways  for  three  mercury  species  are  depicted  in  Figure  6.  Major 
processes  modeled  in  HgSM  include  partitioning  of  Hgll  and  MeHg, 
particulate  settling,  re-suspension  and  burial,  sediment-water  diffusion, 
volatilization,  and  transformations.  Meanwhile,  mercury  species  can  be 
transported  together  with  water,  DOC,  and  suspended  solids.  All 
concentrations  in  source  and  sink  term  equations  of  mercury  state  variables 
are  expressed  in  terms  of  mass  per  unit  volume  of  water  plus  solids  (ng  L1). 
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4.6.1  Elemental  mercury 

Hgo  is  the  predominant  mercury  cycling  product  in  the  water  column. 
Major  processes  involved  in  Hgo  cycling  in  the  water  column  include  Hgll 
reduction  to  Hgo,  demethylation,  and  volatilization.  The  internal  source 
(+)  and  sink  (-)  terms  of  Hgo  computed  in  HgSM  are  listed  in  Table  18. 


Table  18.  HgO  source/sink  terms  and  pathways. 


Source/Sink  term* 

Pathway 

11 

0 

£ 

co 

Co 

V  1  f  HrO  Hg0°  1 

kh/Rtk) 

HgO  volatilization 

—  k\2{T)  ■  HgO 

HgO  oxidation 

+ 1  -33  7°  1  ~  6  — "  (1  -  0.56Q  tkd2lfd_HgI1  +  kdoc2Jdoc_HglI  )Hgll  •  Y21 

*  Opht  Anax  '  ^ 

Hgll  photoreduction 

+ 1  •33-^1~g  "’(l  -  0.56C,  tkdnfd_MeHg  +  kdoc3Jdoc_MeHg )MeHg  •  Y31 

1  Opht  Anax  *  ^ 

MeHg  photoreduction 

*  SSwgo  is  the  sum  of  internal  source  and  sink  terms  of  HgO  in  water  (ng  L 1  d-1). 
Y 2i  is  the  Hgll  photoreduction  yield  coefficient  in  water. 

Y31  is  the  MeHg  photoreduction  yield  coefficient  in  water. 


4.6.2  Inorganic  mercury 

Major  processes  involved  in  Hgll  cycling  in  the  water  column  include  Hgll 
reduction,  methylation,  settling  of  adsorbed  particulates  into  the 
underlying  sediment  layer,  Hgll  diffusion  across  the  sediment-water 
interface,  and  sediment  re-suspension.  Hgll  may  be  adsorbed  by  DOC, 
algae,  POM,  and  suspended  solids.  Both  equilibrium  and  non-equilibrium 
partitioning  of  Hgll  are  allowed  in  HgSM.  Under  the  non-equilibrium 
partitioning  option,  the  internal  source  and  sink  terms  are  computed 
separately  for  potential  partitioning  phases  in  a  similar  manner  as  is  done 
for  contaminants. 

4.6.2. 1  Equilibrium  partitioning 

Under  equilibrium  partitioning  implementation,  the  internal  source  (+) 
and  sink  (-)  terms  of  water  column  Hgll  are  computed  in  HgSM  for  the 
total  concentration  and  are  listed  in  Table  19. 
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Table  19.  Water  column  total  Hgll  source/sink  terms  and  pathways  under  equilibrium  partitioning. 


Source/Sink  term* 

Pathway 

ssHgII  = 

+  A^Hgll 

Atmospheric  Hgll 
deposition 

1 

h 

f  N  \ 

^ spnf pn-Hgll  sapf ap-Hgll  ^  ^ somf  pom-Hgll  ^igll 

V  «=1  7 

Hgll  settling 

1 

+  — 
h 

(  N  A 

^ \V rpnf pn-HgII2  HgII2 

V  n=l  J 

Hgll  re-suspension 

+  y  lfd-Hg,n  +  fdoc-Hg,n  )HgII2  y  -  (fd  HgU  +  fdocHgll  )HgIl} 

Hgll  sediment-water 
transfer 

+  kn(T)  ■  HgO  ■  Yl2 

HgO  oxidation 

~  1  -33 - 0-56Q ){kd2  Jd  _Hg,,  +  kdoc2  Jdoc_HgU )HgII 

*  Opht  Anax  '  ^ 

Hgll  photoreduction 

~  (^23  (r)fd_Hg„  +  kdoc  23  (r)fd_HgII  )HgII 

Hgll  methylation 

+ 1.33-^  -  0.56 CL  tki21fd_MeHg  +  kdocl2fdoc_UMg)MeHg  ■  732 

-*0pht  Anax  ' 

MeHg  demethylation 

*  SShigii  is  the  sum  of  internal  source  and  sink  terms  of  total  Hgll  in  water  (ng  L 1  d1). 
Y12  is  the  HgO  oxidation  yield  coefficient  in  water. 

Y32  is  the  MeHg  demethylation  yield  coefficient  in  water. 


4. 6.2.2  Non-equilibrium  partitioning 

Under  non-equilibrium  partitioning  implementation,  the  state  variables 
for  Hgll  include  concentrations  of  potential  phases.  The  internal  source 
and  sink  equations  of  Hgll  are  developed  in  a  similar  manner  as  for  non¬ 
equilibrium  partitioning  of  contaminants  implemented  in  CSM.  The 
internal  source  (+)  and  sink  (-)  terms  of  water  column  Hgll  are  computed 
in  HgSM  for  dissolved,  DOC  adsorbed,  algae  adsorbed,  POM  adsorbed, 
and  solids  adsorbed  Hgll  and  are  listed  in  Table  20. 


Table  20.  Water  column  Hgll  partitioning  phase  source/sink  terms  and  pathways  under  non¬ 
equilibrium  partitioning. 


Source/Sink  term* 

Pathway 

Dissolved  Hgll 

ssHgl/d  = 

+  A^Hgll 

Atmospheric  Hgll 
deposition 
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Source/Sink  term* 

Pathway 

h 

Hgll  sediment-water 
transfer 

+  kn(T)  ■  HgO  ■  Yl2 

HgO  oxidation 

- 1.33-^- 1~g  "'(1  -  0.56Q  )kd2iHgIId 

*  Opht  Anax  *  " 

Hgll  photoreduction 

-kd23(T)HgIId 

Hgll  methylation 

+  1.33-^1~e  "  (1  - 0.56Q )[fe32 (r )MeHg d )  +  {kdoc32 (T )MeHgdoc )]Y32 

-*0pht  Anax  '  " 

MeHg  demethylation 

-(10'3^(  T)HgIId(qcapApd-HgIlJ-kdaap{T)HgIlJ 

-  (l O'3 kadpom  (T)HgIId  (qcpomPOM  -  Hgll pom )-  kdapom  (T)HgIIpom ) 

-  Z  I1 0_3  KJT)  HgIId  (qcnmn  -  Hgll pn )  -  kdan  ( T)  Hgll pn  ] 

n= 1 

Adsorption-desorption 

DOC  adsorbed  Hgll 

- 

Hglldoc 

+  '-f[HgIIdoc2/<t>-HgIIdoc ] 
h 

Hgll  sediment-water 
transfer 

-1.33-/^1~g  '""'(I - 0.56Q )kdoc2iHgUdoc 

*  Opht  ^max  '  ™ 

Hgll  photoreduction 

—  ^ doc  23  {T)Hgkkdoc 

Hgll  methylation 

Algae  adsorbed  Hgll 

ssHgIIap  = 

Hgll  settling 

+^Ktv(T)HgII^ApJ-HgIlJ-kd^(T)HgII,r 

Adsorption-desorption 

POM  adsorbed  Hgll 

99  - 

Hgllpom  ~ 

pom 

Hgll  settling 

+W-X^(nHgII^rmPOM-HgIIrJ-kimm(T)HgIIf„ 

Adsorption-desorption 

Solids  adsorbed  Hgll 

ssHglIpn  = 
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Source/Sink  term* 

Pathway 

-jVspMgHpn 

Hgll  settling 

+jvrptIHgnPn  2 

Hgll  re-suspension 

Adsorption-desorption 

*  SS Hgiid  is  the  sum  of  internal  source  and  sink  terms  of  dissolved  Hgll  in  water  (ng  L1  d1). 


SSngiidoc  is  the  sum  of  internal  source  and  sink  terms  of  DOC  adsorbed  Hgll  in  water  (ng  L1  d1). 
SSHguap  is  the  sum  of  internal  source  and  sink  terms  of  algae  adsorbed  Hgll  in  water  (ng  L1  d-1). 
SSHgiipom  is  the  sum  of  internal  source  and  sink  terms  of  POM  adsorbed  Hgll  in  water  (ng  L1  d-1). 
SSHgup  is  the  sum  of  internal  source  and  sink  terms  of  solids  adsorbed  Hgll  in  water  (ng  L 1  d1). 


4.6.3  Methylmercury 

Major  processes  involved  in  MeHg  cycling  in  the  water  column  include 
volatilization,  methylation,  photodegradation,  settling  of  adsorbed 
particulates  into  the  sediment  layer,  MeHg  diffusion  across  the  sediment- 
water  interface,  and  sediment  re-suspension  and  burial.  MeHg  may  be 
adsorbed  by  DOC,  algae,  POM,  and  suspended  solids.  Only  equilibrium 
partitioning  of  MeHg  is  modeled  in  HgSM.  The  internal  source  (+)  and 
sink  (-)  terms  of  water  column  MeHg  are  computed  in  HgSM  for  the  total 
concentration  and  are  listed  in  Table  21. 


Table  21.  Water  column  total  MeHg  source/sink  terms  and  pathways  under  equilibrium  partitioning. 


Source/Sink  term* 

Pathway 

ssMeHg  = 

Atmospheric  MeHg 
deposition 

1 

I 

(  N  \ 

^ spn-f pn-MeHg  sapf ap-MeHg  ^ somf pom-MeHg  MeHg 

K  n= 1  ) 

MeHg  settling 

1 

+  — 
h 

f  N  A 

YuV<pnfpn-MeHg  2  MeHg2 

\n= 1  7 

MeHg  re-suspension 

+  ^  ^fd-MeHg2+ fdoc-MeHgl^^gll^  if d-MeHg  +  f doc-MeHg)^e^g\ 

MeHg  sediment-water 
transfer 

v  1  ff  MeHg  1 

*-MeHg  h  ^  d-MeHg  g  Kh  /RTk) 

MeHg  volatilization 

-1.33-^-  '  ''  '  I1  -  0-56Q  )(^3  Jd-MeHg  +  ***31  fdoc-MeHg  K Hg 

*  Opht  Anax  *  " 

MeHg  photoreduction 
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+  (^/23  (T)fd_Hg„  +  kdoc23  (T)  f doc-Hgll  '  ^23 

Hgll  methylation 

— 1.33-^^  ""“"(l - 0.56 CLikd32fd_MeHg  +  kdoc32fdoc_MeHg)MeHg 

-'Opht  "unax  *  ^ 

MeHg  demethylation 

*  SSiweHgis  the  sum  of  internal  source  and  sink  terms  of  total  MeHg  in  water  (ng  L1  d1). 
V23  is  the  Hgll  methylation  yield  coefficient  in  water. 


4.7  Sediment  mercury  mass  balance  equations 

The  major  species  considered  for  sediment  cycling  are  Hgll  and  MeHg, 
Hgo  is  not  included  for  the  sediment  layer  (Figure  6).  Major  processes 
modeled  in  HgSM  include  partitioning  of  Hgll  and  MeHg,  sediment  re¬ 
suspension  and  burial,  sediment-water  diffusion,  and  transformations. 
Because  of  varying  environmental  conditions,  the  rate  of  the  reactions  and 
transformations  can  vary  dramatically  from  those  in  the  water  column. 
Similar  to  the  CSM,  the  sediment  layer  is  assumed  to  have  constant 
properties  including  the  thickness,  volume,  porosity,  and  bulk  density.  In 
the  sediment  layer,  concentrations  of  mercury  state  variables  are 
expressed  in  terms  of  mass  per  unit  volume  of  total  sediments  (ng  L1). 

4.7.1  Inorganic  mercury 

Major  processes  involved  in  Hgll  cycling  in  the  sediment  layer  include 
methylation,  demethylation,  water  column  settling,  pore  water  diffusion 
across  the  sediment-water  interface,  and  sediment  re-suspension  and 
burial.  Hgll  may  be  adsorbed  by  DOC,  sediment  POM  and  solids. 

4. 7.1.1  Equilibrium  partitioning 

Under  equilibrium  partitioning  implementation,  the  source  (+)  and  sink 
(-)  terms  of  sediment  Hgll  are  computed  in  HgSM  for  the  total 
concentration  are  listed  in  Table  22. 


Table  22.  Sediment  layer  total  Hgll  source/sink  terms  and  pathways  under  equilibrium  partitioning. 


Source/Sink  term* 

Pathway 

ssHgH2  = 

1 

+  — 
K 

(  N  \ 

^ spnf  pn-Hgll  sapf ap-Hgll  ^ somJ pom- Hgll 

\n= 1  y 

Hgll  deposition 

l 

h2 

^  N  \ 

Vrpnf pn-Hgll  2  HgII2 

\  n= 1  ) 

Hgll  erosion 
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vb  (x'  ^ 

j  f pn  -  Hgll  2  f pom- Hgll  2  \^S^2 

"2  \n=l  ) 

Hgll  deep  burial 

~  Y  lL-Hg,n  +  .Lc-Hg„2  Wh  I4>  ~  +  fdoc-HgU  hM 

Hgll  sediment-water  transfer 

f  S 04  ^ 

(D  „  S042rmsoJd_HgII2HgII2 

Hgll  methylation 

+  b  dy2-lf d-MeHglHeHgJY^ 

MeHg  demethylation 

*  SSngii2  is  the  sum  of  source  and  sink  terms  of  total  Hgll  in  sediment  (ng  L 1  d1). 


4.7. 1.2  Non-equilibrium  partitioning 

Non-equilibrium  kinetics  of  Hgll  onto  sediment  solids  is  modeled 
identically  to  the  partitioning  of  Hgll  onto  suspended  solids  in  the  water 
column.  However,  the  values  for  the  number  of  binding  sites  for  sediments 
can  be  different,  thus  providing  the  ability  to  reflect  the  difference  in 
adsorption  characteristics  between  sediment  solids  and  solids  suspended 
in  the  water  column.  Under  a  non-equilibrium  partitioning 
implementation,  the  mass  balances  of  Hgll  in  each  phase  must  be  solved 
simultaneously.  The  source  (+)  and  sink  (-)  terms  of  sediment  Hgll  are 
computed  in  HgSM  for  dissolved,  DOC  adsorbed,  POM  and  solids 
adsorbed  Hgll  and  are  listed  in  Table  23. 


Table  23.  Sediment  layer  Hgll  partitioning  phase  source/sink  terms  and  pathways  under  non-equilibrium 

partitioning. 


Source/Sink  term* 

Pathway 

Dissolved  Hgll 

CC  - 

H gild  2 

-MHgnd2/</>-Hgiid\ 

h2 

Hgll  sediment-water 
transfer 

(  S04 

-kso42(T)  - - S042rmsoJd_Hgn2HgII2 

V*Uo4  bt)42  ^ 

Hgll  methylation 

+  bd  32-2  f d-MeHglMeHgjY^ 

MeHg  demethylation 

- (' 0-Hadpom2(T)Hgrrd2{qcpom2POM2  -  HgJIpoml)- kdapom2(T)HgITpom2) 

-  2 1  ®~3  ^adni  (T)HgIId2  {qai2mn2  -  Hgll pn2 )  -  kdan2  H}HgIIpn2  ] 

n= 1 

Adsorption-desorption 

DOC  adsorbed  Hgll 

99 

kJL>Hglldoc2 
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-HHgHdoc2/</>-HgIIdo  J 
h2 

Hgll  sediment-water 
transfer 

POM  adsorbed  Hgll 

SS Hgllpoml  — 

+  Y  (VsapHgIIap  +  Vso,nHgnpom  ) 

Hgll  deposition 

Hgll  deep  burial 

+ 1  ^kadpom2 ( T)HgIId2 (qcpomlPOM 2  -  HgIIpom 2 )-  kdapoml(T)HgIIpoml 

Adsorption-desorption 

Solids  adsorbed  Hgll 

SSHgnpn2  — 

SSHgIIpn2  =  Y  Vspn  UgH  pn 

Hgll  deposition 

-YVrpnHgnpn2 

Hgll  re-suspension 

-v!;vgUpnl 

h2 

Hgll  deep  burial 

Adsorption-desorption 

*  SS Hgiid2  is  the  sum  of  source  and  sink  terms  of  dissolved  Hgll  in  sediment  (ng  L 1  d1). 


SSHgiidoc2  is  the  sum  of  source  and  sink  terms  of  DOC  adsorbed  Hgll  in  sediment  (ng  L1  d-1). 
SSHgiipom2  is  the  sum  of  source  and  sink  terms  of  POM  adsorbed  Hgll  in  sediment  (ng  L 1  d-1). 
SSngiip2  is  the  sum  of  source  and  sink  terms  of  solids  adsorbed  Hgll  in  sediment  (ng  L1  d-1). 


4.7.2  Methylmercury 

Major  processes  involved  in  MeHg  cycling  in  the  sediment  layer  include 
methylation,  biological  demethylation,  pore  water  diffusion  across  the 
sediment-water  interface,  water  column  settling,  sediment  re-suspension, 
and  burial.  MeHg  may  be  adsorbed  by  DOC,  sediment  POM  and  solids. 
The  source  (+)  and  sink  (-)  terms  of  sediment  MeHg  are  computed  in 
HgSM  for  the  total  concentration  and  are  listed  in  Table  24. 
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Table  24.  Sediment  layer  total  MeHg  source/sink  terms  and  pathways  under  equilibrium  partitioning. 


Source/Sink  term* 

Pathway 

OC 

MeHgl 

1 

H - 

h2 

f  N  \ 

^  spnf pn-MeHg  sapfap-MeHg  ^ somf  pom-MeHg  MeHg 

\  n= 1  J 

MeHg  deposition 

1 

h2 

(  N  A 

YjVrpnfpn-MeHg2  MeHg2 
\n= 1  y 

MeHg  erosion 

^  \fd-MeHg2  +  f doc-MeHgl  ) lMeHg2  j (j)  {fd_MeHg  +  fdoc-MeHg  ^MeHg  ] 

MeHg  sediment-water 
transfer 

1 

h2 

(  N  \ 

f pn -MeHg 2  J pom- MeHgl  -MeHg 2 

\n= 1  J 

MeHg  deep  burial 

(  SOA  ^ 

+  *„« (H  - SOA2rmsoJd_HgII2HgII2Y23 

V  ^ sod  ^ 

Hgll  methylation 

_  ^dil-lfd-MeHgl^^Sl 

MeHg  demethylation 

*  SSivieHg2  is  the  sum  of  source  and  sink  terms  of  total  MeHg  in  sediment  (ng  L 1  d1). 


4.8  HgSM  parameters 

The  HgSM  operates  based  on  the  schematization  already  set  up  for  the 
H&H  model  and  the  flows  computed  by  the  H&H  model.  Table  25  provides 
a  summary  of  mercury  input  parameters  that  must  be  specified  when  all 
options  are  turned  on  in  the  HgSM  module.  Literature  values  were  taken 
from  a  compilation  of  several  studies  if  they  are  available.  Most  of  them  are 
model  calibration  parameters  in  real  world  applications.  Similar  to  CSM, 
HgSM’s  parameters  for  the  water  column  and  sediment  layer  can  be 
specified  as  either  a  uniform  constant,  a  single  for  all  water  quality  cells,  or 
spatially-varying  values  by  user  defined  water  quality  regions.  There  are 
three  groups  of  parameters:  global,  water  column  and  sediment  layer.  This 
table  will  be  repeated  for  each  water  quality  region,  allow  the  user  to  define 
the  different  values  for  input  parameters. 


Table  25.  List  of  HgSM’s  input  parameters. 


Symbol 

Definition 

Default 

values 

Approximate 

range 

Units 

Temp 

correction 

Global 

Dm 

Molecular  diffusivity 

- 

n/a 

m2  d-1 

Vm 

Sediment-water  mass  transfer 
velocity 

- 

n/a 

m  d'1 
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Symbol 

Definition 

Default 

values 

Approximate 

range 

Units 

Temp 

correction 

hi 

Sediment  layer  thickness 

0.1  - 
0.15 

m 

Zl 

Sediment  bioturbation  layer 
thickness 

0.05  - 
0.1 

m 

Vss 

Solids  settling  velocity 

- 

n/a 

m  d-1 

Vsom 

Organic  matter  settling  velocity 

- 

n/a 

m  d'1 

kpht(T) 

Aquatic  photolysis  rate 

- 

n/a 

d-1 

Ea 

n/a 

lopht 

Light  intensity  when  kpht  is 
measured 

- 

n/a 

Wm-2 

CXl 

Light  attenuation  adjusting 
coefficient 

1.33 

1.2  - 1.6 

unitless 

Water  column  -  elemental  mercury 

MW 

HgO  molecular  weight3 

200.6 

- 

g  mol1 

S  HgO 

HgO  solubility 

56 

- 

pg  L1 

Vv-Hgo(T) 

HgO  volatilization  velocityf 

0.006 

0.0059  -  0.45 

m  hr1 

0 

n/a 

Kh 

HgO  Henry’s  Law  constant 

0.09' 

- 

Pa  m3  mol1 

HgOo 

HgO  air  concentration 

2.10-3 

n/a 

ng  L1 

kn(T) 

HgO  oxidation  rateb 

-3 

10 

-3  -1 

10  -10 

d-1 

Ea 

n/a 

Y  12 

HgO  oxidation  yield  coefficient 

1.0 

0-2.0 

unitless 

Water  column  -  inorganic  mercury 

MW 

Hgll  molecular  weight3 

271.52 

232.68  (HgS) 
271.52  (HgCh) 

g  mol1 

Kp-Hgll 

Hgll  equilibrium  partition  coefficient 
for  algaec 

- 

5  6 

10  -10 

L  kg 1 

Kp-Hgll 

Hgll  equilibrium  partition  coefficient 
for  si ltb 

2.105 

3  6 

10  -10 

L  kg 1 

Kp-Hgll 

Hgll  equilibrium  partition  coefficient 
for  clayb 

2.105 

3  6 

10  -10 

L  kg 1 

Kp-Hgll 

Hgll  equilibrium  partition  coefficient 
for  solidse 

5.3 

10 

4.2  6.9 

10  -10 

L  kg 1 

bHgll 

Hgll  Freundlich  exponent 

- 

0.4  -  1.2 

unitless 

Kf-Hgll 

Hgll  Freundlich  adsorption  constant*1 

4.5-104- 

2.52-108 

(Mg  g  !)  (Mg  L1)^ 

Kl-Hgll 

Hgll  Langmuir  adsorption  constant11 

- 

6  6 

10  -  7-10 

L  Pg 1 

Cjcn-Hgll 

Hgll  adsorption  capacity  for  solid 

“n”h 

2.8-102  - 
3.58-103 

hgg1 

kadn-Hgll 

Hgll  adsorption  coefficient 

- 

n/a 

L  pg1  d-1 
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Symbol 

Definition 

Default 

values 

Approximate 

range 

Units 

Temp 

correction 

kdan-Hgll 

Hgll  desorption  rate 

1.0 

n/a 

d-1 

kdoc-Hgll 

Hgll  partition  coefficient  for  DOCe 

5.3 

10 

105-3-  105-6 

L  kg 1 

kd21 

Dissolved  Hgll  photoreduction  rateb 

5.10  2 

-3  -1 

10  -5-10 

d-1 

kdoc21 

DOC  adsorbed  Hgll  photoreduction 
rateb 

0.0 

-3  -1 

10  -5-10 

d-1 

Y21 

Hgll  photoreduction  yield  coefficient 

1.0 

0-2.0 

unitless 

kd23(T) 

Dissolved  Hgll  methylation  rateb 

-3 

10 

-5  -2 

10  -5-10 

d-1 

0 

1.14 

kdoc23(T) 

DOC  adsorbed  Hgll  methylation 
rateb 

-3 

10 

-5  -2 

10  -5-10 

d-1 

e 

1.14 

Y 23 

Hgll  methylation  yield  coefficient 

1.07 

0-2.0 

unitless 

Water  column  -  methymercury 

MW 

MeHg  molecular  weight3 

230.66 

(CH32Hg) 

n/a 

g  mol1 

Kp-MeHg 

MeHg  equilibrium  partition 
coefficient  for  algaee 

5 

10 

5  7 

10  -10 

L  kg 1 

Kp-MeHg 

MeHg  equilibrium  partition 
coefficient  for  siltb 

2- 10  5 

3  6 

10  -10 

L  kg 1 

Kp-MeHg 

MeHg  equilibrium  partition 
coefficient  for  calyb 

2- 10  5 

3  6 

10  -10 

L  kg 1 

Kp-MeHg 

MeHg  equilibrium  partition 
coefficient  for  solidse 

5.4 

10 

4.2  6.2 

10  -10 

L  kg 1 

Kdoc-MeHg 

MeHg  equilibrium  partition 
coefficients  for  DOCb 

2-105 

5  6 

10  -10 

L  kg 1 

bMeHg 

MeHg  Freundlich  exponent 

- 

n/a 

unitless 

Kf-MeHg 

MeHg  Freundlich  adsorption 
constant 

- 

n/a 

(Mg  g1)  (Mg  L-1)-13 

Kl-MeHg 

MeHg  Langmuir  adsorption 
constant11 

- 

n/a 

L  Mg 1 

Qcn-MeHg 

MeHg  adsorption  capacity  for  solid 

“n” 

- 

n/a 

Mgg1 

kd31 

Dissolved  MeHg  photoreduction  rate 
into  HgOb 

- 

-3  -1 

10  -5-10 

d-1 

kdoc31 

DOC  adsorbed  MeHg 
photoreduction  rate  into  HgOb 

0.0 

-3  -1 

10  -5-10 

d-1 

Y31 

MeHg  photoreduction  yield 
coefficient 

0.93 

0-2.0 

unitless 

kd32 

Dissolved  MeHg  demethylation  rate 
into  Hgllb 

5-10  2 

-3  -1 

10  -5-10 

d-1 

kdoc32 

DOC  adsorbed  MeHg  demethylation 
rate  into  Hgllb 

0.0 

-3  -1 

10  -5-10 

d-1 
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Symbol 

Definition 

Default 

values 

Approximate 

range 

Units 

Temp 

correction 

Y32 

MeHg  demethylation  yield 
coefficient 

0.93 

0-2.0 

unitless 

Vv-MeHg(T) 

MeHg  volatilization  velocity 

1.9-105 

n/a 

m  d-1 

0 

n/a 

Kh 

MeHg  Henry’s  constant 

4.5-10  6| 

n/a 

Pa  m3  mol1 

MeHgo 

MeHg  air  concentration 

0.0 

n/a 

ngL1 

Sediment  layer  -  inorganic  mercury 

Kp-Hgll  2 

Hgll  equilibrium  partition  coefficient 
for  si ltb 

- 

3  6 

10  -10 

L  kg 1 

Kp-Hgll2 

Hgll  equilibrium  partition  coefficient 
for  clayb 

- 

3  6 

10  -10 

L  kg 1 

Kp-Hgll2 

Hgll  equilibrium  partition  coefficient 
for  solidse 

4.9 

10 

3.8  6 

10  -10 

L  kg 1 

bHgll2 

Hgll  Freundlich  exponent 

- 

n/a 

unitless 

Kf-Hgll2 

Hgll  Freundlich  adsorption  constant 

- 

n/a 

(Mg  g1)  (Mg  L'1)'13 

KlHgll2 

Hgll  Langmuir  adsorption  constant 

- 

51  -  390 

L  Mg 1 

Cfcn-Hgll2 

Hgll  adsorption  capacity  for  solid  "n” 

- 

n/a 

Mg  g'1 

kadn-Hgll2 

Hgll  adsorption  coefficient 

- 

n/a 

L  pg 1  d  1 

kdan-Hgll2 

Hgll  desorption  rate 

0.1 

n/a 

d'1 

kdoc-Hgll2 

Hgll  partition  coefficient  for  DOCb 

- 

4  5 

10  -10 

L  kg 1 

kso42(T) 

Sediment  S04  reduction  rate 

- 

- 

d-1 

e 

n/a 

Ks04 

Half-saturation  constant  for  the 
effect  of  S04  on  methylation 

- 

- 

mg-02  L1 

r mso4 

Ratio  of  sediment  methylation  rate 
and  sulfate  reduction  rate 

- 

- 

L  mg1 

Sediment  layer  -  methymercury 

Kp-MeHg2 

MeHg  equilibrium  partition 
coefficient  for  siltb 

- 

3  6 

10  -10 

L  kg 1 

Kp-MeHg2 

MeHg  equilibrium  partition 
coefficient  for  clayb 

- 

3  6 

10  -10 

L  kg 1 

Kp-MeHg2 

MeHg  equilibrium  partition 
coefficient  for  solidse 

3.6 

10 

2.8  5 

10  -10 

L  kg 1 

Kdoc-MeHg2 

MeHg  equilibrium  partition 
coefficients  for  DOCb 

- 

5  6 

10  -10 

L  kg 1 

b)MeHg2 

MeHg  Freundlich  exponent 

- 

n/a 

Unitless 

Kf~MeHg2 

MeHg  Freundlich  adsorption 
constant 

- 

n/a 

(Mgg-1)(pgL-1)-b 

Kl~MeHg2 

MeHg  Langmuir  adsorption 
constant 

- 

n/a 

L  Mg 1 
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Symbol 

Definition 

Default 

values 

Approximate 

range 

Units 

Temp 

correction 

Qcn-MeHg2 

MeHg  adsorption  capacity  for  solid 

“n” 

- 

n/a 

Mgg1 

kd32-2(T) 

Sediment  MeHg  demethylation  rate 

0.21 

- 

d-1 

Ea  n/a 

a.  ATS  DR  (2005). 

b.  Wool  (etal.  2006). 

c.  Hudson  et  al.  (1994). 

d.  Tsiros  and  Ambrose  (1999). 

e.  Allison  and  Allison  (2005). 

f.  Loux  (2004). 

g.  Miles  et  al.  (2001). 

h.  Chen  et  al.  (2009). 

I.  Lin  etal.  (2012). 

j.  Gilmour  et  al.  (2007). 


As  discussed  in  the  CSM  module,  the  sediment  layer  has  a  constant 
volume  and  thickness  during  a  simulation.  Sediment  particle  density  and 
porosity  are  fixed.  Settling  velocities  of  particles  and  sediment  re¬ 
suspension  rate  parameters  can  be  computed  from  the  GC  module. 
Benthic  sediment  layer  parameters  have  been  discussed  in  the  GC  module 
and  are  listed  in  Table  2. 

4.9  HgSM  outputs 

This  section  summarizes  the  HgSM  outputs  from  a  model  simulation.  The 
fundamental  output  produced  by  HgSM  consists  of  predicted 
concentrations  and  predicted  mass  fluxes  of  mercury  species  at  specified 
locations  within  the  system.  Concentrations  of  the  Hgll  and  MeHg 
computed  for  each  physical  compartment  (water  column  and  sediment 
layer)  include  dissolved,  DOC  adsorbed,  algae,  POM,  and  solids  adsorbed 
phases.  Table  26  lists  the  symbols  and  definitions  of  concentrations  of 
mercury  species  computed  in  HgSM. 
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Table  26.  List  of  concentrations  of  mercury  species  computed  in  HgSM. 


Symbol 

Definition 

Units 

Elemental  mercury 

HgO 

HgO  concentration  in  water 

ng  L1 

Inorganic  mercury 

Hgll 

Concentration  of  total  Hgll  in  water 

ng  L-1 

Hglla 

Concentration  of  dissolved  Hgll  in  water 

ng  L-1 

Hgll  doc 

Concentration  of  DOC  adsorbed  Hgll  in  water 

ng  L-1 

Hgllap 

Concentration  of  algae  adsorbed  Hgll  in  water 

ng  L-1 

Hgllpom 

Concentration  of  POM  adsorbed  Hgll  in  water 

ng  L-1 

Hgllpt 

Total  concentration  of  solids  adsorbed  Hgll  in  water 

ng  L-1 

Hgllpts 

Total  concentration  of  solids  adsorbed  Hgll  in  water 

ngg1 

Hglh 

Concentration  of  total  Hgll  in  sediment 

ng  L-1 

Hglldp2 

Concentration  of  dissolved  Hgll  in  pore  water 

ng  L-1 

Hglldocp2 

Concentration  of  DOC  adsorbed  Hgll  in  pore  water 

ng  L-1 

Hgllpom2 

Concentration  of  POM  adsorbed  Hgll  in  sediment 

ng  L-1 

Hgllpt2 

Total  concentration  of  solids  adsorbed  Hgll  in  sediment 

ng  L1 

Hgllpts2 

Total  concentration  of  solids  adsorbed  Hgll  in  sediment 

ngg1 

Methymercury 

MeHg 

Concentration  of  total  MeHg  in  water 

ng  L-1 

MeHgd 

Concentration  of  dissolved  MeHg  in  water 

ng  L-1 

McHgdoc 

Concentration  of  DOC  adsorbed  MeHg  in  water 

ng  L-1 

MeHgap 

Concentration  of  algae  adsorbed  MeHg  in  water 

ng  L-1 

MoHgpom 

Concentration  of  POM  adsorbed  MeHg  in  water 

ng  L1 

MeHgpt 

Total  concentration  of  MeHg  adsorbed  on  solids  in  water 

ng  L-1 

MeHgpts 

Total  concentration  of  MeHg  adsorbed  on  solids  in  water 

ngg1 

MeHg2 

Concentration  of  total  MeHg  in  sediment 

ng  L-1 

MeHgdp2 

Concentration  of  dissolved  MeHg  in  pore  water 

ng  L-1 

MGHgdocp2 

Concentration  of  DOC  adsorbed  MeHg  in  pore  water 

ng  L1 

M6Hgpom2 

Concentration  of  POM  adsorbed  MeHg  in  sediment 

ng  L-1 

MeHgPt2 

Total  concentration  of  solids  adsorbed  MeHg  in  sediment 

ng  L-1 

MeHgPts2 

Total  concentration  of  solids  adsorbed  MeHg  in  sediment 

ngg1 

Mercury  pathway  fluxes  are  internally  computed  in  HgSM  and  can  be 
reported  in  the  model  outputs.  Table  27  lists  the  symbols  and  definitions 
of  pathway  fluxes  of  mercury  species  computed  in  HgSM.  The  HgSM 
allows  the  user  to  selectively  turn  on  and  off  each  variable  in  the  model 
outputs. 
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Table  27.  List  of  pathway  fluxes  of  mercury  species  computed  in  HgSM. 


Symbol 

Definition 

Units 

Elemental  mercury 

o 

a? 

A 

! 

V 

£ 

HgO  volatilization  in  water 

ng  L-i-d'1 

3? 

A 

6 

3? 

HgO  oxidation  into  Hgll  in  water 

ng  L1  d'1 

Inorganic  mercury 

Hgll~>HgO 

Hgll  photoreduction  into  HgO  in  water 

ng  L1  d'1 

Hgll->MeHg 

Hgll  methylation  into  MeHg  in  water 

ng  L'^-d'1 

Hgll~>Bed 

Hgll  settling 

ng  L1  d'1 

Bed->Hgll 

Hgll  re-suspension 

ng  L1  d'1 

Hgll<->Hgll2 

Hgll  sediment-water  transfer 

ng  L-i-d'1 

Hglh->MeHg2 

sediment  Hgll  methylation  into  MeHg 

ng  L'^-d'1 

Hgll~>Bed 

Hgll  deposition 

ng  L1  d'1 

Bed->HgU 

sediment  Hgll  erosion 

ng  L1  d'1 

Hgll<->Hgll2 

Hgll  sediment-water  transfer 

ng  L1  d'1 

Hglh  burial 

sediment  Hgll  burial 

ng  L'1  d'1 

Methymercury 

Atm<->MeHg 

MeHg  volatilization  in  water 

ng  L'^-d'1 

MeHg->HgO 

MeHg  photoreduction  into  HgO  in  water 

ng  L'1  d'1 

MeHg->Hgll 

MeHg  demethylation  into  Hgll  in  water 

ng  L'^-d'1 

MeHg->Bed 

MeHg  settling 

ng  L-^-d'1 

Bed->MeHg 

MeHg  re-  suspension 

ng  L-i-d'1 

MeHg<~>MeHg2 

MeHg  sediment-water  transfer 

ng  L'1  d'1 

MeHg2->Hgll2 

sediment  MeHg  bacterial  demethylation  into  Hgll 

ng  L-i-d'1 

MeHg->Bed 

MeHg  deposition 

ng  L'1  d'1 

Bed->MeHg 

sediment  MeHg  erosion 

ng  L-i-d'1 

MeHg<~>MeHg2 

MeHg  sediment-water  transfer 

ng  L'1  d'1 

MeHg2  burial 

sediment  MeHg  burial 

ng  L-i-d'1 
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5  Summary 

The  preceding  chapters  describe  three  newly  developed  GC,  CSM,  and 
HgSM  water  quality  modules.  These  modules  are  written  as  DLLs  and 
compiled  as  GC.dll,  CSM.dll,  and  HgSM.dll,  respectively.  Zhang  and 
Johnson  (2016)  present  general  principles  and  protocols  to  integrate  “plug 
in”  water  quality  modules  into  the  lD  HEC-RAS  model.  Following  the  same 
integration  approach  implemented  in  HEC-RAS  for  NSM  modules,  GC, 
CSM,  and  HgSM  modules  have  been  integrated  into  the  HEC-RAS  model. 
Wherein  the  HEC-RAS  is  used  to  describe  physical  processes  of  advection 
and  dispersion,  the  GC,  CSM,  and  HgSM  modules  are  used  to  compute 
speciation,  reactions  and  transformations  in  the  water  column  and 
sediment  layer.  HEC-RAS  obtains  kinetics  source  and  sink  terms  from  each 
water  quality  module  and  performs  the  transport  and  mass  balance 
computation  for  each  water  quality  cell  and  each  state  variable.  HEC-RAS 
has  the  control  over  the  initial  and  boundary  conditions.  The  water  quality 
module  returns  the  rate  change  of  state  variables,  derived  variables,  and 
pathways  back  to  HEC-RAS.  The  HEC-RAS  graphic  user  interface  with  pre- 
and  post-processer  can  be  utilized  to  setup  the  model,  perform  the  runs, 
and  present  and  analyze  the  results.  This  report  is  primarily  intended  for 
use  as  a  technical  reference  for  applying  these  water  quality  modules.  These 
three  water  quality  modules  will  also  be  integrated  into  other  H&H  models 
(e.g.  AdH  and  SRH-2D). 

This  report  provides  theory  and  mathematical  formulations  implemented 
in  GC,  CSM,  and  HgSM  modules.  The  report  describes  the  scientific  basis 
for  three  water  quality  modules  and  how  they  are  formulated.  The  GC  is  a 
supporting  water  quality  module  that  can  be  used  to  model  simple  kinetics 
for  solids  and  user-defined  constituents  in  the  water  column  and  bed 
sediments.  The  scientific  basis  of  the  CSM  and  HgSM  modules  reflects 
empirical  and  theoretical  support  precedence  in  the  literature,  and  in 
widely  used  toxic  chemical  models.  The  GC,  CSM,  and  HgSM  modules 
have  been  tested  and  verified  with  a  variety  of  examples,  and  compared 
against  existing  water  quality  models.  The  results  from  the  comparison 
showed  that  the  computer  codes  are  correctly  implemented  in  these  water 
quality  modules.  Additional  validation  of  these  modules  against  observed 
data  collected  in  aquatic  systems  is  underway.  Future  model  validation 
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will  evaluate  how  well  the  water  quality  module’s  theoretical  foundation 
and  computer  implementation  describe  actual  system  behavior  in  terms  of 
the  degree  of  correlation  between  computed  and  independently  observed 
responses  of  the  aquatic  system. 

The  computer  codes  in  all  of  the  water  quality  modules  have  been  checked 
for  consistency  with  the  formulations  described  in  the  report.  The  variable 
names  in  the  report  correspond  to  those  used  in  the  codes  so  that  the 
mathematical  formulations  and  codes  can  be  compared.  Through  a  series  of 
model  verification  tests,  the  model  formulations  appear  correct.  The 
modularity  forms  the  basis  for  the  flexibility  of  these  “plug  in”  water  quality 
modules  (e.g.  GC,  CSM,  HgSM),  including  the  ability  to  add  and  modify 
processes  for  future  development  and  enhancement.  Despite  the  CSM  and 
HgSM  modules  can  be  used  to  quantify  the  major  processes  of  contami¬ 
nants  and  mercury  species  in  aquatic  systems,  many  of  the  mechanisms  that 
govern  contaminant  partitioning  and  reactions  are  still  poorly  understood. 
Determing  contaminant  fate  in  the  environment  and  understanding  their 
potential  effects  on  human  and  environmental  health  is  a  complex  topic 
requiring  further  research.  These  water  quality  modules  can  be 
continuously  enhanced  and  further  developed  through  revision  of  water 
quality  formulations  and  computer  codes  presented  above  in  this  report. 
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Appendix  A:  Definition  of  Mathematical 
Symbols  used  in  the  GC  Module 


Symbol 

Definition 

Units 

& 

concentration  of  constituent  i 

mg  L1 

dpn 

diameter  of  solid  “n” 

mm 

d* 

dimensionless  solids  parameter 

unitless 

E 

sediment  erosion  rate 

g  cm-2  s1 

Eo 

surface  erosion  rate 

g  cm-2  s1 

g 

gravity  acceleration 

m  S"2 

h 

water  depth 

m 

h2 

sediment  layer  thickness 

m 

i 

index  number  of  user-defined  constituents 

unitless 

koi(T) 

zero-order  decay  rate  of  constituent  i 

mg  L^d'1 

kn(T) 

first-order  decay  rate  of  constituent  i 

d'1 

k(20) 

value  of  the  reaction  rate  at  20  °C 

d'1 

k(T) 

value  of  the  rate  at  local  temperature 

d'1 

k(Twr) 

value  of  the  rate  at  reference  temperature 

d-1 

m 

empirically  determined  constant 

unitless 

mn 

concentration  of  solid  “n"  in  water 

mg  L-1 

mn2 

concentration  of  solid  “n”  in  sediment 

mg  L-1 

n 

index  number  of  solids  classes 

unitless 

N 

number  of  solids  classes 

unitless 

P 

surface  erosion  rate 

g  cm-2  s1 

Pdpn 

deposition  probability  of  solid  “n” 

unitless 

Tw 

water  temperature 

°C 

u 

water  velocity 

m  s*1 

Vb 

sediment  burial  velocity 

m  yr1 

Vdpn 

deposition  velocity  of  solid  “n” 

m  d'1 

Vrpn 

re-suspension  velocity 

m  d-1 

Vs  i 

settling  velocity  of  constituent  i 

m  d-1 

Vs  pn 

settling  velocity  of  solid  “n” 

m  d-1 

Greek 

0 

sediment  layer  porosity 

unitless 

T 

bottom  shear  stress 

N  nr2 
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Symbol 

Definition 

Units 

Tee 

critical  shear  stress  for  erosion 

N  m-2 

Ten 

critical  shear  stress  for  the  initiation  of  movement  of 
noncohensive  solids 

N  m-2 

Tcdl 

lower  critical  shear  stress  for  deposition 

N  m-2 

Tedu 

upper  critical  shear  stress  for  deposition 

N  m-2 

Ppn 

density  of  solid  “n”  in  water 

g  cm'3 

Ppn2 

density  of  solid  "n”  in  sediment 

g  enr3 

Pw 

density  of  water 

g  cm-3 

Ps 

sediment  dry  density 

mg  L-1 

Pb 

sediment  bulk  density 

g  cm-3 

V 

kinematic  viscosity  of  water 

m2  s1 

ao 

empirical  constant 

unitless 

6 

temperature  correction  coefficient 

unitless 
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Appendix  B:  Definition  of  Mathematical 
Symbols  used  in  the  CSM  Module 


Symbol 

Definition 

Units 

Apd 

algal  biomass  in  water 

mg  L1 

As 

surface  water  area 

nr2 

b 

Freundlich  exponent 

unitless 

bpom 

Freundlich  exponent  for  POM  in  water 

unitless 

bpom2 

Freundlich  exponent  for  sediment  POM 

unitless 

bap 

Freundlich  exponent  for  algae  in  water 

unitless 

bpn 

Freundlich  exponent  for  solid  “n”  in  water 

unitless 

bpn2 

Freundlich  exponent  for  sediment  solid  “n” 

unitless 

Ct 

total  oncentration  in  water 

pgL1 

Cion 

ionic  concentration  in  water 

Mg  L1 

cd 

concentration  of  dissolved  phase  in  water 

pgL1 

Cap 

concentration  of  algae  adsorbed  phase  in  water 

pgL1 

Cdoc 

concentration  of  DOC  adsorbed  phase  in  water 

Cpom 

concentration  of  POM  adsorbed  phase  in  water 

pgL1 

Cpn 

concentration  of  solids  “n”  adsorbed  phase  in  water 

pgL1 

Cpts 

total  mass  of  solids  adsorbed  phase  in  water 

mg  kg1 

Cps 

mass  of  solids  “n”  adsorbed  phase  in  water 

mg  kg1 

Cpt 

total  concentration  of  solids  adsorbed  phase  in  water 

pgL1 

Ct2 

total  concentration  in  sediment 

pgL1 

Cd  2 

concentration  of  dissolved  phase  in  sediment 

pgL1 

Cdp2 

concentration  of  dissolved  phase  in  pore  water 

pgL1 

Cdoc2 

concentration  of  DOC  adsorbed  phase  in  sediment 

pgL1 

Cdocp2 

concentration  of  DOC  adsorbed  phase  in  pore  water 

pgL1 

Cpom2 

concentration  of  POM  adsorbed  phase  in  sediment 

pgL1 

Cpn2 

concentration  of  solid  “n”  adsorbed  phase  in  sediment 

pgL1 

Cpts2 

total  mass  of  solids  adsorbed  phase  in  sediment 

mg  kg1 

Cpt2 

total  concentration  of  solids  adsorbed  phase  in  sediment 

Co 

air  concentration  (gaseous)  of  contaminant 

pgL1 

Csd 

solubility 

mg  L1 

DOC 

dissolved  organic  carbon  in  water 

mg-C  L1 
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Symbol 

Definition 

Units 

D0C2 

sediment  dissolved  organic  carbon  in  pore  water 

mg-C  L1 

Db 

biodiffusion  coefficient  representing  particle  diffusivity 

m2  d1 

Dm 

molecular  diffusion  coefficient 

m2  s1 

Ea 

activation  energy 

J  mol'1 

[E] 

concentration  of  environmental  property  driving  generalized 
second-order  reaction  in  water 

mg  L1 

[E2] 

concentration  of  environmental  property  driving  generalized 
second-order  reaction  in  sediment 

mg  L1 

f 

fraction  of  the  chemical  occurring  in  each  of  the  chemical 
species  i 

unitless 

fap 

Freundlich  exponent  for  algae  in  water 

unitless 

fpom 

Freundlich  exponent  for  POM  in  water 

unitless 

fpom2 

Freundlich  exponent  for  POM  in  sediment 

unitless 

fpn 

Freundlich  exponent  for  solid  “n”  in  water 

unitless 

fpn2 

Freundlich  exponent  for  solid  “n”  in  sediment 

unitless 

fa 

fraction  of  dissolved  phase  in  water 

unitless 

fap 

fraction  of  algae  adsorbed  phase  in  water 

unitless 

fdoc 

fraction  of  DOC  adsorbed  phase  in  water 

unitless 

fpom 

fraction  of  POM  adsorbed  phase  in  water 

unitless 

fraction  of  solid  “n”  adsorbed  phase  in  water 

unitless 

fd2 

fraction  of  dissolved  phase  in  sediment 

unitless 

fdoc2 

fraction  of  DOC  adsorbed  phase  in  sediment 

unitless 

fpom2 

fraction  of  POM  adsorbed  phase  in  sediment 

unitless 

foe 

fraction  of  organic  carbon  in  solids 

unitless 

fpn2 

fraction  of  solid  “n”  adsorbed  phase  in  sediment 

unitless 

h 

water  depth 

m 

hi 

sediment  layer  thickness 

m 

[H+] 

hydronium  ions 

mol  L1 

[OH] 

hydroxide  ions 

mol  L1 

lo 

solar  radiation  at  the  water  surface 

Wm2 

lopht 

light  intensity  when  kPbt(T)  is  measured 

Wm2 

j 

species  that  can  be  transformed  into  daughter  products 

unitless 

K, 

equilibrium  constant  for  the  formation  of  the  acid  (Kai),  or 
anionic  species,  or  the  base  (Kbi)  or  cationic  species 

unitless 

Kd 

sorption  or  partition  coefficient 

L  kg 1 

Kf 

Freundlich  adsorption  constant 

(Mg  g1)  (Mg  L'1)'13 

Kfap 

Freundlich  adsorption  constant  for  algae  in  water 

(Mg  g1)  (Mg  L-1)* 

Kfpom 

Freundlich  adsorption  constant  for  POM  in  water 

(Mg  g1)  (Mg  L'1)'6 
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Symbol 

Definition 

Units 

Kfpom2 

Freundlich  adsorption  constant  for  sediment  POM 

(Mg  g-1)  (Mg  L-1)* 

Kfpn 

Freundlich  adsorption  constant  for  solid  “n”  in  water 

(Mg  g1)  (Mg  L  1)'b 

Kfpn2 

Freundlich  adsorption  constant  for  sediment  solid  “n” 

(Mg  g1)  (Mg  L1)15 

K, 

Langmuir  adsorption  constant 

L  Mg 1 

K/ap 

Langmuir  adsorption  constant  for  algae  in  water 

L  Mg 1 

Klpom 

Langmuir  adsorption  constant  for  POM  in  water 

L  Mg 1 

Klpom2 

Langmuir  adsorption  constant  for  sediment  POM 

L  Mg 1 

Klpn 

Langmuir  adsorption  constant  for  solid  “n”  in  water 

L  Mg 1 

Klpn2 

Langmuir  adsorption  constant  for  sediment  solid  “n” 

L  Mg 1 

Kdoc 

equilibrium  partition  coefficient  for  DOC  in  water 

L  kg1 

Kap 

equilibrium  partition  coefficient  for  algae  in  water 

L  kg 1 

Kp 

equilibrium  partition  coefficient  for  solid  “n”  in  water 

L  kg 1 

Kpom 

equilibrium  partition  coefficient  for  POM  in  water 

L  kg 1 

Kpn 

equilibrium  partition  coefficient  for  solid  “n”  in  water 

L  kg 1 

Kdoc2 

equilibrium  partition  coefficient  for  sediment  DOC 

L  kg 1 

Kpom2 

equilibrium  partition  coefficient  for  sediment  POM 

L  kg 1 

Kpn2 

equilibrium  partition  coefficient  for  sediment  solid  “n” 

L  kg 1 

Kow 

octanol-water  partition  coefficient 

unitless 

kadn 

adsorption  coefficient  for  solid  “n”  in  water 

L  pg 1  d1 

kadn2 

adsorption  coefficient  for  sediment  solid  "n” 

L  pg 1  d1 

kadap 

adsorption  coefficient  for  algae  in  water 

L  pg 1  d1 

kdaap 

desorption  rate  for  algae  in  water 

d-1 

kadpom 

adsorption  coefficient  for  POM  in  water 

L  pg  i  d1 

kdapom 

desorption  rate  for  POM  in  water 

d-1 

kadpom2 

adsorption  coefficient  for  sediment  POM 

L  pg  i  d1 

kdapom2 

desorption  rate  for  sediment  POM 

d-1 

kdan 

desorption  rate  for  solid  “n”  in  water 

d-1 

kdan2 

desorption  rate  for  solid  "n”  in  sediment 

d-1 

KAt) 

degradation  rate  for  dissolved  species  in  water 

d-1 

degradation  rate  for  dissolved  species  in  sediment 

d-1 

kip(T) 

degradation  rate  for  algae  adsorbed  species  in  water 

d-1 

Kdoc(T) 

degradation  rate  for  DOC  adsorbed  species  in  water 

d-1 

Kdocl  00 

degradation  rate  for  DOC  adsorbed  species  in  sediment 

d-1 

Kpom(T) 

degradation  rate  for  POM  adsorbed  species  in  water 

d-1 

degradation  rate  for  POM  adsorbed  species  in  sediment 

d-1 
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Symbol 

Definition 

Units 

Kpn  (2”) 

degradation  rate  for  solid  “n”  adsorbed  species  in  water 

d-1 

pn  2  CO 

degradation  rate  for  solid  “n”  adsorbed  species  in  sediment 

d-1 

Ka<J) 

acid  hydrolysis  rate  constant  for  dissolved  species  in  water 

L  mol1  d'1 

Ka^T) 

acid  hydrolysis  rate  constant  for  dissolved  species  in 
sediment 

L  mol-1  d'1 

KJT) 

base  hydrolysis  rate  constant  for  dissolved  species  in  water 

L  mol1  d'1 

Ka,^T) 

base  hydrolysis  rate  constant  for  dissolved  species  in 
sediment 

L  mol1  d'1 

Kn  (T) 

neutral  hydrolysis  rate  for  dissolved  species  in  water 

d-1 

kLi(T) 

neutral  hydrolysis  rate  for  dissolved  species  in  sediment 

d-1 

khadoc  CO 

acid  hydrolysis  rate  constant  for  DOC  adsorbed  species  in 
water 

L  mol1  d'1 

khadoc2  (O 

acid  hydrolysis  rate  constant  for  DOC  adsorbed  species  in 
sediment 

L  mol1  d'1 

^hbdoc  CO 

base  hydrolysis  rate  constant  for  DOC  adsorbed  species  in 
water 

L  mol-1  d'1 

k hbdocl  (r) 

base  hydrolysis  rate  constant  for  DOC  adsorbed  species  in 
sediment 

L  mol1  d'1 

k hndoc  CO 

neutral  hydrolysis  rate  for  DOC  adsorbed  chemical  species  in 
water 

d-1 

k hndocl  CO 

neutral  hydrolysis  rate  for  DOC  adsorbed  chemical  species  in 
sediment 

d-1 

kpht 

direct  near-surface  photolysis  rate  for  dissolved  species 

d-1 

kphtdoc 

direct  near-surface  photolysis  rate  for  DOC  adsorbed  species 

d-1 

Kl 

mass  transfer  velocity  from  the  liquid  film 

m  d'1 

Kg 

mass  transfer  velocity  from  the  gaseous  film 

m  d'1 

Kh 

Henry’s  Law  constant 

Pa  m3  moh1 

ka 

oxygen  reaeration  rate 

m  d-1 

ked 

second  order  rate  for  dissolved  species  in  water 

(mg  L1)-1  d-1 

ked2 

second  order  rate  for  dissolved  species  in  sediment 

(mg  L1)-1  d'1 

kedoc 

second  order  rate  for  DOC  adsorbed  species  in  water 

(mg  L1)-1  d'1 

kedoc2 

second  order  rate  for  DOC  adsorbed  species  in  sediment 

(mg  L1)-1  d-1 

keap 

second  order  rate  for  algae  adsorbed  species 

(mg  L1)-1  d-1 

kepom 

second  order  rate  for  POM  adsorbed  species  in  water 

(mg  L1)-1  d-1 

kepom2 

second  order  rate  for  POM  adsorbed  species  in  sediment 

(mg  L1)-1  d'1 

kepn 

second  order  rate  for  solids  adsorbed  species  in  water 

(mg  L1)-1  d'1 

kepn2 

second  order  rate  for  solids  adsorbed  species  in  sediment 

(mg  L1)-1  d'1 

kj-»(T) 

transformation  rate 

d-1 
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Symbol 

Definition 

Units 

Kal 

ionization  constant  for  anionic  species  in  water 

unitless 

Kal2 

ionization  constant  for  anionic  species  in  sediment 

unitless 

Ka2 

ionization  constant  for  anionic  species  in  water 

unitless 

Ka2  2 

ionization  constant  for  anionic  species  in  sediment 

unitless 

Kbl 

ionization  constant  for  cationic  species  in  water 

unitless 

Kbl2 

ionization  constant  for  cationic  species  in  sediment 

unitless 

Kb2 

ionization  constant  for  cationic  species  in  water 

unitless 

Kb  22 

ionization  constant  for  cationic  species  in  sediment 

unitless 

Lo 

areal  deposition  of  contaminant 

g  nr2  d'1 

mn 

concentration  of  solid  “n”  in  water 

mg  L1 

rrin2 

concentration  of  solid  “n”  in  sediment 

mg  L1 

MW 

molecular  weight 

g  M1 

n 

solids  index 

unitless 

POM 

particulate  organic  matter  in  water 

mg  L1 

POM2 

sediment  particulate  organic  matter 

mg  L1 

pH 

measure  of  the  acidity  or  basicity 

unitless 

Qc 

maximum  amount  adsorbed  by  the  solid 

mg  kg-1 

Cjcap 

maximum  amount  of  contaminant  adsorbed  by  algae  in  water 

mg  kg1 

Qcpom 

maximum  amount  of  contaminant  adsorbed  by  POM  in  water 

mg  kg1 

Qcpom2 

maximum  amount  of  contaminant  adsorbed  by  sediment 

POM 

mg  kg1 

qcr 

maximum  amount  of  contaminant  adsorbed  by  solid  "n”  in 
water 

mg  kg1 

Qcn2 

maximum  amount  of  contaminant  adsorbed  by  sediment 
solid  "n” 

mg  kg1 

R 

universal  gas  constant 

Pa  m3  mol'1  °K1 

Sdec 

contaminant  degradation 

pg  L1 d1 

Shyd 

contaminant  hydrolysis 

Mg  L1 d1 

Spht 

contaminant  photolysis 

pg  L1 d1 

Sexa 

contaminant  generalized  second-order  reaction 

pg  L1 d1 

Strm 

total  contaminant  transformation 

pg  L1 d1 

Svlt 

contaminant  volatilization 

pg  L1 d1 

Satm 

contaminant  air  deposition 

pg  L1 d1 

Sset 

contaminant  settling 

pg  L1 d1 

Sres 

contaminant  re-suspension 

pg  L1 d1 

SSt 

sum  of  internal  source  and  sink  terms  of  total  contaminant  in 
water 

pg  L1 d1 
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Symbol 

Definition 

Units 

SSd 

sum  of  internal  source  and  sink  terms  of  dissolved  phase  in 
water 

pg  L1  d-i 

SSdoc 

sum  of  internal  source  and  sink  terms  of  DOC  adsorbed 
phase  in  water 

pg  L-i  d-1 

SSap 

sum  of  internal  source  and  sink  terms  of  algae  adsorbed 
phase  in  water 

pg  L1  d-i 

SSpom 

sum  of  internal  source  and  sink  terms  of  POM  adsorbed 
phase  in  water 

Mg  L1 d1 

SSpn 

sum  of  internal  source  and  sink  terms  of  solids  adsorbed 
phase  in  water 

Mg  L1 d1 

SSt2 

sum  of  source  and  sink  terms  of  total  contaminant  in 
sediment 

Mg  L1 d1 

SSd2 

sum  of  source  and  sink  terms  of  dissolved  phase  in  sediment 

Mg  L1 d1 

SSdoc2 

sum  of  source  and  sink  terms  of  DOC  adsorbed  phase  in 
sediment 

Mg  L1 d1 

SSpom2 

sum  of  source  and  sink  terms  of  POM  adsorbed  phase  in 
sediment 

Mg  L1 d1 

SSpn2 

sum  of  source  and  sink  terms  of  solids  adsorbed  phase  in 
sediment 

Mg  L1 d1 

sb 

sediment  contaminant  burial 

Mg  L1 d1 

Ssw 

contaminant  sediment-water  transfer 

Mg  L1 d1 

Twk 

water  temperature 

°K 

Tw 

water  temperature 

°C 

Twr 

reference  water  temperature 

°K 

Uw 

wind  speed 

m  s1 

u* 

flow  shear  velocity  along  the  bed,  which  is  approximately  10 
percent  of  the  mean  velocity  of  flow 

m  S"1 

Vb 

sediment  burial  velocity 

m  d-1 

Vm 

sediment-water  transfer  velocity 

m  d-1 

Vms 

user-defined  sediment-water  transfer  velocity 

m  d-1 

Vs  pm 

solids  settling  velocity 

m  d'1 

Vrpn 

re-suspension  velocity 

m  d'1 

Vs  ap 

algal  settling  velocity 

m  d'1 

Vs  0/77 

POM  settling  velocity 

m  d-1 

Vv(T) 

volatilization  velocity  across  the  air-water  interface 

m  d-1 

Yj->i 

transformation  yield  coefficient  from  species  j  to  i 

g-i/g-j 

Zl 

pore  water  diffusion  layer  thickness 

mm 

Z2 

sediment  bioturbation  layer  thickness 

mm 

Greek 

Of  doc 

partition  correlation  coefficient  with  DOC 

unitless 

ERDC/EL  TR-16-8 


119 


Symbol 

Definition 

Units 

Ofap 

partition  correlation  coefficient  with  algae 

unitless 

Ofpom 

partition  correlation  coefficient  with  POM 

unitless 

Ofpn 

partition  correlation  coefficient  with  solids 

unitless 

Amax 

maximium  light  extinction  coefficient 

nr1 

ai 

correction  factor  for  light  attenuation 

unitless 

A 

light  extinction  coefficient 

nrv1 

pb 

sediment  bulk  density 

g  cm-3 

e 

temperature  correction  coefficient 

unitless 

<P 

sediment  layer  porosity 

unitless 

V 

kinematic  viscosity  of  water 

m^d-1 
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Appendix  C:  Definition  of  Mathematical 
Symbols  used  in  the  HgSM  Module 


Symbol 

Definition 

Units 

As 

surface  water  area 

nr2 

Atm<->HgO 

HgO  volatilization  in  water 

ng  L1  d'1 

Atm->HgO 

HgO  air  deposition 

ng  L1  d'1 

Atm<->MeHg 

MeHg  volatilization  in  water 

ng  L1  d'1 

Atm~>MeHg 

MeHg  air  deposition 

ng  L1  d'1 

bpom 

Freundlich  exponent  for  POM  in  water 

unitless 

bpom2 

Freundlich  exponent  for  sediment  POM 

unitless 

bap 

Freundlich  exponent  for  algae  in  water 

unitless 

bpn 

Freundlich  exponent  for  solid  “n”  in  water 

unitless 

bpn2 

Freundlich  exponent  for  sediment  solid  "n” 

unitless 

Bed->Hgll 

sediment  Hgl  1  re-suspension  rate 

ng  L1  d'1 

Bed->MeHg 

sediment  MePIg  re-suspension  rate 

ng  L1  d'1 

DOC 

dissolved  organic  carbon  in  water 

mg-C  L1 

D0C2 

dissolved  organic  carbon  in  sediment  pore  water 

mg-C  L1 

fap 

Freundlich  exponent  for  algae  in  water 

unitless 

fd-Hgll 

fraction  of  dissolved  Hgl  1  in  water 

unitless 

fdoc-Hgll 

fraction  of  DOC  adsorbed  Hgl  1  in  water 

unitless 

fap-Hgll 

fraction  of  algae  adsorbed  Hgll  in  water 

unitless 

fpom-Hgll 

fraction  of  POM  adsorbed  Hgl  1  in  water 

unitless 

fp-Hgll 

fraction  of  solids  adsorbed  Hgll  in  water 

unitless 

fd-MeHg 

fraction  of  dissolved  MePIg  in  water 

unitless 

fdoc-MeHg 

fraction  of  DOC  adsorbed  MePIg  in  water 

unitless 

fap-MeHg 

fraction  of  algae  adsorbed  MePIg  in  water 

unitless 

fpom-MeHg 

fraction  of  POM  adsorbed  MePIg  in  water 

unitless 

fp-MeHg 

fraction  of  solids  adsorbed  MePIg  in  water 

unitless 

fd-Hgll 2 

fraction  of  dissolved  Hgl  1  in  sediment 

unitless 

fdoc-Hgll2 

fraction  of  DOC  adsorbed  Hgl  1  in  sediment 

unitless 

fpom-Hgll2 

fraction  of  POM  adsorbed  Hgl  1  in  sediment 

unitless 

fp-Hgll2 

fraction  of  solids  adsorbed  Hgll  in  sediment 

unitless 

fd-MeHg2 

fraction  of  dissolved  MePIg  in  sediment 

unitless 
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Symbol 

Definition 

Units 

fdoc-MeHg2 

fraction  of  DOC  adsorbed  MeHg  in  sediment 

unitless 

fpom-MeHg2 

fraction  of  POM  adsorbed  MeHg  in  sediment 

unitless 

fp-MeHg 

fraction  of  solidsad  sorbed  MeHg  in  sediment 

unitless 

fpom 

Freundlich  exponent  for  POM  in  water 

unitless 

fpom2 

Freundlich  exponent  for  sediment  POM 

unitless 

fpn 

Freundlich  exponent  for  solid  “n”  in  water 

unitless 

fpn2 

Freundlich  exponent  for  sediment  solid  “n” 

unitless 

HgOo 

air  concentration  of  HgO  (gaseous) 

ngL1 

HgO->Hgll 

HgO  oxidation  yield  into  Hgll  in  water 

ng  L1  d'1 

o 

I? 

A 

3? 

Hgll  photoreduction  yield  into  HgO  in  water 

ng  L1  d'1 

Hgll->MeHg 

Hgll  methylation  yield  in  water 

ng  L1  d'1 

Hglh->MeHg2 

Hgll  methylation  yield  in  sediment 

ng  L1  d'1 

Hgll~>MeHg 

Hgll  methylation  yield  in  water 

ng  L1  d'1 

Hgll~>Bed 

Hgll  settling  in  water 

ng  L1  d'1 

Hgll 2  burial 

sediment  Hgll  burial  in  water 

ng  L1  d'1 

Hgll<->Hgll2 

Hgll  sediment-water  transfer 

ng  L1  d'1 

HgO 

concentration  of  elemental  mercury  in  water 

ng  L1 

Hgll 

concentration  of  total  Hgll  in  water 

ng  L1 

Hglld 

concentration  of  dissolved  Hgll  in  water 

ng  L1 

Hgllpn 

concentration  of  solids  “n”  adsorbed  Hgll  in  water 

ng  L1 

Hgllp 

total  concentration  of  solids  adsorbed  Hgll  in  water 

ng  L1 

Hgllpts 

total  concentration  of  solids  adsorbed  Hgll  in  water 

ng  kg1 

Hgll  2 

concentration  of  total  Hgll  in  sediment 

ng  L1 

Hglld2 

concentration  of  dissolved  Hgll  in  sediment 

ngL1 

HgUdp2 

concentration  of  dissolved  Hgll  in  pore  water 

ng  L1 

Hgll  doc 

concentration  of  DOC  adsorbed  Hgll  in  water 

ngL1 

Hglldocp2 

concentration  of  DOC  adsorbed  Hgll  in  pore  water 

ngL1 

Hgllp2 

total  concentration  of  solids  adsorbed  Hgll  in  sediment 

ngL1 

Hgllpts2 

total  concentration  of  solids  adsorbed  Hgll  in  sediment 

ng  kg1 

h 

water  depth 

m 

h2 

sediment  layer  thickness 

m 

lo 

solar  radiation  at  the  water  surface 

W  nr2 

lopht 

light  intensity  at  which  kPnt(T)  is  measured 

W  nr2 

Kdoc-MeHg 

MeHg  partitioning  coefficient  for  DOC  in  water 

L  kg 1 

KpMeHg 

MeHg  partitioning  coefficient  for  solids  in  water 

L  kg 1 

Kap-  MeHg 

MeHg  Partitioning  coefficient  for  algae  in  water 

L  kg 1 
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Symbol 

Definition 

Units 

Kpom-MeHg 

MeHg  partitioning  coefficient  for  POM  in  water 

L  kg 1 

Kdoc-Hgll 

Hgll  equilibrium  partitioning  coefficient  for  DOC  in  water 

L  kg 1 

Kap-Hgll 

Hgll  equilibrium  partitioning  coefficient  for  algae  in  water 

L  kg 1 

Kpom-Hgll 

Hgll  equilibrium  partitioning  coefficient  for  POM  in  water 

L  kg 1 

Kdoc-Hgll2 

Hgll  equilibrium  partitioning  coefficient  for  sediment  DOC 

L  kg 1 

Kpom-Hgll2 

Hgll  equilibrium  partitioning  coefficient  for  sediment  POM 

L  kg 1 

Kp-Hgll2 

Hgll  equilibrium  partitioning  coefficient  for  sediment  solids 

L  kg 1 

Kdoc-MeHg2 

MeHg  partitioning  coefficient  for  sediment  DOC 

L  kg 1 

Kpom-MeHg2 

MeHg  Partitioning  coefficient  for  sediment  POM 

L  kg 1 

Kp-MeHg2 

MeHg  partitioning  coefficient  for  sediment  solids 

L  kg 1 

Kfpom-Hgll 

Hgll  Freundlich  adsorption  constant  for  POM  in  water 

(Mg  g1)  (Mg  L'1)'13 

Kfpom-Hgll2 

Hgll  Freundlich  adsorption  constant  for  sediment  POM 

(Mg  g1)  (Mg  L-1)-13 

Kfap-Hgll 

Hgll  Freundlich  adsorption  constant  for  algae  in  water 

(Mg  g1)  (Mg  L-1)-13 

Kfp-Hgll 

Hgll  Freundlich  adsorption  constant  for  solid  “n”  in  water 

(Mg  g1)  (Mg  L-1)-13 

Kfp-Hgll2 

Hgll  Freundlich  adsorption  constant  for  sediment  solid  “n” 

(Mg  g1)  (Mg  L-1)-13 

Kfpom-MeHg 

MeHg  Freundlich  adsorption  constant  for  POM  in  water 

(Mg  g'1)  (Mg  L_1)'b 

Kfpom-MeHg2 

MeHg  Freundlich  adsorption  constant  for  sediment  POM 

(Mg  g'1)  (Mg  L_1)'b 

Kfap-MeHg 

MeHg  Freundlich  adsorption  constant  for  algae  in  water 

(Mg  g 1)  (Mg  L'1)-*3 

Kfp-MeHg 

MeHg  Freundlich  adsorption  constant  for  solid  “n”  in  water 

(Mg  g1)  (Mg  L-1)-13 

Kfp-MeHg2 

MeHg  Freundlich  adsorption  constant  for  sediment  solid  “n” 

(Mg  g1)  (Mg  L_1)"b 

Klpom-Hgll 

Hgll  Langmuir  adsorption  constant  for  POM  in  water 

L  Mg 1 

Klpom-Hgil2 

Hgll  Langmuir  adsorption  constant  for  sediment  POM 

L  Mg1 

Klap-Hgll 

Hgll  Langmuir  adsorption  constant  for  algae  in  water 

L  Mg1 

Klp-Hgll 

Hgll  Langmuir  adsorption  constant  for  solid  “n”  in  water 

L  Mg 1 

Klp-Hgll2 

Hgll  Langmuir  adsorption  constant  for  sediment  solid  "n” 

L  Mg 1 

Klpom-MeHg 

MeHg  Langmuir  adsorption  constant  for  POM  in  water 

L  Mg 1 

Klpom-MeHg2 

MeHg  Langmuir  adsorption  constant  for  sediment  POM 

L  Mg 1 

Klap-MeHg 

MeHg  Langmuir  adsorption  constant  for  algae  in  water 

L  Mg 1 

Klp-MeHg 

MeHg  Langmuir  adsorption  constant  for  solid  “n”  in  water 

L  Mg 1 

Klp-MeHg2 

MeHg  Langmuir  adsorption  constant  for  sediment  solid  “n” 

L  Mg 1 

Khl-HgO 

Henry’s  Law  constant  of  HgO 

Pa  m3  mol1 

Kw-MeHg 

Henry’s  Law  constant  of  MeHg 

Pa  m3  mol1 

Ks04 

half-saturation  constant  for  the  effect  of  sulfate  on 
methylation 

mg-02  L'1 

kso42(T) 

sediment  S04  reduction  rate 

d-1 

kdm32 

sediment  MeHg  demethylation  rate 

d-1 

kadn 

Hgll  adsorption  coefficient  for  solid  “n”  in  water 

L  pg1  d  1 
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Symbol 

Definition 

Units 

kdan 

Hgll  desorption  rate  for  solid  “n”  in  water 

d-1 

kadap 

Hgll  adsorption  coefficient  for  algae  in  water 

L  pg 1  d1 

kadpom 

Hgll  adsorption  coefficient  for  POM  in  water 

L  pg 1  d1 

kdaap 

Hgll  desorption  rate  for  algae  in  water 

d-1 

kdapom 

Hgll  desorption  rate  for  POM  in  water 

d-1 

kadn2 

Hgll  adsorption  coefficient  for  solid  "n”  in  sediment 

L  pg 1  d1 

kadpom2 

Hgll  adsorption  coefficient  for  sediment  POM 

L  pg-1  d1 

kdan2 

Hgll  desorption  rate  for  solid  "n”  in  sediment 

d-1 

kdapom2 

Hgll  desorption  rate  for  sediment  POM 

d-1 

ki2(T) 

HgO  oxidation  yield  rate  in  water 

d-1 

kd21 

Hgll  dissolved  Hgll  photoreduction  rate  in  water 

d-1 

kdoc21 

DOC  adsorbed  Hgll  photoreduction  rate  in  water 

d-1 

kd23(T) 

dissolved  Hgll  methylation  rate  in  water 

d-1 

kd23-2(T) 

dissolved  Hgll  methylation  rate  in  sediment 

d-1 

kdoc23(T) 

DOC  adsorbed  Hgll  methylation  rate  in  water 

d-1 

kd31 

dissolved  MeHg  photoreduction  rate  in  water 

d-1 

kdoc31 

DOC  adsorbed  MeHg  photoreduction  rate  in  water 

d-1 

kd32 

dissolved  MeHg  demethylation  rate  in  water 

d-1 

kdoc32 

DOC  adsorbed  MeHg  demethylation  rate  in  water 

d-1 

kd32-2(T) 

dissolved  MeHg  demethylation  rate  in  sediment 

d-1 

Lngii 

Hgll  areal  air  deposition  rate 

mg  nr2  d_1 

1-MeHg 

MeHg  areal  air  deposition  rate 

mg  nr2  d'1 

MeHg->HgO 

MeHg  photodegradation  yield  into  HgO  in  water 

ng  L'1  d'1 

MeHg~>Hgll 

MeHg  demethylation  yield  into  Hgll  in  water 

ng  L'1  d'1 

MeHg2~>Hgll2 

MeHg  demethylation  yield  into  Hgll  in  sediment 

ng  L'1  d'1 

MeHg~>Bed 

MeHg  settling  in  water 

ng  L'1  d'1 

MeHg2  burial 

sedimate  MeHg  burial  in  water 

ng  L'1  d'1 

MeHg<->MeHg2 

MeHg  sediment-water  transfer 

ng  L1 d1 

MeHgo 

air  concentration  of  MeHg 

ng  L1 

MeHg 

concentration  of  total  MeHg  in  water 

ng  L1 

MeHgd 

concentration  of  dissolved  MeHg  in  water 

ng  L1 

MGH^doc 

concentration  of  DOC  adsorbed  MeHg  in  water 

ng  L1 

MeHgp 

total  concentration  of  solids  adsorbed  MeHg  in  water 

ng  L1 

MeHgpts 

total  concentration  of  solids  adsorbed  MeHg  in  water 

ng  kg1 

MeHg2 

concentration  of  total  MeHg  on  in  sediment 

ng  L1 

MeHgdp2 

concentration  of  dissolved  MeHg  in  pore  water 

ng  L1 
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Symbol 

Definition 

Units 

MoH^docp2 

concentration  of  DOC  adsorbed  MeHg  in  pore  water 

ng  L1 

MeHgP2 

total  concentration  of  solids  adsorbed  MeHg  in  sediment 

ng  L1 

MeHgPts2 

total  concentration  of  solids  adsorbed  MeHg  in  sediment 

ng  kg1 

Qcn-Hgll 

maximum  amount  of  Hgll  adsorbed  by  solid  “n”  in  water 

Pgg1 

Qcn-Hgll2 

maximum  amount  of  Hgll  adsorbed  by  sediment  solid  “n” 

Mgg1 

Qcap-Hgll 

maximum  amount  of  Hgll  adsorbed  by  algae  in  water 

Pgg-1 

Cfcpom-Hgll 

maximum  amount  of  Hgll  adsorbed  by  POM  in  water 

Mgg1 

Qcpom-Hgll2 

maximum  amount  of  Hgll  adsorbed  by  sediment  POM 

Mgg1 

Qcn-Hgll 

maximum  amount  of  Hgll  adsorbed  by  suspended  solids 

Pgg1 

Qcn-Hgll2 

maximum  amount  of  Hgll  adsorbed  by  sediment  solids 

Pgg-1 

Qcap-MeHg 

maximum  amount  of  MeHg  adsorbed  by  algae  in  water 

Pgg1 

Qcpom-MeHg 

maximum  amount  of  MeHg  adsorbed  by  POM  in  water 

Mgg1 

Qcpom-MeHg2 

maximum  amount  of  MeHg  adsorbed  by  sediment  POM 

Mgg1 

Qcn-MeHg 

maximum  amount  of  MeHg  adsorbed  by  solid  “n”  in  water 

Pgg1 

Qcn-MeHg2 

maximum  amount  of  MeHg  adsorbed  by  sediment  solids 

Pgg1 

r mso4 

ratio  of  sediment  methylation  rate  and  sulfate  reduction  rate 

L  mg1 

S042 

sediment  pore  water  sulfate 

mg-02  L1 

SSwgo 

sum  of  internal  source  and  sink  terms  of  HgO  in  water 

ng  L1  d'1 

SShsii 

sum  of  internal  source  and  sink  terms  of  total  Hgll  in  water 

ng  L1  d'1 

SSHglld 

sum  of  internal  source  and  sink  terms  of  dissolved  Hgll  in 
water 

ng  L1  d'1 

SSHglldoc 

sum  of  internal  source  and  sink  terms  of  DOC  adsorbed  Hgll 
in  water 

ng  L1  d'1 

SSHgllap 

sum  of  internal  source  and  sink  terms  of  algae  adsorbed 

Hgll  in  water 

ng  L1  d'1 

SSHgllpom 

sum  of  internal  source  and  sink  terms  of  POM  adsorbed  Hgll 
in  water 

ng  L1  d'1 

SSHgllp 

sum  of  internal  source  and  sink  terms  of  solids  adsorbed 

Hgll  in  water 

ng  L1  d'1 

SSHgll2 

sum  of  source  and  sink  terms  of  total  Hgll  in  sediment 

ng  L1  d'1 

SSHglld2 

sum  of  source  and  sink  terms  of  dissolved  Hgll  in  sediment 

ng  L1  d'1 

SSHglldoc2 

sum  of  source  and  sink  terms  of  DOC  adsorbed  Hgll  in 
sediment 

ng  L'1  d'1 

SSHgllpom2 

sum  of  source  and  sink  terms  of  POM  adsorbed  Hgll  in 
sediment 

ng  L'1  d'1 

SSng//p2 

sum  of  source  and  sink  terms  of  solids  adsorbed  Hgll  in 
sediment 

ng  L'1  d'1 

sum  of  internal  source  and  sink  terms  of  total  MeHg  in  water 

ng  L'1  d'1 

sum  of  source  and  sink  terms  of  total  MeHg  in  sediment 

ng  L'1  d'1 
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Symbol 

Definition 

Units 

Vb 

sediment  burial  velocity 

m  d-1 

Vm 

sediment-water  transfer  velocity 

m  d'1 

Vrpn 

re-suspension  velocity  for  solid  “n” 

m  d'1 

Vspn 

settling  velocity  for  solid  “n” 

m  d'1 

Vsap 

algal  settling  velocity 

m  d-1 

Vsom 

POM  settling  velocity 

m  d-1 

Vv-Hgo(T) 

volatilization  velocity  of  HgO 

m  d-1 

Vv-MeHg(T) 

volatilization  velocity  of  MeHg 

m  d'1 

Yl2 

HgO  oxidation  yield  coefficient  in  water 

unitless 

Yll 

Hgll  reduction  yield  coefficient  in  water 

unitless 

Y23 

Hgll  methylation  yield  coefficient  in  water 

unitless 

Y31 

MeHg  photodegradation  yield  coefficient  in  water 

unitless 

Y32 

MeHg  demethylation  yield  coefficient  in  water 

unitless 

Greek 

Amax 

maximum  light  extinction  coefficient 

nrr1 

e 

temperature  correction  coefficient 

unitless 

<P 

sediment  layer  porosity 

unitless 
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